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→CAAthrAfused 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

ACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTG  >  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)
                                                                 cATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATg                                                                                                                                                                         >  1:58181/1‑251 (MQ=255)
                                                                      cataGGGATCTTCGGAAAGCTCAATCTCTGTCGGGGTAAATAACGAAGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCAt                                                                                                                                                                    <  2:125033/248‑1 (MQ=255)
                                                                         aGGGATCTTCTGCACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCCGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGGGCGCTGCTGGAGCAACTGAAGCGTCAACAAAgct                                                                                                                                                                                                                                               >  2:240974/1‑174 (MQ=255)
                                                                         aGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAgct                                                                                                                                                                                                                                               <  1:240974/174‑1 (MQ=255)
                                                                          gggATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGCTCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGGAGCGTCCACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACCCGAAGGCAATACCCACCAATGTGCATGGcc                                                                                                                                                                >  2:463057/1‑251 (MQ=255)
                                                                           ggATCTTCTGAAAGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGATATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCt                                                                                                                                                               >  1:112668/1‑251 (MQ=255)
                                                                                 tctGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATc                                                                                                                                                          <  1:253589/250‑1 (MQ=255)
                                                                                    gAACGCTCAATCGCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGTCTAAATCTg                                                                                                                                                        <  1:107304/249‑1 (MQ=255)
                                                                                     aaCGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGa                                                                                                                                                      >  1:666637/1‑250 (MQ=255)
                                                                                      aCGCTCAATCTATGTCGTGGGAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATACTGGTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGAGACGGTGGCGTTGACGGTACGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGATGAAGAATAAACATATCGCCTTAAGTGTCTGCGGTGTTGACAACTCGACGGACCTGCTCACCAATgggcagggcctacaagtggta                                                                                                                                                     >  2:420680/1‑230 (MQ=255)
                                                                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACt                                                                                                                                                 <  1:227325/249‑1 (MQ=255)
                                                                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTgg                                                                                                                                               <  1:256849/251‑1 (MQ=255)
                                                                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTg                                                                                                                                                <  1:51793/250‑1 (MQ=255)
                                                                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCGGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTTTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACt                                                                                                                                                 <  2:374563/249‑1 (MQ=255)
                                                                                                   gTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAAAATAAACATATCGACTTACGTGTCTGCGGGGTTTACAACTCGAAGGCACTGCTCACACATGTGCATGGCCTAAATCTGGAAAACTGCCAGga                                                                                                                                          >  2:576554/1‑248 (MQ=255)
                                                                                                     cGTGGTAAATAACGATGATGCGACCACTGGCTGGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCGGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTgaag                                                                                                                                                                                                                                       <  2:429854/154‑1 (MQ=255)
                                                                                                     cGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTgaag                                                                                                                                                                                                                                       >  1:429854/1‑154 (MQ=255)
                                                                                                        ggTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCTACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTgg                                                                                                                                  <  1:626377/251‑1 (MQ=255)
                                                                                                             aTAACGATGAAGCGACCACTGGCATGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGGGGCGTTGGCGGTGCGATGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTATCTGCGGTGTTGCCAACCCGAAGTACCTGCTCACCAAATGGACTGGACCAAATCTGGAAAACTGGCAGgaaaaaacggcccc                                                                                                                              >  2:444113/1‑239 (MQ=255)
                                                                                                               aaCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGGACTGGCGCAAGc                                                                                                                           <  1:501096/251‑1 (MQ=255)
                                                                                                                 cGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCa                                                                                                                         <  1:290569/251‑1 (MQ=255)
                                                                                                                      atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAAt                                                                                                                                                           >  2:3764/1‑212 (MQ=255)
                                                                                                                      atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAAt                                                                                                                                                           <  1:3764/212‑1 (MQ=255)
                                                                                                                      atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                                                                    <  1:268555/251‑1 (MQ=255)
                                                                                                                           aCCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATTGCCTAAATCTGGAAAACTGGCCGGAAGAACTGGCGCAAGCCAAAGAGCCGtt                                                                                                               >  2:20377/1‑251 (MQ=255)
                                                                                                                              aCTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTaa                                                                                                            <  1:20377/251‑1 (MQ=255)
                                                                                                                               cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACaaa                                                                                                                                                                                                                                                  >  2:600620/1‑117 (MQ=255)
                                                                                                                               cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACaaa                                                                                                                                                                                                                                                  <  1:600620/117‑1 (MQ=255)
                                                                                                                                 gtcgttctcgTTACTCTGCAGATGCTGTTCAAGTCCGATCAGGTTATCGACGTGTTTGTGATTGGCGTCGGTGTCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACTTGTCTTCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCGAAATCTTGAAAACTGGCAGGAAGAACTTGCGTAAGCCAAATAGCCGTTTAAt                                                                                                           <  2:34004/241‑1 (MQ=255)
                                                                                                                                    gTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCgg                                                                                                      <  1:249673/251‑1 (MQ=255)
                                                                                                                                          gTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTAGCGTCGGTGGCGTTGGCGGTGCGATGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACAAATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAAACTGGAAAACTGGCAGGAAGAACTGGACCAAGCCAAAGAACCGTTTAATATCGGGCGCtt                                                                                                >  2:477947/1‑251 (MQ=255)
                                                                                                                                              cTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATCTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGACAAAGAGCCGTTTAATCTAGGGCGCTTAAtt                                                                                            >  2:163745/1‑251 (MQ=255)
                                                                                                                                                catctGATGCTGGTCAATTCCGATCAGGTTATCGAAGTGTTTGGGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAATAATAAACATGTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCGTACCAATGTGCATGGCCTAAATCTGGAAAACGGGCAGGAAGAACTGGCGCAAGCTATAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                           <  2:655808/250‑1 (MQ=255)
                                                                                                                                                catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGGCGGTGGCGTTGGCGGTGCGTTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                           <  1:388932/250‑1 (MQ=255)
                                                                                                                                                catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCg                                                                                          <  1:247855/251‑1 (MQ=255)
                                                                                                                                                catcaGAGGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCg                                                                                          <  1:617089/251‑1 (MQ=255)
                                                                                                                                                    aGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTc                                                                                      <  1:609608/251‑1 (MQ=255)
                                                                                                                                                     ggtgtggTTCAATACCGAGCAGGTTTTCGAAGTGGTTGTGATTGGCTGCGGTGGCGTGGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACGGGCAGGAAGAACTGGCGCAAGCCAAAGAGCTGTTTAATCTCTGGCGCTTAATTCGCCTCg                                                                                     <  2:400836/245‑1 (MQ=255)
                                                                                                                                                                  aCCGATCAGGTTCTCGTAGTGTTTGTGATTCACGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAGAGCTGGCTGGGGAATAACCATATCGGCTTACGTGTCTGCGGTGTTGCCAACTCGCAGGCACTGCTCCCCACTGAGCATGGCCTAAATCTGGaa                                                                                                                                                     >  2:73224/1‑174 (MQ=255)
                                                                                                                                                                  aCCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGaa                                                                                                                                                     <  1:73224/174‑1 (MQ=255)
                                                                                                                                                                         aGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGAGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgctg                                                                 >  1:41477/1‑251 (MQ=255)
                                                                                                                                                                            ttATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGCATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACAAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTAGAGCAAGCAAAAGCGCAGTTTTATATCGGGCGCTTTATTCGCCTCGTGAAAGAATATCATCTGATGAAc                                                              >  2:29525/1‑251 (MQ=255)
                                                                                                                                                                            ttATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAc                                                              >  1:439176/1‑251 (MQ=255)
                                                                                                                                                                             tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                           <  1:125305/205‑1 (MQ=255)
                                                                                                                                                                             tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                           >  2:125305/1‑205 (MQ=255)
                                                                                                                                                                              aTCGAAGTGTTTGGGATTGGCGGCGGGGGAGTAGGCGGTGCGCTTCTGGTGCAACTGAATCGTCAACAAAGCGTGCTGATGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCATTGCTCACCAATGTGCATGGCCTAAATTTGGAAAACTGGCAGGAAGAACTGGCGCAAGACAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAcc                                                             <  2:375278/250‑1 (MQ=255)
                                                                                                                                                                                cGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCg                                                           <  1:156721/250‑1 (MQ=255)
                                                                                                                                                                                     tgtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCCAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGCAAGAATATCATCTGCTGAACCCGGTCAtt                                                     >  2:563405/1‑251 (MQ=255)
                                                                                                                                                                                      gtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAttg                                                    <  1:545447/251‑1 (MQ=255)
                                                                                                                                                                                            gATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACt                                              >  2:176529/1‑251 (MQ=255)
                                                                                                                                                                                            gATTGGAGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCCCCCCGTGAAAAAATATCATCTTCTGACCCCGGTCATTGTTGACt                                              >  2:280962/1‑251 (MQ=255)
                                                                                                                                                                                             aTTGGTGTCTGTGGCGTTGGGGGTTCGTTGATGGAGCAACTGAGGCGTCAATAACGCTGGCTGAAGAAGAAACATATCGTCTGACGTGTCTGCTGTGTTGCCAACTCGAAGGCACAGCACACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACGGGCGCAAGCCAAATAGCCGTTTAATCACGGGCGCGTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGGCATTGTTGACt                                              <  2:543366/250‑1 (MQ=255)
                                                                                                                                                                                                              tGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGACTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGCACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGccg                            >  2:261741/1‑249 (MQ=255)
                                                                                                                                                                                                                            tgGAGCAACTGAAGCGTCAACAAAGCTGGCTGCAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGt                                                                                                                                                                        >  2:352448/1‑97 (MQ=255)
                                                                                                                                                                                                                            tgGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGt                                                                                                                                                                        <  1:352448/97‑1 (MQ=255)
                                                                                                                                                                                                                               aGCAACTGAAGCGTCAACAAAGCTGTCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGcc           <  1:611862/251‑1 (MQ=255)
                                                                                                                                                                                                                               aGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGcc           <  1:344392/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATgacgg         >  1:89258/1‑247 (MQ=255)
                                                                                                                                                                                                                                     tGAAGCGTCAACAAAGCTGTCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGACCAATATGCCgccctc     >  2:673986/1‑246 (MQ=255)
                                                                                                                                                                                                                                        aGCGTCAACAAATATGGCTGAAGAATAATCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACGGGCGCAAGCCAACGAGACGTTTAATCTCGGGCGCTTAATTAGCCTCGTGAACGAATATCATCTGCTGAACCCGGTCCTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGGTCAATATGCCGACtttctt  >  2:378470/1‑247 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
ACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTG  >  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: