Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,224,762 G→A G312D (GGC→GAC)  metH → cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,224,7620GA97.7% 270.7 / ‑6.6 87G312D (GGC→GAC) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (42/43);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGG‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAAC  >  NC_000913/4224517‑4225006
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
aTTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑cacca                                                                                                                                                                                                                                                  <  1:190367/251‑1 (MQ=255)
 ttttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGCACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑caccac                                                                                                                                                                                                                                                 >  1:230177/1‑251 (MQ=255)
  ttttcAACTCATTGCTCCCCGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGTGCAAGCGGGTTTTCTCAATATCGTCGGCGGCGGCTGTGA‑caccac                                                                                                                                                                                                                                                 <  1:240984/246‑1 (MQ=255)
      cAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGCCGGCGGCCGCTGCGA‑CACCACGCcac                                                                                                                                                                                                                                            >  1:398277/1‑251 (MQ=255)
       aaCTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCcac                                                                                                                                                                                                                                            >  2:101538/1‑250 (MQ=255)
           cATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACat                                                                                                                                                                                                                                       >  1:64666/1‑251 (MQ=255)
              tGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATAtt                                                                                                                                                                                                                                    >  2:399191/1‑251 (MQ=255)
                cgcCACGCCGGAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGCACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGAACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCACGACGCCGACACGATGGCAAAACAGATACGTGCATGGGCGCAAGCGGGTTTTCTCAATCTCGTCGGCGGCTGCTGTGA‑CACCCCGCCACAACATATTGc                                                                                                                                                                                                                                  >  2:353943/1‑251 (MQ=255)
                 gcCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTga                                                                                                                                                                                                                                  >  1:301886/1‑249 (MQ=255)
                     cGCCGAAGCTCTGACCTTTGGCCTGACCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGGCACCGCGCACCCGAACGCCGGGCTACCCATCGCCTTTTGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATTGCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGa                                                                                                                                                                                                                             <  1:273938/251‑1 (MQ=255)
                      gCCGAAGCTCTGACCTTTTGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGAt                                                                                                                                                                                                                            <  2:468914/251‑1 (MQ=255)
                       ccGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGCTGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGCACGCCGGGCTACCCAACGCCTTTGGTG‑AGTCCGATCTCGACGCCGACACGATGGCAAAGCAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATCTCGTCGGCGGATGCCGTGA‑CACCCCGCCACTACGTATAGCAGCgtgg                                                                                                                                                                                                                           >  2:222399/1‑248 (MQ=255)
                            gcgctGACCGTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCg                                                                                                                                                                                                                      <  1:537908/248‑1 (MQ=255)
                               ctgtccTTTTGCCTGAACTGTGCGCTGGGGGCCGATGAACTGCGCCAGTCCGTGCAGTAGCTGTCGCGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGGCGt                                                                                                                                                                                                                     <  2:423477/245‑1 (MQ=255)
                                      ttGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAg                                                                                                                                                                                                              <  1:181671/249‑1 (MQ=255)
                                               aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGCGTCGGGCAGTAGAAGGCTTAgcc                                                                                                                                                                                                   >  2:392469/1‑250 (MQ=255)
                                                  gtgCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGcgccg                                                                                                                                                                                                >  1:101828/1‑251 (MQ=255)
                                                    gcgcTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCcgcg                                                                                                                                                                                              <  1:13814/251‑1 (MQ=255)
                                                             cccGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGcc                                                                                                                                                                                     >  2:465238/1‑251 (MQ=255)
                                                             cccGATGAACTGCGCCAGTACGTGCAGGAGCTGCCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTGGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATGCGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACAAATTGCAGCGAGGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCGAACTGcc                                                                                                                                                                                     <  2:505816/251‑1 (MQ=255)
                                                                     aCTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTc                                                                                                                                                                             <  2:151792/251‑1 (MQ=255)
                                                                        gcgcCAGTACGGGCAGGAGCTGTCACGGATTGCGGAATGCTTCGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTccc                                                                                                                                                                           <  2:522793/250‑1 (MQ=255)
                                                                         cgcCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGt                                                                                                                                                                         >  2:130125/1‑251 (MQ=255)
                                                                         cgcCAGTACGTGCAGGAGCTGTCACGGATTGCGGAAGGCTCCGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGt                                                                                                                                                                         <  2:178464/251‑1 (MQ=255)
                                                                             aGTACGTGCAGGAGATGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCCGTAGAAGGATTAGCGCCGCGACAACTGCCGGCAATTCCCGTaccc                                                                                                                                                                     >  2:144336/1‑248 (MQ=255)
                                                                                 cgggCAGGGGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATCCGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTTGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑cc                                                                                                                                                                <  2:198066/248‑1 (MQ=255)
                                                                                 cGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑cc                                                                                                                                                                <  1:38978/251‑1 (MQ=255)
                                                                                    gcgggAGCTGTCACGGATTGCGGAATGCTTCGTCACGGCGCCCCCGAACCCCGGGCTCCCAAACGCCTTGGGGG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATAGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTTTGA‑CCCCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATGAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGtt                                                                                                                                                             <  2:401915/248‑1 (MQ=255)
                                                                                    gCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGtt                                                                                                                                                             <  2:481051/251‑1 (MQ=255)
                                                                                      aGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTg                                                                                                                                                           >  1:165275/1‑251 (MQ=255)
                                                                                      aGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTg                                                                                                                                                           >  1:460800/1‑251 (MQ=255)
                                                                                      aGGAGCTGTCACGGATGGCGGAATGCGACGCCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGGG‑AGTACGATCTCGACGCCTACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATAGTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTg                                                                                                                                                           <  2:42898/251‑1 (MQ=255)
                                                                                        gAGCTGTCAAGGCTTGCGGAATGCCACGTCAACGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGCTGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTc                                                                                                                                                         >  2:313808/1‑251 (MQ=255)
                                                                                         aGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTcc                                                                                                                                                        <  1:82019/251‑1 (MQ=255)
                                                                                          gCTGGCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACTCGATGGCAAAACAGATACGGGAATGGGCGCAAGCGGGTTTTCTCAATAGCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATAGTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCg                                                                                                                                                       <  1:488927/251‑1 (MQ=255)
                                                                                               cACGGATTGCGGAAGGCTACGTCACCGCGCCCCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CCCCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGGAGCCTG‑CCGTTTGTCCGGCCt                                                                                                                                                   <  2:545068/250‑1 (MQ=255)
                                                                                                aCGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCTCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTgg                                                                                                                                                 <  1:85201/251‑1 (MQ=255)
                                                                                                aCGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTgg                                                                                                                                                 >  1:445482/1‑251 (MQ=255)
                                                                                                 cGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTgg                                                                                                                                                 >  1:201882/1‑250 (MQ=255)
                                                                                                      tGCGGAATGCTACGTCACCCCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGGGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCg                                                                                                                                            <  2:466365/250‑1 (MQ=255)
                                                                                                        cGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGGG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTGTCTCAAGATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTg                                                                                                                                         <  2:177526/251‑1 (MQ=255)
                                                                                                         ggAATGCTACGTAACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCTGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGCAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTg                                                                                                                                         >  1:266437/1‑250 (MQ=255)
                                                                                                          gAATGCTACGTCACCGCGCCCCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGaa                                                                                                                                       <  1:69748/251‑1 (MQ=255)
                                                                                                           aaggCTACTTCACCGCGCGCCCGAACGCCGGGCGACCCAACGCCTT‑GGTGGAGTACGATCTCGCCGCCGACACGATGGCAAAACAGATACGGGAAGGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCT‑TGGACACCGCGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACGGCCGGAAATGCCCGTAGCCTG‑CCGTTTGGCCGGCCTGGAGCCGCTGAAc                                                                                                                                      <  2:28803/248‑1 (MQ=255)
                                                                                                                   gTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCGGTA‑CCTGGCCGTTTG‑CCGGCCTgggagcgccggacccttggcg                                                                                                                              >  2:277952/1‑231 (MQ=255)
                                                                                                                     cACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGa                                                                                                                              >  2:416581/1‑249 (MQ=255)
                                                                                                                             ccccGACCGCCGGGCTACCCAACGCCTTGGGTG‑AGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAGCGTCGGCGTCTGCTGTGA‑CGCCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTg                                                                                                                    <  1:68116/250‑1 (MQ=255)
                                                                                                                                cGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGttt                                                                                                                 <  1:323280/251‑1 (MQ=255)
                                                                                                                                cGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGttt                                                                                                                 >  2:107039/1‑251 (MQ=255)
                                                                                                                                  aaCGCAGGGCTACCCAACGCCTTTGGTG‑CGTACGATCTCGCCGCCGACACGATGGCAAAACAGATCCGTGAATGGGCGCAACCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGCACATTCGCGAAGATAGCCTTTTtgt                                                                                                               >  2:268911/1‑251 (MQ=255)
                                                                                                                                     gCCGGGCTCACCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGCCACGATGGCCACACAGCTACGTGCATGGGCGCACGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGACGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGa                                                                                                             >  2:285845/1‑250 (MQ=255)
                                                                                                                                        gggCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGt                                                                                                         <  1:130125/251‑1 (MQ=255)
                                                                                                                                         ggCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGACTGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGCAGATAGCCTGTTTGTGAACggg                                                                                                        >  2:50376/1‑249 (MQ=255)
                                                                                                                                             aCCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGt                                                                                                     <  1:388101/250‑1 (MQ=255)
                                                                                                                                             aCCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGt                                                                                                     <  2:101828/250‑1 (MQ=255)
                                                                                                                                             aCCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTg                                                                                                    >  2:164521/1‑251 (MQ=255)
                                                                                                                                               ccAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAAGGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGaa                                                                                                  <  2:150447/251‑1 (MQ=255)
                                                                                                                                                             g‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGGGAACGCACCAACGTCAc                                                                                    >  1:147155/1‑251 (MQ=255)
                                                                                                                                                             g‑AGTACGATCCCGACGCCGACACGATGGCAAAACAGATACGTGAAGGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCAc                                                                                    <  2:165654/251‑1 (MQ=255)
                                                                                                                                                               agtgcGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAGCGTCGGCGGCTGCGGTTA‑CGCCACGCCACAACATATTGCAGCGATGAGACGTGCAGTAGAAGGATGAGCGCCGCGCGAACGGCTGGAAATTCCCGTAGCCGG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCAcc                                                                                   <  1:548575/247‑1 (MQ=255)
                                                                                                                                                               aGTACGATCGCGACGCTGACAGGATGGCAAAACGGATACGTGAAGGTGCGCAGGCGGGTTTTCTCAACAGCGGCGGCGGCTGGTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGGAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCATCAACGTCAc                                                                                    <  2:539571/250‑1 (MQ=255)
                                                                                                                                                               aGTACGAGTTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGTGCCGCGCAAACTGCCGGAAATGCCCGTAGCCTG‑CCGTTTGGGCGGTCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCAcc                                                                                   <  1:347956/251‑1 (MQ=255)
                                                                                                                                                               aGTACGAGATCGACTCCGACACGATGGCAAAACAGAAACGTGAATGGGCGCAATCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGTTGAACGCTCCAACGTCAcc                                                                                   <  2:372376/251‑1 (MQ=255)
                                                                                                                                                                          ggcGCCGACACGATGGCAAAACAGTAACGTGAAGGGGCGCAAGCGGGTTTCCTCAATATCGTCGGCGGCGGCTGTGA‑CACCACGCCACAACATATGGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGGTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGGGAACGCACCAACGTCGCCGGTTCCGCTaa                                                                        <  2:223151/249‑1 (MQ=255)
                                                                                                                                                                            cGCCGACACGATGGCAACACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACCTATTGCAGCGGTGAGTCGTGCAGTAGAGGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGCGCAGCTGAACATTGGCGACGATAGCCTGTTTGTGCACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGt                                                                      >  1:532309/1‑251 (MQ=255)
                                                                                                                                                                             gCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCAACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGCTTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGACGAAGATAGCCTGTTTGTGACCGTGGGTGAACGCACCAACGCAACCGGTTCCGCTAAGtt                                                                     >  1:87605/1‑251 (MQ=255)
                                                                                                                                                                              ccGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CCCCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTc                                                                    >  1:524477/1‑251 (MQ=255)
                                                                                                                                                                                     gATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGCACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGcc                                                             >  1:13359/1‑251 (MQ=255)
                                                                                                                                                                                      aTGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAGCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                            <  1:244571/251‑1 (MQ=255)
                                                                                                                                                                                        ggCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTACGTTCAAGCGCCTg                                                           >  2:132858/1‑250 (MQ=255)
                                                                                                                                                                                        ggCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTg                                                           >  1:76272/1‑250 (MQ=255)
                                                                                                                                                                                        ggCAAAACAGATACGTGAAGGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTg                                                           <  2:445482/250‑1 (MQ=255)
                                                                                                                                                                                           aaaaCAGATACGTGAATGGGCGCAAGCGGGTTTTCTCCATATCGTCGGCGGCTGATGTGA‑CACCCCGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTTAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATc                                                        >  2:247898/1‑250 (MQ=255)
                                                                                                                                                                                              aCAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaag                                                    <  1:250293/251‑1 (MQ=255)
                                                                                                                                                                                                   tACGTGAGTGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGGAGACGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGCACGCACCAACGTCACCGGTTCCGCTACGTTCAAGCGCCTGATCAaagaag                                                 >  2:220128/1‑249 (MQ=255)
                                                                                                                                                                                                    aCGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATAGTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAgaga                                              <  1:33257/251‑1 (MQ=255)
                                                                                                                                                                                                        gAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgcg                                                               >  1:194082/1‑230 (MQ=255)
                                                                                                                                                                                                        gAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgcg                                                               <  2:194082/230‑1 (MQ=255)
                                                                                                                                                                                                             ggCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATCGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCg                                     >  2:174185/1‑251 (MQ=255)
                                                                                                                                                                                                                      gggTTTTCTCAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGGGAACGCACCAACGTCACGGGTTCCGCTAAGGTCACGCGCCTGATCAAAGAAGAGAAATACAGCGAGGAGCTgg                            >  1:280354/1‑251 (MQ=255)
                                                                                                                                                                                                                            tctcAATATCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcg                      >  1:528968/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAATATCGTCGGCGGCTGCGGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCACGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgcgc                   <  1:430161/251‑1 (MQ=255)
                                                                                                                                                                                                                                   atCGTCGGCGGCTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTc                >  1:15606/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ggcggcTGGTGTTG‑CAACACAACACCACATATTGCAGCGATGAGCCGTGCAGTAGAAGGATTAGCGCCGCGCAACCTGCCGGAAATTCCCGTAGCCTG‑CCGATTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGCAAGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCACAGCAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAggg         >  2:262893/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gcggcgGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCACGTAGCCTG‑CCGATTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGAGCTGGATGTCGCGCGTCAACAGGTg        >  2:179442/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cggcTGCTGTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTgg       <  1:366292/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gctgTGA‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACAGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAAc  >  1:230873/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
ATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTG‑AGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGG‑CACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTG‑CCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAAC  >  NC_000913/4224517‑4225006

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