Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,502,7610GC26.5% 103.0 / 13.3 49A224G (GCG→GGG) insGIS4 family putative transposase
Reads supporting (aligned to +/- strand):  ref base G (11/25);  new base C (13/0);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 9.51e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTGACCTGCGGGTAGAGCGCCGGGTT‑CCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTG  >  NC_000913/4502527‑4503007
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
tcccctccGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCCGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCAt                                                                                                                                                                                                                                                <  2:448685/247‑1 (MQ=255)
  tcctccGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑C‑GGCTCCCGGCCtttt                                                                                                                                                                                                                                             >  1:499940/1‑250 (MQ=255)
     tccGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGCGGTATCATCCAGTGACGGTG‑GTCTCCCCCC‑A‑GGccccctgccattgac                                                                                                                                                                                                                                           >  2:207742/1‑236 (MQ=255)
      ccGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑c                                                                                                                                                                                                                                          <  1:173861/249‑1 (MQ=255)
       cGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGGG‑TTTTCCCCCC‑C‑GGCTccccggacttaata                                                                                                                                                                                                                                          >  1:147261/1‑237 (MQ=255)
         gggaaggctatgtcggccaccaacaaCGTCAGCAGATGTCAGCCTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCCGTCACTTCTTTCCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTGTCATCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTcc                                                                                                                                                                                                                                     <  2:147261/226‑1 (MQ=255)
                             cATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGGTTTTCCCCCC‑G‑GGC‑CCAGGCATTTAA‑AAGTCCCCGTGAG‑GAATAAACttttt                                                                                                                                                                                                                >  1:86367/1‑249 (MQ=255)
                             cATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTAt                                                                                                                                                                                                                <  1:435765/251‑1 (MQ=255)
                                cGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCa                                                                                                                                                                                                             <  1:145802/251‑1 (MQ=255)
                                   cGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCCCC‑ACGGC‑CCAGGCACTTTA‑CAGTCCCAGTTAG‑TAATAAACCTTATCCAtt                                                                                                                                                                                                           >  1:535040/1‑250 (MQ=255)
                                    gTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCCCC‑C‑GGCTCCCAGCATTTAA‑CAACCCCCGTGGG‑TAATAACCttttttcattc                                                                                                                                                                                                          >  1:385723/1‑242 (MQ=255)
                                       agcagATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGt                                                                                                                                                                                                      <  1:48178/251‑1 (MQ=255)
                                           gATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTGCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTcg                                                                                                                                                                                                                                     <  1:212875/217‑2 (MQ=255)
                                                 aGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCCCC‑C‑GGCTCCAAGCATTTAA‑AAGTCCCCGTGGG‑TAATAACCTTTTTCCATTAACGTCAGAGGGttt                                                                                                                                                                                            >  1:424090/1‑250 (MQ=255)
                                                         ccTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCgg                                                                                                                                                                                     <  1:504201/250‑1 (MQ=255)
                                                                gcgTCCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCt                                                                                                                                                                            <  1:155293/251‑1 (MQ=255)
                                                                      cGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCCCCCA‑GGC‑CCAGGCCATTAA‑CAGACCCCGTCGG‑ACATAACCTTTTTCCATTTAAGTCAGGGGGTTATAGCCGG‑TGTTTTCTATTAGt                                                                                                                                                                      >  1:150327/1‑251 (MQ=255)
                                                                                gCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCaa                                                                                                                                                            <  1:29523/251‑1 (MQ=255)
                                                                                      aGTCACTTCATTTCCCCGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTTCCCCC‑A‑GGCTCCCGGCATTTTA‑CAATCCCCGAGAG‑GAATAACCTTTTTCCCTTAACGTCAGAGGGTTATCGCCGG‑GTTTTTTTTTAAGTTTCTCAGCAAACTCAt                                                                                                                                                      >  1:504414/1‑251 (MQ=255)
                                                                                          aCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTTCCCCCC‑C‑GGCTCCCGGCATTTAA‑CAATCCCAGTGAG‑TAATAACCTTTATTCATTTAAGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCCAACTCAttttt                                                                                                                                                  >  1:94312/1‑250 (MQ=255)
                                                                                                   cccAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGttct                                                                                                                                          <  1:500059/250‑1 (MQ=255)
                                                                                                            ggCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTTCCCCCC‑A‑GGCCCCAGGGCTTTAA‑CAGTCCCCGTGAG‑GAATAAACTTTATTCATTAACGTCAGAGTGTTTTCGCCGG‑GTTTTTCTTTAAGTTGCTCAGCAAACTCATTTTTCCTGTTCTTCATC‑Gttcc                                                                                                                               >  1:28897/1‑248 (MQ=255)
                                                                                                               cacttttgtacggccgggtggCTGTTTTTCACCTTCACCATATCATCGCCTTTCCCCATTTTTCTTATCTCTTCGTATTGCGCGCCCTTTCTTAGAGGTATCATCCAATGGCGGTG‑TTTTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCCTTCTCCATTAACGTCAGAGTGTTAGCGCCGG‑TTTGTTCTTTAAGTTGCTCAGCAAGCTCATTTTCTCTGTTCTTCATC‑GTGCCGaa                                                                                                                            <  2:182398/232‑1 (MQ=255)
                                                                                                                        cGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGc                                                                                                                   <  1:215249/251‑1 (MQ=255)
                                                                                                                               gCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTca                                                                                                              <  1:371899/249‑1 (MQ=255)
                                                                                                                                  ggCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTTCCCCC‑A‑GGCTCCCGGCCTTTAA‑CAGTCCCCGTGAG‑TAATAACCTTTATTCATTTACGTCAGGGAGTTATCGCCGG‑TTTTTTCTTTTAATTGCTCAGCCAACTCCTTTTTTCCGTTCTTCATC‑GTGCCGAAGGCttccgccggcatc                                                                                                            >  1:358017/1‑236 (MQ=255)
                                                                                                                                   gCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGAt                                                                                                        <  1:67504/251‑1 (MQ=255)
                                                                                                                                    cTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑GTTTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGGCCCCGTGAG‑GAAAAACCTTTTTCCATTAACGTCAGAGTGTTATTGCCGG‑GTTGTTCTATAAGTTGCTCAGCAAGCTCCTTTTTGCTTTTCTTCATC‑GTGCCCAAGGccgccaccgtccgccaatg                                                                                                       >  1:102770/1‑233 (MQ=255)
                                                                                                                                      ggTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATgg                                                                                                      <  1:496108/250‑1 (MQ=255)
                                                                                                                                             cgCTTCCCCAGATGATCCCCTTTACCCAGTTTTCGGATCTCTTCTTATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCAGTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAg                                                                                             <  2:253602/250‑1 (MQ=255)
                                                                                                                                                  cACCAGCTGAGCGCCGTTACCCAGTTTTCTGATCTCGTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCCCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTGCCa                                                                                        <  2:188860/251‑1 (MQ=255)
                                                                                                                                                        atgatCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCATTTCTGAGAGGTATCATCCAGTGGCGGTG‑GTTTACCCCC‑A‑GGGCCCCGGGATTTTA‑AAATCCCAGGGGGGAAAAAACCTTT‑TCCATTAAAGACAGGATGTTAACACCGGGTTTGTT‑TATAAATTGCTAAACAAACTCATTTTCGCTGATCTTAATCAGGGCCGAAGAC‑GCAGCCCTCAACAGATAGC‑TGGTCAGTTcgctcttgc                                                                                  >  1:250699/1‑243 (MQ=255)
                                                                                                                                                         tgatCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCCGGCATTTTA‑CAGTCCCCGGGAG‑GAATAACCTTTATCCATTAAAGTCAGAGTGTTATTGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTTGCTGTTTTTTATC‑GTGCCGCAGGCTGCAGCCGGCAGCAGGTGGC‑TGGTCAGTTCCATCttgcc                                                                                 >  1:122160/1‑247 (MQ=255)
                                                                                                                                                            tCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑C‑GGCTCCAGGCATTTAA‑CAGTCCCCGGGAG‑GAATAAACTTTATCCATTAAAGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTTTTCATC‑GTGCCGAAGGCTGACGCCGTCAGCAAATGGC‑TGGGCAGTTTCCTCTGGCAGAc                                                                              >  1:405773/1‑251 (MQ=255)
                                                                                                                                                             cGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTA‑TTTTTCCGCCCA‑GG‑TCCCGGCATTTAAACA‑TCCCAGTGAG‑TAATAACCTTTATCCATTAAAGTCAGAGGGTTAATGCCGG‑GTTTTTCTTTAAATTGATCAGCCAGCTCATTTTCGCTGGTCTTCATC‑GTGCCGAAAGATGCAGCCGGTAGCAGATGGCGTGG‑CACTTTCATCTGGCCGAcc                                                                             >  1:24822/1‑251 (MQ=255)
                                                                                                                                                                 tttACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGGACCAGGCATTTAA‑AAAAACCGATGAG‑GAAAAAACGTTATCCATTAAAGACAGAGAGTTATCGCCGG‑GTTGGTCTATAAGGTGGCCAGAAAAACCAGATTCGCTGTTTTTCAAC‑GTGCCGAAAGGGGCAGCCGTCAGCAGAAGGG‑AGGTCAGTTTCAACTGGCAggacactt                                                                         >  1:319493/1‑244 (MQ=255)
                                                                                                                                                                    aCCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTg                                                                      <  1:102130/250‑1 (MQ=255)
                                                                                                                                                                              cTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑C‑GGCTCCCGGCCTTTAA‑AAATCCCAGGGAG‑GAATAACCTTTTTCCCTTAAAGTCAGAGTGTTATCCCCGG‑TTTGTTTTTTAAGTTGCTCAGCCAGCTCATTTTTGCTTGTCTTCCAC‑GTTCCCCAGGCTGCAGCCGTCAGCAAATGGC‑TGGTCAGTTTCCTCTTGCCGACCA‑TTTTTACCTGCggggg                                                           >  1:570605/1‑249 (MQ=255)
                                                                                                                                                                                gATCCCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTGACCTGCGGGTaga                                                         <  2:122160/251‑1 (MQ=255)
                                                                                                                                                                                   cgcttcaatttGCGCTCCCTTTCGTAGAGGTAGCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑TGCTCCAGGCATTGAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAATGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGTCCA‑TTTTGACCTGCGGGTAGAgc                                                       <  2:302435/241‑1 (MQ=255)
                                                                                                                                                                                     ctTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTGACCTGCGGGTagag                                                        <  1:488532/247‑1 (MQ=255)
                                                                                                                                                                                                   cccTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCCCC‑C‑GGCTCCAGGCATTTAA‑CAGTCCCCGTGAG‑TAATAACCTTTATCCATTAAAGTCAGAGTGTTATCGCCGG‑GTTGTTCTTTAAGTTTCTCAGCAAGCTCATTTTTGCCGTTCTTCATC‑GTGCCGAAAGCTGCAGCCCTCAGCAGATGGC‑TGGTCAGTTCCCTCTTGCCGACCC‑TTTTTACCTGCGGGTAGAGCGCCGGGTTTCCCggatg                                      >  1:83291/1‑248 (MQ=255)
                                                                                                                                                                                                     cTTTCTGAGAGGTATCATCCAGTGGCGGTG‑GTCTCCCGCC‑A‑GGCTCCCGGCATTTTA‑CAATCCCCGGGAG‑GAAAAACCTTTTTCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCCGCAAGCTCATTTTTGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATAGC‑TGGTCAGTTCCATCTGGCAGACCAATTTTG‑CCTGCGGGGAGGACGCCGGGGT‑CCCGGCAtgtg                                   >  1:399589/1‑251 (MQ=255)
                                                                                                                                                                                                           gagaGGTATCATCCAGTGGCGGTT‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTGACCTGCGGGTAGAGCGCCGGGTT‑CCCGGCATGTGTCt                                <  1:192390/248‑1 (MQ=255)
                                                                                                                                                                                                             gagGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑C‑GGGTCCCGGCATTTTA‑CAGTCCCCGTGAG‑GAAAAACCTTTATCCATTAACGTCAGAAGGTTATCGCCGG‑TTTGTTCTATAAATTGCTCAGCCAGCTCATTTTTGCTGTTCTTCATC‑GTGCCGAAAGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTTACCTGCGGGGAGAGCGCCGGGTT‑CCCGGCCTGTGTCTGGCgg                           >  1:444580/1‑251 (MQ=255)
                                                                                                                                                                                                                  gtcgtcCAGGGGCGGCG‑TTCCCCCGCC‑A‑GGATCCAGGCATGTAA‑CAGTCCCAGTGAG‑TAATAACCTGGATCCATTAACGGCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCGCAGCAAGCTCATTTGCGCTGTGCTGCATC‑GGGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGTCAGACCG‑TTTTGACCTGCGGGGAGAGCGCCGGGTT‑CCTGGCATGTGTCTGTCGGGGGaa                      <  2:295574/247‑1 (MQ=255)
                                                                                                                                                                                                                             gCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTGACCTGCGGGTAGAGCGCCGGGTT‑CCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGt           >  1:534915/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ctcCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGACGACCA‑TTTTGACCTGCGGGTAGAGCGCCGGGTT‑CCCGGCATGTGTCTGGAGGGGGAAGGCTGAATCGTTCTCTGgt   >  2:467808/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cccGCC‑A‑GGCTTCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCATCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTGACCTGCGGGTAGAGCGCCGGGTT‑CCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGgtg  <  2:311557/250‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
TCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCC‑A‑GGCTCCAGGCATTTAA‑CAGTCCCAGTGAG‑TAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGG‑TTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATC‑GTGCCGAAGGCTGCAGCCGTCAGCAGATGGC‑TGGTCAGTTCCATCTGGCAGACCA‑TTTTGACCTGCGGGTAGAGCGCCGGGTT‑CCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTG  >  NC_000913/4502527‑4503007

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