Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,646,970 G→T Q224H (CAG→CAT)  gor → glutathione reductase (NADPH)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,646,9700GT97.4% 252.9 / NA 76Q224H (CAG→CAT) gorglutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (42/32);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCAGCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGG  >  NC_000913/3646733‑3647215
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
cccGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACg                                                                                                                                                                                                                                          <  2:113905/251‑1 (MQ=255)
                gTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGGGATTAACGGCCTCGGCGCGGAAAAGCCTCTGTTTGTGCGTAAACCTGCGCCGCTGCGCAGCTTCGACCCGATGATTTTCGAAACGCTGGTCGCGGTGATGAACGCCG‑AAGGCCCGCCTCTGCAC‑ACCAACGCCATACCGCAAGCggg                                                                                                                                                                                                                          >  2:177083/1‑250 (MQ=255)
                   gAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGt                                                                                                                                                                                                                       <  1:124598/251‑1 (MQ=255)
                   gAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGGGATTTACGGCCTCGGCGCGAAAACGCCTCTTTTTGTGCGTAAAAATGGGCCCCTGCGCAGCTTCGACCCGATGGTTTCCGAAACGCTGGTCGAAGGGGTGAACGCCG‑AAGGCCCGCATCTGCCC‑ACCCACGCCCTCCCCAAAGCGGtggg                                                                                                                                                                                                                       >  2:345549/1‑248 (MQ=255)
                    aaTACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTg                                                                                                                                                                                                                      <  2:129800/251‑1 (MQ=255)
                        cGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGaaaa                                                                                                                                                                                                                  <  2:380830/251‑1 (MQ=255)
                            aTTGATTCTGATGGCTTCTTCGCCCTTCCTGCGTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATAc                                                                                                                                                                                                              <  2:195881/251‑1 (MQ=255)
                              tGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGGTGTTGGCCCGGGTTACATCGCCGTTGAGCTGGGGGGCGGGATTAAAGGGCCCGGCGCGAAAACGCATCTGTTTGTGCGGAAACATGCGCCGCCGCGCCGCTTCGGCCCGCTGATTTCCGAAACGCCGGTCGTAGGGAGGACCGCCG‑GAGGCCCCCATTTGCAC‑ACAAACCCCACCCCGCAAGCggggaggaaaactacgc                                                                                                                                                                                                            >  2:49757/1‑236 (MQ=255)
                               gATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGCGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGAATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCCGGTCCGAGGGATGAACGCCG‑GAGGCCCCCATCTGCAC‑ACCAACGCCATCCAGCAAGCGGTAGCGAAAAATACCg                                                                                                                                                                                                            >  1:51559/1‑250 (MQ=255)
                                  tCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGGTTTACCGCCTCGGCGCGAAAACGCATCTGTTTTGGCGGAAACATGCCCCCCTGCGCCGCCTCGGCCCGCTGATTTCCTAAACGCCGGTCGAAGTGATGAACGCCG‑AAGGCCCCCATCTGCCCCACCA‑CGCCCTCCCGAAAGCGGTAGTGAAAAATACCGATgg                                                                                                                                                                                                        >  2:869/1‑251 (MQ=255)
                                        ggCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCt                                                                                                                                                                                                  <  1:489423/251‑1 (MQ=255)
                                                      tCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTTAGCCGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTCAACATGCGCCCCTGCGCGGCCTCGACCCGATGGTTTTCGAAACCCTGGTCGAAATGAGGAACGCCG‑GAAGCCCGCATTTGCAC‑ACCAACGCCATCCCCAAAACGGGAGGGAAAAATACCGATgtggacctgacgcgggcgctgg                                                                                                                                                                                    >  2:248348/1‑230 (MQ=255)
                                                      tCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACgctggagctg                                                                                                                                                                                     <  1:134734/250‑1 (MQ=255)
                                                       ccTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAAAGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCCTCCCGGAAGCGGTAGGGGAAAATACCGAAGGTAGCCCGACgctgggggtggg                                                                                                                                                                                   >  1:239164/1‑244 (MQ=255)
                                                         tGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGGGCTGGCGGGGGTGATTAAACGCCCCGGCGGGAAAAAGCATCTGTTTGTGCGGAAACATGCGCCGCTGCGCCGCTTCGGCCCGATGATTTTCCAAACGCCGGGCGAAGTGGTGAACGCCGGAAGCCCCGC‑TCTGCCC‑ACCAAAGCCATCCCGAAAGCGGGAGGGAGAAATACCGATGGGAGCCCGACactggagctgga                                                                                                                                                                                   >  2:433807/1‑249 (MQ=255)
                                                                ccAGAGCGCGTGGCGGTTGTTGGCGCGGGTTTCATCGCCGTTGAGCTGGCGGGGGGGATTAACGGCCTCGGCGCGCAAACGCCTCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGGCCCGCTGATTTCCGGAACGCTGGTCGGAGTGATGAACGCCG‑AAGGCCCGCCTCTGCAC‑ACCCAAGCCATCCCGAAAGGGGTAGTGCAAAAAACCGCTGGGAGCCTGACGCTGGAGGTGGAAGATggg                                                                                                                                                                            >  2:84898/1‑248 (MQ=255)
                                                                         gTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCAGCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTGCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGGGAAAc                                                                                                                                                                 <  2:276141/251‑1 (MQ=255)
                                                                              ggTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGAATAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGCTGGTAGCCTGCCGCTGGAGCTGGAAGATGGGCGCCGTGAAACGGTgg                                                                                                                                                            >  1:40893/1‑251 (MQ=255)
                                                                                 tgttgGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGAGTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGAt                                                                                                                                                          <  2:113909/250‑1 (MQ=255)
                                                                                       cgcgGGTTACCTCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGACCGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTgg                                                                                                                                                   >  1:172799/1‑250 (MQ=255)
                                                                                         cgGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCATCGGCCAGAAAACGCATCTGTTTGTGCGTAAAAATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑GAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGGGAAAAATACCGATGGTAGGCTGACGCTGGAGCTGGAAGATGGTCGCAGTTACACGGTGGATTGCCTGAtt                                                                                                                                                 >  1:320835/1‑251 (MQ=255)
                                                                                                cATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAAGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATAACGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTAGCAGTGAAACGAGGGATTGCCTGATTTGGGCgg                                                                                                                                          >  2:5273/1‑250 (MQ=255)
                                                                                                    gCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTgg                                                                                                                                      <  1:869/251‑1 (MQ=255)
                                                                                                      cGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTc                                                                                                                                    <  1:177083/251‑1 (MQ=255)
                                                                                                       gTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTcg                                                                                                                                   >  1:115449/1‑251 (MQ=255)
                                                                                                         tGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGCTGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTcgcg                                                                                                                                 >  1:22636/1‑251 (MQ=255)
                                                                                                         tGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTcgc                                                                                                                                  <  2:308644/250‑1 (MQ=255)
                                                                                                                   ggcgTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCaa                                                                                                                       <  2:115449/251‑1 (MQ=255)
                                                                                                                        gATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGac                                                                                                                   >  1:443307/1‑250 (MQ=255)
                                                                                                                         aTTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGAGTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCAGCCCGAAAGCGGTAGTGAAAAATACCGAGGGTAGCCGGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGacaa                                                                                                                 <  2:267937/251‑1 (MQ=255)
                                                                                                                            aaCGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGATGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAAAGCCG‑AAGACCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTCGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCATGCCAATGACAACAt                                                                                                              >  1:481671/1‑251 (MQ=255)
                                                                                                                             aCGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAAAGGTGGATTGCCTGATTTTGGCGCTTGGTCGCGAGCCTGCCAATGACAACATc                                                                                                             >  1:17154/1‑251 (MQ=255)
                                                                                                                              cGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGACATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATc                                                                                                             >  2:387799/1‑250 (MQ=255)
                                                                                                                               ggCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCCCATCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTTCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGATCTGGAAGATGTTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCaa                                                                                                           <  1:308943/251‑1 (MQ=255)
                                                                                                                                         cgAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTGGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGcc                                                                                                 <  2:486006/251‑1 (MQ=255)
                                                                                                                                                      gTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGACTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTaaaa                                                                                    <  2:123286/251‑1 (MQ=255)
                                                                                                                                                          gtgCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTa                                                                                 <  1:345549/250‑1 (MQ=255)
                                                                                                                                                             cGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGa                                                                             <  1:130992/251‑1 (MQ=255)
                                                                                                                                                               tAAACATGCGCCGCTGCGCCGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGGCGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTACCGaa                                                                            >  2:137392/1‑250 (MQ=255)
                                                                                                                                                                aaaCATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGaaa                                                                           >  2:163632/1‑250 (MQ=255)
                                                                                                                                                                aaaCATGCGCCGCTGCGCAGCTTCGACCCGATGATTGCCGAAACGCTGGTCGAAGTGAGGAACGCCT‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGAGTTGGGCGATTGGTCGCGAGCCTGTCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGaaaa                                                                          <  1:518690/251‑1 (MQ=255)
                                                                                                                                                                 aaCATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCGGGAAGATGGTCGCAGTGAAACGGTGGATTGCCGGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGaaaaa                                                                         <  2:186082/251‑1 (MQ=255)
                                                                                                                                                                 aaCATGCGCCGCTGCGCAGCTTCGACCCGATGATCTCCGAAACGCTGGTCGAAGTTATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCATCCTGGAAGCCGCTGGCGTTAAAACTAACGaaaaa                                                                         <  1:22719/251‑1 (MQ=255)
                                                                                                                                                                   cATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAgg                                                                       >  1:84702/1‑251 (MQ=255)
                                                                                                                                                                       cgccgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGCTTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCtat                                                                   >  2:498658/1‑251 (MQ=255)
                                                                                                                                                                       cgccgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCtat                                                                   >  2:42312/1‑251 (MQ=255)
                                                                                                                                                                        gccgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCtat                                                                   >  1:80494/1‑250 (MQ=255)
                                                                                                                                                                          cgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCTAGCCTGCCAATGACAACATCAACCTGGAAGCCGCGGGCGTTAAAACTAACGAAAAAGGCTATAtc                                                                <  1:498658/251‑1 (MQ=255)
                                                                                                                                                                          cgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtc                                                                >  1:337301/1‑251 (MQ=255)
                                                                                                                                                                          cgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtc                                                                <  1:301720/251‑1 (MQ=255)
                                                                                                                                                                            cTGCGCCGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtcgt                                                              >  2:293455/1‑251 (MQ=255)
                                                                                                                                                                              gcgcAGCGTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGGGAAAAATACCGATGGTAGCCTGACGCGGGAGCTGGAAGAGGGGCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCGGGCGTTAAAACTAACGAAAAAGGCTATAtcgtcg                                                            <  1:207215/251‑1 (MQ=255)
                                                                                                                                                                                   gCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGACG‑AACGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGAAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAAACTAGAAGCCGCTGGCGTTAAAACACACGAAACAGGCTATAtcgtcgccgct                                                       >  1:233964/1‑246 (MQ=255)
                                                                                                                                                                                     ttCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGATCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAACAAGGCTATATCGTCGTCGATa                                                      >  1:2148/1‑250 (MQ=255)
                                                                                                                                                                                      tCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATa                                                      <  1:163632/249‑1 (MQ=255)
                                                                                                                                                                                          cccGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATc                                                <  1:385591/251‑1 (MQ=255)
                                                                                                                                                                                           ccGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGCTTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATc                                                >  2:421302/1‑250 (MQ=255)
                                                                                                                                                                                            cGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCa                                               >  2:115393/1‑250 (MQ=255)
                                                                                                                                                                                              atgatTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCACTGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATCTCGTCGTCGATAAATATCaaa                                             >  2:476389/1‑250 (MQ=255)
                                                                                                                                                                                                gatTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAac                                           >  1:395452/1‑250 (MQ=255)
                                                                                                                                                                                                 atTTCCGCAACGCTGGTCGAAGTGATGAACGCCG‑ACGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAacac                                         >  2:548497/1‑251 (MQ=255)
                                                                                                                                                                                                   ttCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACAcc                                        >  1:123228/1‑250 (MQ=255)
                                                                                                                                                                                                     ccGAAACGATGGCCGAAGTGATGAACGCCG‑AAGGCACGCATCTGCAC‑ACCAACGCCATCACGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGTCGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTTGGCGATTGGGCGCCAGCCTGCCAAAGACAACAGCATCCTGGAAGCCGCTGGCATTAAAACTAACGAAAAAGGATATATCGTCGGCGATAAATATCAAAACCCCaa                                      <  2:552088/250‑1 (MQ=255)
                                                                                                                                                                                                        aaaCGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAACATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGCTTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATCTCAAAACACCAATAtt                                  >  2:525629/1‑251 (MQ=255)
                                                                                                                                                                                                               ggTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGt                            <  2:392444/250‑1 (MQ=255)
                                                                                                                                                                                                               ggTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGt                            >  1:518331/1‑250 (MQ=255)
                                                                                                                                                                                                               ggTCGAAGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGAAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGACAAAGGCTATATCGTCGTCGATAAATATCACAACACCAATATTGAAGGt                            >  1:85562/1‑250 (MQ=255)
                                                                                                                                                                                                                    aaGTGATGAACGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGGTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACGCCAATATTGAAGGTATTTAc                      <  1:548497/251‑1 (MQ=255)
                                                                                                                                                                                                                           gaACGGCG‑AAGGCCCGCAGCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCGGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGATGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTgg               <  1:389464/251‑1 (MQ=255)
                                                                                                                                                                                                                           gaACGCCG‑AAGGCCCTCAGCTGCAC‑ACCAACGTCATCCAGAAGGCGGGAGTTAAAAATACCTGTGTTAGCCTGACGCTGGATCTGTAAGATGGTCTCAGTGTAACGGTGGAGTGCCTGATTTGGGTGATTGGTCGCGAGCCTGCCAATTACAACATCAACCTGGAAGCCGCTGGTGTTAAAACTAACGAAAAAGGCTATGGCGTCTTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTgg               <  1:273545/251‑1 (MQ=255)
                                                                                                                                                                                                                             aCGCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGc             <  1:514855/251‑1 (MQ=255)
                                                                                                                                                                                                                               gCCG‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGGTGGAAGATGGTCGAAGTGAAACGGTGGATTGCATGATTTGGGCGATTGGTCGCGAGCCTGCCAAGGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGa           <  1:298199/251‑1 (MQ=255)
                                                                                                                                                                                                                                  g‑AAGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGCTTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATa         >  1:271770/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aaGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAac       >  2:14640/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aGGCCCGCATCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAac       >  2:436206/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ccGCATCTGAAC‑ACCCACGCCATCCAGAAGGCGGTAGCGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACggg  >  2:322418/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
CCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCG‑AAGGCCCGCAGCTGCAC‑ACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGG  >  NC_000913/3646733‑3647215

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