Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,164,123 C→G A162G (GCA→GGA)  btuB → cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,164,1230CG97.7% 282.3 / ‑3.1 87A162G (GCA→GGA) btuBcobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base G (38/47);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGC‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATG  >  NC_000913/4163879‑4164368
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
ggTCAACTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGATCCCGGAACGGAAATTTCAGG‑AGGGTg                                                                                                                                                                                                                                                  <  1:362407/251‑1 (MQ=255)
  tCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTggg                                                                                                                                                                                                                                                >  1:228126/1‑251 (MQ=255)
  gcAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTggg                                                                                                                                                                                                                                                <  1:911961/250‑1 (MQ=255)
    aGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAAACCGGAACGGAAATTTCAGG‑Agggggggg                                                                                                                                                                                                                                               >  1:773431/1‑250 (MQ=255)
        ctcATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGTCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGGGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAg                                                                                                                                                                                                                                            <  2:267093/249‑1 (MQ=255)
           aTCTATTTTTATTCGCGGTACAAAAGCCAGTCATGTGATGGTGTTAATTGATGGAGTACGCCTGAATCTGGCGGGGGTGAGTTGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCCGCGTGTTGAATATATCCGTGGGCCGCGCTCAGCTGTTTATGGGTCCGAGGCAATAGGCGGGGAGG‑TGAGTATCATCACAAAGCG‑CGCTGAACTCGGACCAGAAATTTAAGG‑AGGGTGGGGaaacact                                                                                                                                                                                                                                        >  1:679635/1‑245 (MQ=255)
            tCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGGTCCGGTGCAATAGGCGGGGGTG‑TGAATATCAGTGCGAAGAG‑CGATGACAACGGAACGGAAATTTTAGGGAGGGTGGGGagcagt                                                                                                                                                                                                                                         >  2:534381/1‑244 (MQ=255)
             cTATTTTTATTCTCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGGGGGGAAGCAAT‑AGt                                                                                                                                                                                                                                    >  1:169534/1‑251 (MQ=255)
               aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCTTGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTTAATATATCCGTGGGCCGCGCCCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATAGCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGtt                                                                                                                                                                                                                                   <  1:593675/250‑1 (MQ=255)
               aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑Ag                                                                                                                                                                                                                                     <  1:822565/248‑1 (MQ=255)
               aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGAGGGGAAGCAAT‑AGttt                                                                                                                                                                                                                                  >  1:267490/1‑250 (MQ=255)
                tttttATTCGCGGTACAACTGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAAAATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCCGG‑AGGGTGGGGAAGCAAT‑AGtttt                                                                                                                                                                                                                                 >  1:59985/1‑249 (MQ=255)
                 ttttATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATc                                                                                                                                                                                                                                >  1:468486/1‑251 (MQ=255)
                   ttATGCGCGGTACAAATGCCAGTTATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTGCCCGATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGTTCCGCTGTGTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGATGCAAT‑AGTTATCa                                                                                                                                                                                                                               <  2:749685/250‑1 (MQ=255)
                      ttCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGGAATTTCAGG‑AGGGGGGGGAAGCAAT‑AGTTATCAGAAc                                                                                                                                                                                                                           >  1:325827/1‑251 (MQ=255)
                        cgcgGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGGAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACt                                                                                                                                                                                                                          >  1:49908/1‑250 (MQ=255)
                             tACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGGTGAACCCGGAACGGGAATTTCAGG‑GGGGGGGGGAAGGAAT‑AATTATCACAAATatgggt                                                                                                                                                                                                                    >  2:945693/1‑248 (MQ=255)
                                   tGCCAGACATGTGTTGGAGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCAAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTTTGGTTCCGATGCAATCGGAGGGGGGG‑GGGAACTTATCACTA‑GCGGCGATGAACCCGGAACGCGAATTTCAGG‑TGGGGGGGGAAGCAATCAGTTCT‑AAAACTATGATGTatctt                                                                                                                                                                                                              >  2:149009/1‑250 (MQ=255)
                                     ccAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACg                                                                                                                                                                                                            <  1:51784/251‑1 (MQ=255)
                                       aGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCCGAACTATGATGTCTCTACgc                                                                                                                                                                                                           >  1:633550/1‑250 (MQ=255)
                                          cATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAAACCGGAACGGAAATTTCAGG‑GGGGTGGGGAAGCAAA‑AGGTTTCAGAACTAAGATGtttttacacagc                                                                                                                                                                                                        >  2:165061/1‑239 (MQ=255)
                                            tgtgTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGcaac                                                                                                                                                                                                     >  1:229370/1‑251 (MQ=255)
                                             gtgtTGGTGGTAATTGATGGCGGCCGCCTGAATCTGGCGGGGGTGAGTGGTTCGGCCGACCTGAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CTATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                                                                                                                                                                                    <  2:768808/251‑1 (MQ=255)
                                                  ggTGTTAATTGATTGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                                                                                                                                                                                <  2:733536/250‑1 (MQ=255)
                                                       tAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGCTGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCCGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGACCTATGATGTCTCTACGCAACAACAACTGGGGGa                                                                                                                                                                                           >  1:213856/1‑250 (MQ=255)
                                                        aaTTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt                                                                                                                                                                                          <  1:634150/250‑1 (MQ=255)
                                                            gATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                                                                                                                                                                                     >  1:623385/1‑251 (MQ=255)
                                                               ggCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                                                                                                                                                                                  <  2:169534/251‑1 (MQ=255)
                                                                      gCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGGATATAATCACGACGCG‑CGATGGACCCGGGACGGAAATTTCAGG‑AGGGTGGGGAAGCAAA‑AGGTATCAGAACTATGATGTCTCTACGCAGAAAAAACTGGGGGGAAAGACACGGGTAACg                                                                                                                                                                           >  2:349467/1‑251 (MQ=255)
                                                                          gAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCCCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGAG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt                                                                                                                                                                         <  1:534381/249‑1 (MQ=255)
                                                                           aaTCTGGCTGGGGTGTGGGGTTCTGCCGCCCTTATCCTGCTCCCTTGTGCGCTTGTCCAGCGTGTTGAAGATATCCGGGGGCCGCTCTCAGCTGTTTATGGTGCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACTACGTG‑CGATGAACTCGGAACGGAAATTGCAGG‑AGGGTGGGGAAGCAAT‑AGTAATCATAACTATGATGTCTCTACGCAGCAACAACGTGGGGAGAAGACACGGGTAACGCTGtt                                                                                                                                                                      <  2:558778/251‑1 (MQ=255)
                                                                                ggCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGCTGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCCAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGGAAAGACACGGGTAACGCTGTTGGGCg                                                                                                                                                                 >  2:330393/1‑251 (MQ=255)
                                                                                ggCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTGGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGAAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg                                                                                                                                                                   <  1:72567/249‑1 (MQ=255)
                                                                                  ggggggTGAGTGGTTCTGCCGACCTGAGCCAGTTCCCTATTGCGCTTGTCCAGCTTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAt                                                                                                                                                               <  2:315918/250‑1 (MQ=255)
                                                                                         gAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGcc                                                                                                                                                         <  1:259083/250‑1 (MQ=255)
                                                                                             ggttaggCCGACATGAGCCAGTTCGCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATAc                                                                                                                                                    <  2:151510/245‑1 (MQ=255)
                                                                                                   gCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATgg                                                                                                                                              <  2:371722/251‑1 (MQ=255)
                                                                                                    ccGACCTTAGGCAGTTGCCTATTGGCCTTTTCCAGGGTGTTGAAAATAATCGTGGGACGCGGCACGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGTGA‑TATCATCACGACGCG‑GGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGGCAT‑AGTTAACAGAACTATTGTGTCTCTACGCAGCAACCACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGt                                                                                                                                             <  2:193379/251‑1 (MQ=255)
                                                                                                      gACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGtt                                                                                                                                            >  1:409505/1‑250 (MQ=255)
                                                                                                           tAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTatgatg                                                                                                                                      >  1:403461/1‑251 (MQ=255)
                                                                                                           tAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTatgat                                                                                                                                       <  1:346373/250‑1 (MQ=255)
                                                                                                            aGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTTAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGAtgt                                                                                                                                     <  1:733173/251‑1 (MQ=255)
                                                                                                              ccAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGC‑AGGGTGTGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGAtgttt                                                                                                                                   >  1:247312/1‑250 (MQ=255)
                                                                                                               cAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGAtgttgt                                                                                                                                  <  1:47187/251‑1 (MQ=255)
                                                                                                                    cccTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGcc                                                                                                                              <  1:986381/250‑1 (MQ=255)
                                                                                                                    cccTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCt                                                                                                                             <  1:952834/251‑1 (MQ=255)
                                                                                                                    ccatATTGCGCGTGTCCAGCGGGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCt                                                                                                                             <  1:18656/248‑1 (MQ=255)
                                                                                                                     ccTATTGCGCTTGTCCAGCGTGTTGCATATACCCGTGGGCCCCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGACTATCATCACCACGTG‑CCATGACACCGGAACGGCAATTTCAGG‑AGGGGGGGGCAGGCAT‑CGTTATTCGAACTATGAACTCTCTACACCACAAACACTGGGGGCTAAGTCCCGGGCATCGCAGTTGGGCGATTATGCCCATACTCATGGTTATTAtgtttgttccc                                                                                                                             >  2:120594/1‑243 (MQ=255)
                                                                                                                     ccTATTGCGCTTGGCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑GGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTa                                                                                                                            <  1:782173/251‑1 (MQ=255)
                                                                                                                        aTTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGAAACGGGTAACGCTGTTGGGCGATTATGCCCATACTGATGGTTATGATGTTGTTGCCTATgg                                                                                                                         >  1:6477/1‑251 (MQ=255)
                                                                                                                            cgcTTGTCCAGCGTGTTGAATATATCCGTTGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑TGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAAt                                                                                                                     <  1:904516/251‑1 (MQ=255)
                                                                                                                                   ccAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGa                                                                                                               <  1:753255/250‑1 (MQ=255)
                                                                                                                                   ccAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGG‑TGAATATCATCACGACGCG‑CGCTGAACCAGGAACGGAAATTTCAGG‑AGGGGGGGGAAGAAAT‑AGTTATCAGAACTATGATGTCTCTACGGCGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGGTGCCTATGGTAATACCGGa                                                                                                               >  1:464519/1‑250 (MQ=255)
                                                                                                                                     aGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACg                                                                                                            >  1:399974/1‑251 (MQ=255)
                                                                                                                                     aGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGG‑GGAATATCATCACGACGCG‑CGCTGCACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACAAGGGTAACGCTGTTGGGCGATTATGCCCCTACTCCTTGTTATGATGTTGTTGCCTATGGTAATACCGGaaac                                                                                                            >  2:398263/1‑249 (MQ=255)
                                                                                                                                        gtgtTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAAACCGGAACGGAAATTTCAGG‑AGGGGGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTTTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTTTGATGTTGTTGCCTATGGTAATACCGGAACGCa                                                                                                          >  1:509871/1‑250 (MQ=255)
                                                                                                                                             gAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa                                                                                                                                                             <  2:51097/194‑1 (MQ=255)
                                                                                                                                             gAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa                                                                                                                                                             >  1:51097/1‑194 (MQ=255)
                                                                                                                                               atatatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGcaga                                                                                                  <  1:610619/251‑1 (MQ=255)
                                                                                                                                                     ccGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGAt                                                                                             <  1:667358/250‑1 (MQ=255)
                                                                                                                                                        tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACg                                                                                         >  1:304415/1‑251 (MQ=255)
                                                                                                                                                         gggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGa                                                                                        <  1:540839/251‑1 (MQ=255)
                                                                                                                                                           gtcgtgcTCCGCTATTTAGGGGTTCGATGCAATAGGAGGGGTGG‑TGAATATGATCACGACGCG‑CGATGAACCCGGAACGTAAATTTCAGG‑AGGGGGGGGAAGCTAT‑AGTTATCAGAACTATGATGTCGCGACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGTTAATACCGGAACGCAAGCGCAGACAGATAACTATg                                                                                      <  2:696212/246‑1 (MQ=255)
                                                                                                                                                           gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATg                                                                                      <  1:476030/251‑1 (MQ=255)
                                                                                                                                                               cgctgcgctgTTTATGGGTCCGATGCAAGAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGTTTAACCCGGAACGGAAATTTCAGG‑AGGGTGGTGAAGCAAT‑AGTTATCAGAACTATGATTTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAGGCCCATACTCATGGTTATGATGTTGTTGCCTATTGTAATACCGGAACACAAGCGCAGACAGATAACGATGtttt                                                                                  <  2:272447/251‑1 (MQ=255)
                                                                                                                                                               cgcTGCGCTTTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTGCAGG‑AGGGTGGGGAAGCAAT‑AGTGATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGTGTAACGCTTTTGGGCGATTATGCCCATACTCATGGTTATTATTTTGTTGCCTATGGTAATCCCGGAACGCAAGCGCAGACAGATAACGATGGttt                                                                                  <  1:552582/251‑1 (MQ=255)
                                                                                                                                                                 cTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGttttt                                                                                >  1:316868/1‑251 (MQ=255)
                                                                                                                                                                    cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTaa                                                                             <  1:370069/251‑1 (MQ=255)
                                                                                                                                                                        gTTTATTGTTCCGATGCAATAGGCGGGGCGG‑TGAATATCATCACTCCGCG‑CGATGAACCCGGAACGGATATTTCAGG‑AGGGTGTGGAAGCACT‑AGTTATCAGCACTATAGTGTCTTTACGCAGCAACAACTGCGGGATAAGACACGGGTAACGCTGATGGGCGATTATGCCCATAATCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAGCGCACGGGCAGTCAGATAACGATGGTTTTTTAAGTa                                                                          >  2:376330/1‑250 (MQ=255)
                                                                                                                                                                            aTGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAAc                                                                      <  1:699424/250‑1 (MQ=255)
                                                                                                                                                                            aTGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAAc                                                                      <  2:228126/250‑1 (MQ=255)
                                                                                                                                                                                  ccGATGCAATTGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCACATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTAt                                                               <  1:178386/251‑1 (MQ=255)
                                                                                                                                                                                   cGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATg                                                              >  1:980281/1‑251 (MQ=255)
                                                                                                                                                                                      tGCAATAGACGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGTACGGCAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCCACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGTTTATGCCCATACTCATGGGTATGATGTTGTTGCCTATGGGTATACCGGAACCCAAGCGCAGACCACTAACGATGGTTTTTTAAGTAAAAACCTTTATTgcg                                                           >  2:699973/1‑251 (MQ=255)
                                                                                                                                                                                          aTAGGCGGGGGGG‑TGAATATCATCTCGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCt                                                       <  2:805358/251‑1 (MQ=255)
                                                                                                                                                                                              gCGGTGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAg                                                   <  2:633550/251‑1 (MQ=255)
                                                                                                                                                                                               cGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGc                                                  >  1:467962/1‑251 (MQ=255)
                                                                                                                                                                                                        tgAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCCGC‑CGGGTGCGGCAGCAAT‑AGTTATCAGAACTATGATGGCCCGACGGAGCAACCACTGGGGGATAAGACACGCGTAACGCTGTTGGGCGATTAGGCCCATCCTCACGGTTATGATGGTGTTGCCTATGGTAATACCGGTCAGCAAGCGCAGGCAGATATCGATGGTTTTTTAAGTAAAACGCTGTATGGCTCGCTGGAGCATAACttt                                          >  1:473233/1‑251 (MQ=255)
                                                                                                                                                                                                          aaTATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTAc                                        >  1:203657/1‑251 (MQ=255)
                                                                                                                                                                                                           atatCAGCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACt                                       <  1:859528/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTGAACCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCttt                     <  1:828911/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTGAAACCGGAACGGCAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCttt                     >  1:275215/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaCCCGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTtgtg                  >  2:906006/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccGGAACGGAAATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTgcgc               >  2:187821/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aacggaaATTTCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTgcgagtct           <  2:623385/251‑6 (MQ=255)
                                                                                                                                                                                                                                                 ttCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCggctatggctag   <  1:328782/250‑2 (MQ=255)
                                                                                                                                                                                                                                                  tCAGG‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAGGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATTCCGGAACTCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGgctatggctatg  <  1:561686/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
GGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGG‑TGAATATCATCACGACGCG‑CGATGAACCCGGAACGGAAATTTCAGC‑AGGGTGGGGAAGCAAT‑AGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATG  >  NC_000913/4163879‑4164368

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: ‑