Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,844,988 A→T 19.3% intergenic (+64/+55) apt → / ← recR Adenine phosphoribosyltransferase/Recombination protein RecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,844,9880AT19.3% 285.3 / 25.9 173intergenic (+64/+55)apt/recRAdenine phosphoribosyltransferase/Recombination protein RecR
Reads supporting (aligned to +/- strand):  ref base A (45/91);  new base T (20/14);  total (67/106)
Fisher's exact test for biased strand distribution p-value = 9.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

GGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGG  >  ACB143/4844854‑4845123
                                                                                                                                            |                                                                                                                                         
gggcgggTCGGTGGCACGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTAt                                                                                                                                          <  1:387440/130‑1 (MQ=255)
    gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAAAC‑GCGCCTACGCAAG‑CCGC‑GGATGttt                                                                                                                                      >  2:61322/1‑136 (MQ=255)
    gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAAACC‑GCACCTACGAAAG‑CCGC‑GGATGtgt                                                                                                                                      >  2:119748/1‑134 (MQ=255)
       tCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCGCCAACGCAAG‑CCGC‑GGATGtttttt                                                                                                                                   >  4:304887/1‑134 (MQ=39)
        cGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTa                                                                                                                                  <  5:173232/136‑1 (MQ=255)
          gTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAgg                                                                                                                                <  4:138655/136‑1 (MQ=255)
          gTTGCAAGCGGCCGGGGGGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAgg                                                                                                                                <  6:169223/136‑1 (MQ=255)
              aaaGCGGACGGGGTGCCTACATTCTGCCTTACCGAGTTTTCCCTCAGCGCGTATTGATTGGGCATTTGGGGCTGAATTTCAGCACAAT‑‑CG‑CCGGTAAAACC‑CCGCCGAAGAAAA‑CCGC‑GGATTTTTTTTATgcg‑cc                                                                                                                           >  2:263309/2‑135 (MQ=255)
                   ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTAGCG‑GGGGTAAACCC‑‑CTCGAACAAAAG‑CCGC‑GAGTTTTTTTTTGGCg‑gggtgt                                                                                                                       >  4:15870/1‑133 (MQ=255)
                   ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑GCAACGAACGAAT‑CCGC‑GGATGGTTTTTAGGCG‑CCGGgtt                                                                                                                      >  1:378834/1‑136 (MQ=255)
                   ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑AG‑CGGGGTAACCC‑GCGCCCAACGAAG‑CCGC‑GGTTGTGTTTTTGGGG‑CGGGgtt                                                                                                                      >  7:21234/1‑136 (MQ=255)
                                      tGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAAACC‑CCACCAACCAAAA‑CCGCGGTTTGTTTTTTGGGGC‑GGGGTTTACCCCGCAATGGgccccc                                                                                                    >  7:168780/1‑136 (MQ=255)
                                         ctcaCCGAGTTTCCCCTCAGCGAGTATTGCGTGGGCTTTTGGGGCTGCTTTCCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTCCGAACG‑CCGC‑GTATGTTTTGTCGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                <  4:84109/133‑1 (MQ=38)
                                           gACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑GCTGCTACGAAAG‑CCGC‑GGATGTTTTGTTTGCG‑CGGGGTTTCCCGCGCAAGGGCCGCCaacagg                                                                                              >  2:109702/1‑132 (MQ=37)
                                           gACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAAACC‑CCGCCAACCAAAT‑CCGC‑CGATGTTTTTTAAGCG‑CGGGTTTAACCGCCAAATGGCCGCCAACCgg                                                                                              >  3:44475/1‑136 (MQ=255)
                                             ccGAGTTTTCCCTCTCCGAGTATTGAGTGGGCTTGTGGGGCGGCTTGGCAGCCAATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAC‑CCGC‑GTATGTTTTGTAGGCC‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGcc                                                                                            <  2:279331/136‑1 (MQ=37)
                                              cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GGATGTTTTGTAGGCG‑CCGGGTTACCCGCGAATGGGCCGCCCAGCGGccc                                                                                           >  1:252567/1‑136 (MQ=255)
                                              cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑CCGCCAACGAAAG‑CCGC‑GGATGTTTTGTAGGCG‑CGGGGTTACCCGCGAATGGGCCGCCCAGCGGccc                                                                                           >  7:95833/1‑136 (MQ=38)
                                               gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCTCCGACGCAAC‑CCGC‑GGATGGTTTTTAGGCG‑CGGGGTTTCCCGCGAATGGGGCCCCCAGCGGcccc                                                                                          >  5:392701/1‑136 (MQ=37)
                                               gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCGCCAAACCAAG‑CCCC‑CGATGTTTTTTTTGCG‑CGGGTTTTCCCCCCAAATGGCCCCCCACCGGgccc                                                                                          >  5:302619/1‑136 (MQ=16)
                                               gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑CCTCCGACGAAAG‑CCGC‑GTTTGTTTTTTAAGCG‑CGGGGTTAACCGCGAAAGGGCCGCCCAGCGGccc                                                                                           >  1:48927/1‑135 (MQ=37)
                                               gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑CCGCCAACGAAAA‑CCGC‑GGATTTTTTTTATGCG‑CGGGGTTTACCGCGAATTGGCCGCCCAGCGGgcc                                                                                           >  5:363973/1‑135 (MQ=255)
                                               gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGcccc                                                                                          >  2:309734/1‑136 (MQ=255)
                                               gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAAACC‑GCTCCAACCAAAC‑CCGC‑GGATGTTTTTTAGGCG‑CGGGTTTACCCCCGAATGGGCCGCCCACCCGcccc                                                                                          >  8:188238/1‑136 (MQ=37)
                                               gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCACATT‑‑CG‑CGGGTTAACCC‑GCTCCGACGAAAG‑CCGC‑GGATGTTTTGTAGGCG‑CGGGTTTTACCGCGCATGGGGCGCCCAGCCGcccc                                                                                          >  5:91497/1‑136 (MQ=255)
                                                     tCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑CCTCCAACCAAAG‑CCCC‑GGATTTTTTTTACGCC‑CCGGTTTTACCCCCAATGCGGCGCCCAGCCGCCCCAAAGa                                                                                     >  2:220827/1‑135 (MQ=19)
                                                        cTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTTGGCG‑CGGGTTTACCCGCGAAAGGGCCGCCCAGCGGCCCCAATGAAt                                                                                   >  5:303689/1‑134 (MQ=255)
                                                             gtaGTTTTGAGTGGGCTTTGGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAACCCC‑GCTCCTACGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                                           <  6:264861/134‑1 (MQ=38)
                                                             cGATTATTGCGGGGGCTTTTGGGGGTGCTTTGCAGCCCATT‑‑CG‑CGGGTAACCCC‑GAGCCTCCGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                                           <  8:253223/136‑1 (MQ=37)
                                                             cGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCCCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                                           <  7:299185/136‑1 (MQ=255)
                                                              gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑GCGGCTAAGAAAG‑CCGC‑GTATTTTTTTTATGCC‑CGGGGTTTCCCGCGGATTGGCCCCCAAGCGGCCCCCATGAATCAGAGC‑g                                                                           >  8:80066/1‑135 (MQ=255)
                                                              gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑CCTGCGTCGCAAG‑CCGC‑GGTTGTTTTGTAGGCG‑CGGGGTTTCCCGCGAATGGGCCCCCCACCGGCCCCAATGAATAAGcgc‑g                                                                           >  7:35531/1‑135 (MQ=255)
                                                              gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GGATTTTTTGTTGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATAAGAGC‑g                                                                           >  5:183453/1‑135 (MQ=255)
                                                                  aTTGAGTGCGGTTTTGGGGCGTCTTCCCA‑CCAATTG‑CC‑CGGTAAAACCC‑GCTCAAACTAAAG‑CGGC‑TTTTGTTTTGTCGCCG‑GGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATg                                                                      <  4:395567/136‑1 (MQ=255)
                                                                       gggggCTTTTGCGCCTGCATTGAAGACACTT‑‑CG‑CTGCTAACCCC‑CCTCCTACCAAAC‑CCGC‑TTAGGTTTTGTAGCGG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                                 <  6:104991/134‑1 (MQ=19)
                                                                        tGGGGTTTTGGGGCGGCTTTGCTGCCCATT‑‑GG‑GGGGTAAACCC‑GCGCCTACGAACG‑CGGC‑GTATGTTTTGATGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                                 <  5:181469/135‑1 (MQ=37)
                                                                        tGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCCGGGTAAACCC‑‑CGCCATCGCAAGCCCGG‑GTTTGTTTT‑TTGGCGCCGGGTTTGCCCCCCAAGGGGCC‑CCCACCGGGCCCAATTAATACGCGC‑GCAATTGGCCgg                                                                >  8:156446/1‑136 (MQ=255)
                                                                        tGGGCTTTTGGGGCTGCTTTGCAGCCCATG‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                                 <  3:50639/135‑1 (MQ=255)
                                                                        tGGGCTTTTGGGCCTGGTTGGCCGCCATTG‑‑CG‑CGGGTAACCCC‑GCTCATCCCAACC‑CCGC‑TTATGTTTTGTCGGCG‑CGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCgg                                                                <  8:265870/136‑1 (MQ=255)
                                                                        gggggTTTGGGGGCTGGTTTCCCGCCCTTG‑‑CG‑CGGGTAACCCC‑GCTCCTACCAACG‑CCGC‑TTATTTTTTGTAGGGG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                                 <  7:233511/134‑1 (MQ=25)
                                                                        ggggCTTTGGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                                 <  2:263323/134‑1 (MQ=255)
                                                                          ggCTTTTGGGGCTGCTTTCCAGCCAATT‑‑CG‑CGGGAAACCCC‑GCTCCTACCAACC‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGgc                                                               <  4:340065/135‑1 (MQ=38)
                                                                             ttttGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCGCCTACGAAAG‑CCGC‑GGATGTTTTTTAGGCG‑CGGGGTTACCCGCGAATGGGCCGCCCAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGccc                                                           >  7:105491/1‑136 (MQ=38)
                                                                             ttttGGGGCTGCTTTCCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACCAAAC‑CCGC‑GTTTGTTTTGTCGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGccc                                                           <  6:174153/136‑1 (MQ=37)
                                                                              tttGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCGGCTTCGAAAG‑CCGC‑GGTTGTTTTGTAGGCG‑CGGGGTTTCCCGCGACTGGGGCGGCCAGCGGGCCCCAAGGTTAAGTGC‑GGAATGGGGCGGCGcc                                                            >  7:261184/1‑134 (MQ=16)
                                                                              tttGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑GCGCCGTCGAAAC‑CCCC‑GGTTTTTTTGTTGGCG‑CGGGGTTTCCCGCGCATGGGGCGGCAAGCGGGCCCAATGGATAAGCGC‑GGAATTGGCCGGGGcccc                                                          >  6:14700/1‑135 (MQ=255)
                                                                               ttGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGTAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGGATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCt                                                          <  7:7712/135‑1 (MQ=255)
                                                                               ttGGGGCTGCTTGGCCGCCCATT‑‑CG‑CGGTGAACCCC‑GGGCCTACGAAAG‑CGCG‑GTATGTTTTGTAGGCG‑CGGTTTTTCCCGGGAATGGGCCGGCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCt                                                          <  8:192966/135‑1 (MQ=255)
                                                                                tGGGGCGGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GAGCCAACGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTg                                                         <  7:372988/135‑1 (MQ=38)
                                                                                   ggCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCTGCTACGAAAG‑CCGC‑GGTTGTTTTGTAGGCG‑CGGGGTTACCCGCGCATGGGCCGCCCAGCGGCCCCAATGGATCAGAGC‑GATATGGGCCGGCGCCCGGCa                                                       >  3:304530/1‑134 (MQ=37)
                                                                                   ggCTGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑CCGCCGACGGAAG‑GCGCGGTTTGTTTTTTTGGGC‑GGGGTTTACCGCGGATTGGGCCCCCACCCGG‑CCCACAGAATGAGCGCGGAAATGGGCC‑GCGccccggaa                                                      >  4:270630/1‑130 (MQ=255)
                                                                                    ggggctTTGCAGCCCATT‑‑CG‑CGGGAAAACCC‑GCTCATACGAAAG‑CCGC‑GTTTGTTTTGAAGGCG‑CGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  4:353100/133‑1 (MQ=255)
                                                                                    gctgctTTGCAGCCCTTT‑‑CG‑CGGGTAACCCC‑GCTCCTACGAAAC‑CCGC‑GTATGTTTTGTCGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  6:370130/136‑1 (MQ=255)
                                                                                    gctgctTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  1:309734/136‑1 (MQ=255)
                                                                                    gctgctTTGCAGCCCATT‑‑CG‑CGGGAAAACCC‑GCGCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  7:279688/136‑1 (MQ=255)
                                                                                    gctgctTTCCAGCCCTTT‑‑CG‑CGGGTAACCCC‑GCTCCTCCGAAAG‑CCGC‑TTTTGTTTTGTCGGGG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  8:36434/136‑1 (MQ=38)
                                                                                    gcggctTTGCAGCCCATT‑‑CG‑CGGCTAACCCC‑GATCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  1:101953/133‑1 (MQ=255)
                                                                                    gcggctTTGCAGCCAATG‑‑GG‑CGGGAAAACCG‑GCTCCTACGAACG‑CGGC‑GTATGTTTTGTCGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  5:229131/133‑1 (MQ=38)
                                                                                    gcggatttGCAGCCCATT‑‑CG‑CGGTAAACCCC‑GCTCCTCCGAACG‑CCGC‑GTATGTTTTGTCGGCG‑GGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  7:154047/131‑1 (MQ=255)
                                                                                     ctgctTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACTAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAgg                                                   <  1:55620/136‑1 (MQ=255)
                                                                                         tttGCAGCCCATT‑‑CG‑CGGGAAACCCC‑GCTCCAACGAACG‑CCGC‑GTATGTTTTGTCGGCG‑CGGGTTTCCCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                               <  8:121688/136‑1 (MQ=255)
                                                                                          ttggcccccATT‑‑GG‑GGGGTAAACCCGGCTCCACCGAACG‑CCGC‑TTTTGTTTTGTAGGGG‑GGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                               <  3:197687/130‑1 (MQ=255)
                                                                                           tGCAGCCCATT‑‑CG‑CGGGTAACCCC‑GCTCCTACTAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                               <  8:232609/134‑1 (MQ=255)
                                                                                           tGCAGCCCATT‑‑CG‑CGGGTAACCCC‑GCGCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                               <  1:136977/134‑1 (MQ=255)
                                                                                           tGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                               <  8:128026/134‑1 (MQ=255)
                                                                                           tGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGAAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                               <  7:18316/134‑1 (MQ=255)
                                                                                           tGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCGCCTACGAAAG‑CCGC‑GTTTGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                               <  2:8185/134‑1 (MQ=255)
                                                                                           tGCAGCCCATT‑‑CG‑CGGGAAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTCGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATTGGCCGGCGCCCTGCAAAGGCATg                                               <  2:256620/134‑1 (MQ=255)
                                                                                                        tcccgGGTACACCG‑GCTCCACCGACCG‑CCGG‑GTTTGTGTGGTCGGCG‑CGGGTATCCCGGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTa                                     <  4:102381/131‑1 (MQ=255)
                                                                                                        cg‑cgGGTAACCCC‑CCCCCTACTAAAG‑CCGC‑GTATGTTTGGTCGGCG‑CGGGTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccg                                  <  3:229514/136‑1 (MQ=255)
                                                                                                             ggTAACCCC‑GCTCCTAGGAATG‑CCGC‑GTATGTTTTGTAGGCG‑GGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt                              <  5:365665/136‑1 (MQ=255)
                                                                                                             ggTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt                              <  2:250873/136‑1 (MQ=255)
                                                                                                               tAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc                            <  6:357226/136‑1 (MQ=255)
                                                                                                                aaaCCC‑GCGCCTGCGAAAG‑GCGC‑GGATGTTTTGTAGGCG‑CGGGGTTTGCCGCGCATGGGGCGCCAAGCGGCCCCAATGAATCAGAGG‑GAAATGGGCCGGCGCCCTGCAAAGGCAGGGGCCCGGGTGCCGCCCTccc                           >  2:68144/1‑135 (MQ=255)
                                                                                                                            tACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATTAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGc                >  2:65320/1‑136 (MQ=255)
                                                                                                                            tACGAAAG‑CCGC‑GTATGTTTTGTAGGAG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCCCCATCTCCAGc                >  3:170879/1‑136 (MQ=255)
                                                                                                                                aaaG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgc            >  6:70173/1‑136 (MQ=255)
                                                                                                                                 aaG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                                    <  4:115773/95‑1 (MQ=255)
                                                                                                                                 aaG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgcc           <  1:11176/136‑1 (MQ=255)
                                                                                                                                 aaG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgcc           <  6:196789/136‑1 (MQ=255)
                                                                                                                                  accccGC‑TTATTTTTTGGCGGCG‑CGGGTTACCCCGCTAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgcc           <  1:187734/133‑1 (MQ=255)
                                                                                                                                   g‑CCGC‑GTATGTTTTGTAGGCG‑CGGGATTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgccgc         <  3:305955/136‑1 (MQ=255)
                                                                                                                                    cccGC‑GTCTGTTTGGTAGGCG‑CGGGTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgccgc         <  4:188644/135‑1 (MQ=255)
                                                                                                                                      cgc‑gTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCgccgcc        >  3:353566/1‑135 (MQ=255)
                                                                                                                                      cgc‑gTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGccc       <  1:257085/136‑1 (MQ=255)
                                                                                                                                       gc‑gTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCgccgcc        >  8:315012/1‑134 (MQ=255)
                                                                                                                                       gc‑gTATGTTTTGTAGGCG‑CGGGTTTACCCGAGAATGGGCCGCAAAGCGGCCCCAATGCATCAGAGC‑GAAAATGGCCGGCGCCCTGAAAAGGCATGGGCCAGGGTACCTCCGGACACAATATCCAGCTAgccgcc        >  3:324131/1‑134 (MQ=25)
                                                                                                                                        c‑gTATGTTTTGTCGGCG‑CGGGTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAg     <  4:399455/136‑1 (MQ=255)
                                                                                                                                          gTATTTTTTGTCGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt    <  4:193156/136‑1 (MQ=255)
                                                                                                                                          gTATTTTTTGTAGCCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt    <  3:235289/136‑1 (MQ=255)
                                                                                                                                          gTATGTTTTGTCGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt    <  1:243968/136‑1 (MQ=255)
                                                                                                                                          gTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt    <  6:271097/136‑1 (MQ=255)
                                                                                                                                          gTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt    <  6:179699/136‑1 (MQ=255)
                                                                                                                                          gTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt    <  8:34148/136‑1 (MQ=255)
                                                                                                                                          gTATGTTTTGACGGCG‑CGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt    <  5:203202/136‑1 (MQ=255)
                                                                                                                                            aTGTTTTGTAGGCG‑CGGGTTTACCCGCGAATTGGCCCCCAAGCGGCCCCAAAGAATCAGAGC‑GAAAAGGGCCGGCGCCCTGCAAAGGCAAGGGCCACGGTACCGCCGTCCCCCATATCCAGCTCGCCGCCCAGttg  >  8:286093/1‑134 (MQ=255)
                                                                                                                                            aTGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCt                                                          <  3:133490/80‑1 (MQ=255)
                                                                                                                                            aTGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt    <  5:121764/134‑1 (MQ=255)
                                                                                                                                            aTGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTgg  >  7:62877/1‑136 (MQ=255)
                                                                                                                                            aTGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCCCCAAGCGGCCCCAATGAAACAGAGC‑GAAAAGGGCCGGCGCCCTGCAAACGCATGGGCCAGGGGACCGCCGTCCACCATCTCCACCACACCCCCCAGTgg  >  2:383049/1‑136 (MQ=38)
                                                                                                                                            |                                                                                                                                         
GGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGG  >  ACB143/4844854‑4845123

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