Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,805,237 A→C 29.2% intergenic (+164/‑585) dsbG → / → CJ019_04513 Thiol:disulfide interchange protein DsbG/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,805,2370AC29.2% 71.8 / 7.4 98intergenic (+164/‑585)dsbG/CJ019_04513Thiol:disulfide interchange protein DsbG/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/46);  new base C (4/15);  total (29/69)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGTGATTGG  >  ACB143/4805103‑4805368
                                                                                                                                          |                                                                                                                                         
gTGCAGCAGGCGACGTGGGGCAGGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAg                                                                                                                                          <  1:395366/136‑1 (MQ=255)
   aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGG‑C‑g‑g                                                                                                                                    >  6:72325/1‑136 (MQ=255)
   aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑G‑‑GGGG‑C‑g‑g                                                                                                                                    >  5:220371/1‑136 (MQ=255)
   aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑G‑‑GGAC‑C‑g‑g                                                                                                                                    >  6:243655/1‑136 (MQ=255)
    aGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCC‑T‑A‑‑GGCG‑C‑c‑gg                                                                                                                                   >  2:116256/1‑136 (MQ=255)
           gACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑AC‑GGGC‑C‑G‑GGTTTAcc                                                                                                                             >  4:346898/1‑136 (MQ=255)
            aCGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑TGGG‑C‑GCGGGTTAccg                                                                                                                            >  1:14677/1‑135 (MQ=38)
            aCGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGCG‑CGG‑GGTTTAccc                                                                                                                            >  7:89914/1‑136 (MQ=255)
            aCGTGGGGCATGGCAAATTATACGCCGGTGTTCGCGGGTGACCCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGAGGTCGGCCG‑T‑AC‑GGGG‑G‑G‑GGTTTGccc                                                                                                                            >  3:303070/1‑136 (MQ=255)
                ggggCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCC‑T‑A‑‑GGGG‑CGG‑GGTTTACCcccg                                                                                                                         >  1:384078/1‑135 (MQ=38)
                      tGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGG‑CGG‑GGGTTACCGCCGAAGGGGc                                                                                                                  >  8:304955/1‑136 (MQ=37)
                          aCCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGCGCC‑G‑GGTTTACCCCCGCATGGCCCgcc                                                                                                              >  6:162048/1‑134 (MQ=38)
                               cTACGCCGGTGTTCGCTGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑AGTGGTG‑G‑G‑GGTTTACCCCGGACAGGGGCgccgagg                                                                                                          >  4:234537/1‑130 (MQ=255)
                               cTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑G‑‑GGGG‑C‑G‑GGTTT‑CCCGCCCAAgggccgg‑gccaggg                                                                                                       >  6:189069/1‑132 (MQ=37)
                               cTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCCGCGCGGTCGGCCG‑TGG‑‑GGGG‑C‑G‑GGGTTGCCCCCCAAGGGCCCggccagcc                                                                                                         >  6:118411/1‑130 (MQ=37)
                                  cGCCGGTGTTCGCGGGGTAACCCGCTCCTACGGGGCTCGCGCAAACGGCTGGGTTTTGCGCCCGGCCCGTGAGCGCACGGCCGGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGc                                                                                                   <  5:293775/136‑1 (MQ=37)
                                          ttCGCGGGGGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTc                                                                                            <  7:379339/135‑1 (MQ=255)
                                              cgGGTTAACCCGCTCCTACGGGGCGCGCGCTAACGGCGGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTccacca                                                                                       <  7:195932/136‑1 (MQ=255)
                                                 gttaccccGCTCCTACGGGGCGCGGGATACCGGGTGGGTTTTGCGCCCGGCCCGTTATGGCACGGCGTGCGGGGTGGGCCT‑T‑T‑‑GGTG‑G‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGcc                                                                                    <  7:73580/131‑1 (MQ=16)
                                                 gttaaCCCGCTCCTACGGGGCGCGGGCTAAGGGGTGGGTTTTGCGCCTGGCCCGTTATCGCACGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGTG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGcc                                                                                    <  6:354757/133‑1 (MQ=21)
                                                             gtACGGGGCGCGCGCTAAGGGGGGGGTTTTGCGCCCGGCCCGTGATCGCAAGCCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgc                                                                        <  2:21900/135‑1 (MQ=21)
                                                             cTAGGGGGCTCGCGCTACCGGCTGGGTTTGGCGCCAGGCCCGTGAGCGCACGGCCGGCGCGGGCGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgc                                                                        <  8:261215/136‑1 (MQ=21)
                                                                gggggCGCGCGCTACGGGCTGGGTTTTGCGCCTGCCCGGTTAGCGCACGGCGGGCGGGGTCGGCGG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgc                                                                     <  3:355377/135‑1 (MQ=18)
                                                                  gggccggcgcTATCGGGTGGGGTTTGCGCCTGGCCCGTGATCGCATGGCCGGGGGGGCGGGCCT‑T‑A‑‑GTTG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  7:230265/130‑1 (MQ=17)
                                                                  gggCTCGGGCAAACGGCGGGGTTTTGCGCCAGGCCGGTCATCGCACGGCCTGCGCGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  6:181743/136‑1 (MQ=18)
                                                                  gggCTCGCGCTACCGGCGGGGTTTTGCGCCAGGCCGGTGATCGCAAGGCCTGCGCGGGCGGCCG‑T‑A‑‑GGCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  7:259190/136‑1 (MQ=21)
                                                                  gggCTCGCGCTACCGGCGGGGTTTTGCGCCAGGCCGGTGATCCCACGGCCTGCCGGGTCGGCCG‑T‑A‑‑GGAG‑G‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  6:438/136‑1 (MQ=21)
                                                                  gggCGCGCGCTAACGGGTGGGTTTTGCGCCAGGCCCGTGATCGCACGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  7:371438/136‑1 (MQ=255)
                                                                  gggCGCGCCCTAAGGGCGGGGTTTTGCGCCAGGCCCGGTAGCGCACGGCGGGCGGGGTGGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  8:43353/136‑1 (MQ=18)
                                                                   ggCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑G‑‑GGGGCC‑G‑GGGTTACCC‑CGCAAGGGGCGG‑GCCCGCGGGGCCAATTAACACCACCCGTTCCCCGGTgcccccc                                                                   >  5:111645/1‑135 (MQ=255)
                                                                    gcgcgcgCTAACGCCTGGGGTTGGCGCCAGGCCCGTGATCGCAAGGCCGGCGCGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccgcc                                                                 <  7:128954/133‑1 (MQ=21)
                                                                    gCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGGG‑C‑G‑GGGTTAACCGCGCAAGGGCCGG‑GCCGGC‑GGGGCCAATTCAACCAGCCGGTCACCCGTgccgccccc                                                                 >  1:251729/1‑136 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑G‑‑GGGG‑G‑G‑GGTTTTCCCCGGAAGGGGCCGC‑CCGGGG‑GGGCCAATTTGC‑CCAGCCGGGGACCGGTggcgccggccccg                                                             >  5:80154/1‑130 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGGCC‑G‑GGTTTACCGCCCAACGGGCGGG‑CCGGGG‑GGGCCAATTACC‑CCACCCGGGAACCCTTgccgcccccccag                                                             >  5:256827/1‑132 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGGCC‑G‑GGTTTACCC‑CGCAATGGGCGG‑CCCAGG‑TGGCCAATTTACCCCACCCGTTCAACGGTgccgcccccccag                                                             >  5:297717/1‑132 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCCCCCAAGGGGCGG‑ACCTGC‑TGGCCCATTTACCCCCGCCGTTGCCCGTTGCCCCCCCCCCTg                                                             >  8:273346/1‑136 (MQ=21)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGGTTACCCCCGAAAGGGCCCG‑CCCAGC‑TGGCCCATTTACCCCAGCCGGTCCCCGGTTCCCCCGCCCCTg                                                             >  8:95058/1‑136 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGGTTACCCCCCAAATGGGCGG‑CCCAGC‑GGGCCAATTTCCCCCCGCCGGTCCCCCTTgccgccgccccag                                                             >  7:307852/1‑132 (MQ=21)
                                                                          gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑TGAG‑CGG‑GGTTTACCC‑CGCAAGGGGCGG‑CCCATC‑GGGGCAATTTTCACCCGGCGGTGACCGGTggcgccgccccagtt                                                           >  7:220033/1‑130 (MQ=255)
                                                                          gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGGCC‑G‑GGTTTACCCCGCAAGGGGCCGG‑CCGGGG‑GGGCCAATTTCCCCCGCCCGTGACCCGTT‑CCGCCGCCCCAGTa                                                           >  2:211403/1‑136 (MQ=255)
                                                                          gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGCG‑CGG‑GGTTTACCCCGCAAGGGGCCGG‑CCCGCG‑GGGCCAATTGCC‑CCCGCCGGTGACCGGTgccgccgccccagtt                                                           >  4:304833/1‑130 (MQ=255)
                                                                          gcTAACGGCGGGGTTTTGCGCCAGGCCGGTTATCGCACGGCCGGCGCGGTCGGCGG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                           <  2:47952/136‑1 (MQ=18)
                                                                          gacagcGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                           <  2:322679/131‑1 (MQ=255)
                                                                           aaaaCGGCTGGGTTTTGCGCCAGGCCCGTGATCGCACGGCCGGCGGGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                          <  2:58528/134‑1 (MQ=21)
                                                                                     tgTTTTGCGCCAGGCCGTTTTGCGAAAGGCCTGCGCGGTGGGCCG‑T‑A‑‑GGTG‑G‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  6:85619/135‑1 (MQ=14)
                                                                                     ggTTTTGGGCCCGGCCCGTGATCCCAAGCCCTGGGCGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  7:40322/136‑1 (MQ=18)
                                                                                     ggTTTTGCGCCAGGCCGGTGAGCGCAAGGCCGGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  4:287217/136‑1 (MQ=255)
                                                                                     ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  8:89059/136‑1 (MQ=255)
                                                                                     ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  8:312986/136‑1 (MQ=255)
                                                                                     ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCGGCGGGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  2:384078/136‑1 (MQ=255)
                                                                                     ggTTTTGCCCCCGGCAGTTGAGCCCAAGCCCTGCGCGGTGGGCCT‑T‑A‑‑GGAG‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  5:219933/136‑1 (MQ=14)
                                                                                     ggTTTTCCGACCGGACCGTGATCGCACGGCCTGCGCGGTCGACCG‑T‑A‑‑GGGG‑C‑G‑GTTTTACCAGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCATCCGTTCACCGTTGACGCCGCCACTGTATTCGC‑GGTTaa                                               <  8:113411/136‑1 (MQ=17)
                                                                                     ggTTTGGCGCCAGGCCCGTGAGCGCAGGCCCTGCGGGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  6:93660/136‑1 (MQ=18)
                                                                                     ggTTTGGCCCCTGCCCCGTTATCCCACGGCCGGCGCGGTGGGCCT‑T‑A‑‑GGCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  2:307684/136‑1 (MQ=17)
                                                                                     ggTTTGCCGCCCGCCCGGTGACCGCACGGCCGGCGGGGTCGGCGG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                               <  6:314782/136‑1 (MQ=17)
                                                                                      gTTTTGCGCCTGGCCCGTGATCGCATGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑G‑G‑GGTTTACCCGCGAAAGGGCGGG‑ACGGGC‑TGGCAAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaaa                                              <  7:188788/136‑1 (MQ=18)
                                                                                      gTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑C‑GGTTTTCCCCCGAAAGGGCCCG‑GCCGGC‑GGGCCAATTTACACCCGCCCTTGACCCGTGCCGCCGCCACAGTATTAGC‑GGGGaa                                               >  8:312917/1‑135 (MQ=255)
                                                                                         ttGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGG‑CGG‑GGTTTACCC‑CGCAAGGGGCGGCCCCGGG‑GGGCCAATTGCC‑CCCGCCGGTGACCGGTTGCGCCGCCCCTGTTTTAGCGGGGTAAAcc                                            >  2:241283/1‑136 (MQ=255)
                                                                                           gccccAGGCCGTTTATCCCACGGCCGGCGGGGTCGCCTT‑T‑T‑‑GCCG‑C‑G‑GTTTAACCGGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCg                                          <  7:209060/132‑1 (MQ=14)
                                                                                             gcCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                           <  7:153932/117‑1 (MQ=255)
                                                                                             gcCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                           >  8:153932/1‑117 (MQ=255)
                                                                                               ccggCCGGTGATCTCACGGCCGGCGGGGCCGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                                     <  3:329533/134‑1 (MQ=17)
                                                                                               ccggCCCGTGATCGCAAGGCCGGCGCGGTCGGCGG‑T‑C‑‑GGAG‑C‑G‑GTTTAACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                                     <  2:372510/134‑1 (MQ=17)
                                                                                               cAGGCCCGTGATCGCCCGGCCTGCGCGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                                     <  4:306195/136‑1 (MQ=255)
                                                                                               cAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                                     <  2:292270/136‑1 (MQ=255)
                                                                                               cAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                                     <  8:106918/136‑1 (MQ=255)
                                                                                               cAGGCCCGTGACCGCAAGCCCGCCGCGGGCGGCGT‑T‑C‑‑GCCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGATCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                                     <  5:100578/136‑1 (MQ=17)
                                                                                               acgccccGTGACCGCACGCCCTGCGGGGCCGGCCT‑T‑A‑‑GGAG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                                     <  5:113202/132‑1 (MQ=17)
                                                                                                  gCCCGTGATCGCACGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACg                                  <  8:195716/136‑1 (MQ=255)
                                                                                                    cggttaTCCCACGGCCTGCGCGTGCGGCGT‑T‑A‑‑GGAT‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCTTTCACCGTTGCCTCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCt                                <  8:157855/131‑1 (MQ=14)
                                                                                                    ccGTGATCGCACGGCCGGCGCGGTCGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCt                                <  7:75202/136‑1 (MQ=255)
                                                                                                    ccGTCAGCGCACGGCCGGCGCGGTGGGCGG‑T‑A‑‑GCAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCt                                <  3:278435/136‑1 (MQ=17)
                                                                                                      gTGATCGCACGGCGGGCGGGGTCGGCCG‑T‑A‑‑GTAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGa                              <  4:267657/136‑1 (MQ=17)
                                                                                                         agcGCACGGCCGGCGCGGGCGCCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCg                           <  4:294835/134‑1 (MQ=14)
                                                                                                         aTCGCACGGCCGGCGGGGGCGGCGG‑T‑A‑‑GCCG‑G‑G‑GTTTAACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATc                            <  2:89742/135‑1 (MQ=17)
                                                                                                          tCGCATGGCGGGCGGGGGGGGCGG‑T‑T‑‑GGCG‑G‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCg                           <  5:196628/135‑1 (MQ=14)
                                                                                                          gcGCACGGCCGGCGGGGGCGGCCT‑T‑T‑‑GGAG‑G‑G‑GGTTAACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCg                           <  6:122572/134‑1 (MQ=17)
                                                                                                           cGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCgg                          <  6:358066/135‑1 (MQ=255)
                                                                                                           cGCAAGGCCTGCGCGGGCGGCGG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCgg                          <  5:179880/135‑1 (MQ=255)
                                                                                                             catggCCGGCGGGGTGGGCCG‑T‑T‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCggggg                       <  2:6294/133‑1 (MQ=17)
                                                                                                             cacggCCGGCGGGGGCGGCGG‑A‑A‑‑GTTG‑G‑G‑GGTTACCCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCggggg                       <  4:41124/133‑1 (MQ=17)
                                                                                                             cAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCggggg                       <  1:233371/136‑1 (MQ=255)
                                                                                                               tggCCTGCGGGGTCGGCCG‑T‑A‑‑GGTG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTg                     <  2:245645/135‑1 (MQ=255)
                                                                                                               cgggcTGCGCGGTGGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTg                     <  5:127484/132‑1 (MQ=255)
                                                                                                               aGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTg                     <  1:248870/136‑1 (MQ=255)
                                                                                                                ggCCTGCGCGGGCGGCCT‑T‑T‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTgg                    <  6:161160/136‑1 (MQ=255)
                                                                                                                gcccTGCGGGGTCGGCCG‑T‑A‑‑GCAG‑G‑G‑GTTTAACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTgg                    <  1:214868/134‑1 (MQ=17)
                                                                                                                     gggcgGTCGGCCG‑T‑A‑‑GGCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGaca               <  3:246591/134‑1 (MQ=255)
                                                                                                                      cgcgGACGGCCG‑T‑A‑‑GTAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGaca               <  5:68428/135‑1 (MQ=255)
                                                                                                                        gggTGGGCGG‑T‑A‑‑GCCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt            <  2:245430/135‑1 (MQ=255)
                                                                                                                        gggTCGGCGG‑A‑A‑‑GTCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt            <  2:75684/135‑1 (MQ=17)
                                                                                                                        gggTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt            <  2:111005/135‑1 (MQ=255)
                                                                                                                        cgGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt            <  8:33648/136‑1 (MQ=255)
                                                                                                                         ggTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                          <  8:374611/90‑1 (MQ=255)
                                                                                                                         ggTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                          >  7:374611/1‑90 (MQ=255)
                                                                                                                            gggCGG‑T‑A‑‑GGTG‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGTGTGGAGACACATTggg        <  2:100838/135‑2 (MQ=17)
                                                                                                                                    gtttggAG‑G‑G‑GTTTAACCCGGGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGTGAtt    <  2:396059/132‑1 (MQ=17)
                                                                                                                                       tggAG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGTGATTgg  <  8:4488/134‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGTGATTGG  >  ACB143/4805103‑4805368

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