Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,697,739 T→G 18.6% V428G (GTG→GGG)  bcsE → cellulose production protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,697,7390TG18.6% 177.6 / 4.0 89V428G (GTG→GGG) bcsEcellulose production protein
Reads supporting (aligned to +/- strand):  ref base T (17/50);  new base G (18/4);  total (35/54)
Fisher's exact test for biased strand distribution p-value = 5.76e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCACTACCTTGCTGTCAATGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACG  >  NC_000913/3697506‑3697976
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
tatCACTACCTTGCTGTCAATGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGTTGAATAACCCTCTATTACCTGCCCACGGTAAAGGAGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCGGTG‑TCGCCCTAACCGTACCGGCGAAATCATGACCATTGGG‑GGGAAGCGGCAGGGG‑gtgt                                                                                                                                                                                                                                                 >  2:453509/1‑239 (MQ=255)
           tGCTGTCAATGACCCAGCCGCTCAAACTGCGTGGTTTCCAGACGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGCAGCTGTG‑TCGCCCTAACCGTACCGGCGAGATACTGACCATTGCC‑GGTAATCGCCTGGGG‑CAGtttacaccact                                                                                                                                                                                                                                       >  1:359822/1‑230 (MQ=255)
           tGCTGTCAATGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGACCAACCCTGGACCCTGTG‑TCCCCCTAACCGAACCGCCGTAATCATGACCATTGCC‑GGCAATAGGAGTGTG‑GTGttttgtccaatc                                                                                                                                                                                                                                      >  1:640583/1‑230 (MQ=255)
           tGCTGTCAATGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGGATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATAACCATTGGC‑GGGAAACGGCGCGGG‑CGGTTTCTCCCAttt                                                                                                                                                                                                                                      >  1:506195/1‑238 (MQ=255)
            gCTGTCAATGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCCCAtttt                                                                                                                                                                                                                                     >  2:406743/1‑237 (MQ=255)
             ttgttaattccccttcctcttaaaCTGGGTGGTTTCCAGAAGTGGGATTTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTTCCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCAGTACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTg                                                                                                                                                                                                                                    <  1:179010/218‑1 (MQ=255)
             cTGTCAATGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCCGGGG‑GGGTTTTTTTCATTTTg                                                                                                                                                                                                                                    >  2:723671/1‑239 (MQ=255)
                   aTGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGGTGGGG‑CTGTGTTTCTCttttttttgggt                                                                                                                                                                                                                              >  2:624920/1‑230 (MQ=255)
                   aTGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGGGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGAt                                                                                                                                                                                                                              <  2:980973/239‑1 (MQ=255)
                     tgcccAGCCCTTCAATCTGCGTGGTTTCCAGAATTGGGATGGGTTCTTGAATGCCGTCAACAACATGATGAATAACCCTCTATTACCGGCCCCCGGTAAAGGGCTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTa                                                                                                                                                                                                                            <  1:255120/237‑1 (MQ=255)
                       ccctcccGCTAAAACTGCGTTGTTTCCAGAAGGGGGATTTGTTCTGTAAGGCCGGCAAAAACATGATGAATAACCCTCTATTCCCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAAc                                                                                                                                                                                                                          <  1:935865/234‑1 (MQ=255)
                       cccAGCCGCTCAAACGGCGTGGTTTCCAGAAGTGCGGTGTTTTCTGGACGGCCGTCACAAACAGTATGATCAACCCTCTATTACCTGCCCATGAGAAAGGCGTTCTGGCTCCCCGACGTCCGGTTCCGGGTATCCGCGTTGAACAGGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGCATTAAc                                                                                                                                                                                                                          <  1:801869/239‑1 (MQ=255)
                        ccAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAAc                                                                                                                                                                                                                          <  2:345184/238‑1 (MQ=255)
                          aGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATTTCATGACCATTGGC‑GGTAATCCGCCGGGG‑CTGttttttctccttctttggggtaaaga                                                                                                                                                                                                                        >  1:685452/1‑215 (MQ=255)
                          aGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGGAATCGGCTGGGG‑CTGGTTCTTTCCTTCTTTCGGTTTAACGAt                                                                                                                                                                                                                       >  2:986372/1‑239 (MQ=255)
                          aGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTCCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGAAGCTGTG‑TCGCCCTAACCGTACCGGCGGTATCCATACCATTGGCGGGTAACCGGCGGGGG‑GTGttctttttcttttaacaacccggag                                                                                                                                                                                                                         >  1:115539/1‑215 (MQ=255)
                             cGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGAt                                                                                                                                                                                                                       <  2:685452/236‑1 (MQ=255)
                               cTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGGAATCGGGCGGGG‑CGGGTTTTCTCCTTTTTTCGGGTTTATGCTCTggg                                                                                                                                                                                                                  >  1:345184/1‑238 (MQ=255)
                                      tGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑cgcg                                                                                                                                                                                                          <  2:164744/239‑1 (MQ=255)
                                        cGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTATTCGTCTGGGG‑CTTGTTCTCTCATTTTGTTGGATTTTCGATCCGGATAc‑ctc                                                                                                                                                                                                           >  2:845526/1‑234 (MQ=255)
                                         gTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCCGCTGGTG‑CTGTTTTCTTCATTTTGTCGGGTTAACGATATGGGTAc‑ccggttt                                                                                                                                                                                                       >  2:744723/1‑233 (MQ=255)
                                            gTTTCCAGAAGGGGGATGTGGTCTGTAATGCCGTCAACAACATGATGAATTACCCTCTATTTCCTGCCCAGGGTTAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAt                                                                                                                                                                                                    <  2:647812/239‑1 (MQ=255)
                                                     aGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCGCGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTc                                                                                                                                                                                          <  2:168498/239‑1 (MQ=255)
                                                       tGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGAACCGGCGAAATCATGACCATTGGC‑GGGAATCGGCTGGGG‑CTGGTTTTCTCATTCTTTTGGATTTAAAAATTCGATAC‑CGCCGTGAAT‑AACATTTTtcc                                                                                                                                                                                        >  1:662269/1‑239 (MQ=255)
                                                        gggagtgtttCTGTAATGCCGTCAACCACATTATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCa                                                                                                                                                                                       <  2:810728/231‑1 (MQ=255)
                                                         ggATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTGCCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCa                                                                                                                                                                                       <  2:562809/238‑1 (MQ=255)
                                                                   tGTAATGCCGGCAACAACATGATGAATAACCCTCTATTACCTGCCCCCGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCAt                                                                                                                                                                                      <  2:549467/229‑1 (MQ=255)
                                                                     tAATGCCGTCAAAAAAAAGATGAATACCCCTCTATTACCTGCCCACGGTAAAGGCGTTCGGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGc                                                                                                                                                                          <  1:234268/238‑1 (MQ=255)
                                                                          ccGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTACAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGGG‑GCGTTTCTTCCCTTCTGTCGGGTTTACCGACTGGGTAC‑CGGCGTGGTT‑AATATTTTCTCAtttctcactcggcgcaat                                                                                                                                                                     >  1:229524/1‑223 (MQ=255)
                                                                          ccGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATc                                                                                                                                                                                                                      <  1:621877/192‑1 (MQ=255)
                                                                          ccGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATc                                                                                                                                                                                                                      >  2:621877/1‑192 (MQ=255)
                                                                             tCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATAACCATTGGC‑GGTAATAGGCTGGGG‑CCGGTTTTTTCATTTTTTCGGGTTAACGCACTGGGAACGCCCGTGGAAC‑CATATTTTCACAttccatcaggccgttccgct                                                                                                                                                                   >  1:549467/1‑221 (MQ=255)
                                                                               aacaacATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGGG‑CTGTTTTTCTCATTCTGTCGGGTTAAAGAACTGGATAC‑CCCGTTGAAT‑CATATTTTCCCCTTGCCTACttggtccatttttt                                                                                                                                                                 >  2:234268/1‑225 (MQ=255)
                                                                                cccacattttatattaCCCTTTTTTCCCTGCCCCCGGTAAAGGCGTTTTGGTTGCCCTGCGTCCGTTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATCCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCa                                                                                                                                                              <  1:406743/224‑1 (MQ=255)
                                                                                  aacaTGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCaaa                                                                                                                                                            <  2:506195/239‑1 (MQ=255)
                                                                                         tgaATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTTATCGGCCGGTG‑CGGTTTTTCCCTCTCTTTCGGGATTACGAACTGGGTAC‑CCCCTTTTATCCATTTTTTCCC‑TTACCTACTAGGCG‑CCTTTTTTCaaaaccaatg                                                                                                                                                     >  1:29334/1‑232 (MQ=255)
                                                                                           taTAACCGTCTAGTTCCTGCCCACGGTTAGTGCGTTGAGGGTGCAGTACGGCCCGTCTGGGGTATCCGAGTTGAACAAGCGCTGTCTGTGTG‑TCGCCCTAACCGTACCGGCGATATCAGGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTGTCTCTCATTCTGTCGGATTAACGATTTGGATAC‑CGCGTGGAAT‑CACATTTGCCCATTGCCCACT‑GGCGACACTTTCTCAACCCCCATGGt                                                                                                                                                   <  1:300522/238‑1 (MQ=255)
                                                                                            attccccTCTATTCCCTCCCCGCCGTAAAGGTCTTCTGGTGGCCCTAAGTCCGTTACCGGGTATCCGCGTTAAACACGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCTATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGAGCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATGGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTc                                                                                                                                                  <  1:859893/235‑1 (MQ=255)
                                                                                             tAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAAGCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCATACT‑GGCGACATTTTCTCAAACCGTATGGTc                                                                                                                                                  <  2:191875/238‑1 (MQ=255)
                                                                                                          ccTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTTTGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTCAT‑GGCGACATTTTTTCAAACCGTATGGTCTGGTTTGaaattgt                                                                                                                                    >  2:259730/1‑234 (MQ=255)
                                                                                                             gCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTTGGG‑CTGTTTCTTCCATTCTGTCGGAGTAACGATCTGGGTAA‑CCCGTTTGAT‑CATCTTTTTCCCTTGCCTACC‑TGCGACCATTTTTCAAACCCGAAGTGCTGGTTTTAAgatgaaaa                                                                                                                                 >  1:360512/1‑236 (MQ=255)
                                                                                                              cccACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTTTCGGGTTAACGATCTGGGAAC‑CCCGTTGTAT‑CATATTTTCCCATTGCCTACT‑TGCGGCATTTTCTCCAACCCCATGGGCTGGGTTTAAGATGTTCaa                                                                                                                                >  2:580907/1‑239 (MQ=255)
                                                                                                              cccACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGCTAC‑CGCGTTGAAT‑CATATTTTCCCATTTCCTACT‑GGGGACATTTTCTCAAACCGTATGGGTCGGTTTGAAAATGTTCaa                                                                                                                                >  2:868072/1‑239 (MQ=255)
                                                                                                                ctgtgTAAAGGCGTTTTTTTCCCCCTTAGGCCGGTGTCGGGTTTCCGCGATGGACATCCCTTGACGCTGGG‑TCGCCCTAAGCGTACCTGCGATATCATGGCCATTGGC‑TGTAATCCGCTGGTG‑CGGTTTCGCTCATTCTGTCGGATTAGCGATCTGGTTAT‑CTCGTTGAAA‑CATATTGTACCATTGCCTACA‑GGCGACATTATCTCAAACCTTATGGTCTGGTTTGAAGATGATCAAAt                                                                                                                              <  1:453509/235‑1 (MQ=255)
                                                                                                                cACGGTAAAGGCGTTCTGGTTCCCCTGCGGGCGGTACCGGGTTACCGCGGTGAACAAGCCCTGGCGGTGTG‑TCGCCCTAACCGTACCGGTGAGATCAGTACCATTGGC‑GGTAAGCGGCTGGGG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAAt                                                                                                                              <  1:192994/239‑1 (MQ=255)
                                                                                                                     tAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCaa                                                                                                                                                             <  2:842636/203‑1 (MQ=255)
                                                                                                                      aaaGGCTTTCTGGTTGCCCTTCTGCCGGTTCCGGGTTTCCGGGTTGAACAAGCCCCGACGCTGTG‑TCCCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTg                                                                                                                         <  1:358212/238‑1 (MQ=255)
                                                                                                                        aGGCTTTTTGGTTGCCCTTCGGCCGGTACCGGGTATCCCCGTTGAACAAGCCCTGACGCTTTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCg                                                                                                                      <  1:723671/239‑1 (MQ=255)
                                                                                                                             ttCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTGCCTGTTT‑TCTCATTTTTTCGGGTTTACGGTCTGGATAA‑CGCGTTTGAA‑CCTATTTTCCCCTTTCCTACC‑TGCGACATTTTTCCAAACCGTCTAGTCTGGGTTTAAgtagactaaaaaaaggcagcgagc                                                                                                                 >  1:159718/1‑216 (MQ=255)
                                                                                                                             ttCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGGTTAACGATCTGGGTAC‑CGCGTTTAAT‑CATATTTTCCCATTGCCTACT‑TGCGCCATTTTCTCCAACCGTATGGTCTGGTTTTACGAGTGTCAAATCAGTGCCGGG‑Ctt                                                                                                                >  2:346227/1‑238 (MQ=255)
                                                                                                                             ttCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAAc                                                                                                                                                                                                                          >  2:675523/1‑137 (MQ=255)
                                                                                                                             ttCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAAc                                                                                                                                                                                                                          <  1:675523/137‑1 (MQ=255)
                                                                                                                                              cccgtccccggtttcccttttGAACAAGCCCTGGCCCTGTG‑TTGCCCTAACCGTTCCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCgcttgct                                                                                               <  1:868072/220‑1 (MQ=255)
                                                                                                                                              ccGGGACCGGGTTCCCGGGTGGAACAAGCCCTGAGGC‑GTGTTCGCCCTAACCGTGCCGTCGATAACATGCCCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCCGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCgcttgct                                                                                               <  1:624920/239‑1 (MQ=255)
                                                                                                                                              ccGGGACCCGGTATCCGCGTTGAACAATCCCTGACGTCGTG‑TCGCACTAACCGTACCGGAGATGACATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATCC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCGCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCgcttgct                                                                                               <  1:986372/239‑1 (MQ=255)
                                                                                                                                               cGGTACCGGGGAGCCGCGTTCAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCgcttgct                                                                                               <  2:229524/238‑1 (MQ=255)
                                                                                                                                                  tACCGGGTATCCGCGTTGAACAAGCCCGGACGCTGGT‑GCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTGCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGc                                                                                                                        <  1:200971/211‑1 (MQ=255)
                                                                                                                                                   tccGGGTATCCGCGTTGAACAAGCCCTGACGCTTTG‑TCGCCCTAACCGTACCGGCGATATCAGGACCATTGGC‑GGTAATCGGCTGGGG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTGGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑ccc                                                                                        <  2:181521/238‑1 (MQ=255)
                                                                                                                                                   tccGGGGATCCGCTTTGAACAAGCCCTTGCGCTGTT‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑ccc                                                                                        <  1:681549/238‑1 (MQ=255)
                                                                                                                                                     ggagatgtttataagggacaGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCa                                                                                      <  2:356001/220‑1 (MQ=255)
                                                                                                                                                      gggTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGt                                                                                                                                                        <  2:238123/175‑1 (MQ=255)
                                                                                                                                                                  tGAACAAGCCCTTGCCGTGTG‑GCGCCCTAACCGTACCGGCGATATCATGACCAGTGGC‑GGTAATCGGCTGGTG‑CTGGTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGc                                                                         <  2:29334/238‑1 (MQ=255)
                                                                                                                                                                              gACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTTTCGGATTAAc                                                                                                                                                                                                                          >  2:261048/1‑88 (MQ=255)
                                                                                                                                                                              gACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAAc                                                                                                                                                                                                                          <  1:261048/88‑1 (MQ=255)
                                                                                                                                                                              gACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCtt                                                            <  2:359822/239‑1 (MQ=255)
                                                                                                                                                                                     tg‑tCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTTTCTCCTTTTTTCGGGTTAACGATCTGGATAC‑CCCGTTGGAT‑CACATTTTCCCATTGCCTACT‑TGCGCGATTTTTTTCAACCGTATGGCCTGGTTTTGAAATAGTCAAATAAAGACCGTG‑CTCGTTCAGAAGCGCTTGCTCTG‑CCCCGAAAA‑AAGggggggtcttttttacacctctt                                                       >  1:356001/1‑219 (MQ=255)
                                                                                                                                                                                     gg‑tCGCCCTACTCGTGCGGGTGATATCATGACCTTTGGC‑GGTAATCGGCTGGTG‑CTGGTTCTCTCAGTCGGGCGGATTAACGATGTGGATAC‑CGCGTGGAAT‑CATATTTGCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCa                                                     <  1:430898/238‑1 (MQ=255)
                                                                                                                                                                                      ttgcGCCCTAACCGTACCGGCGATTTCATGGCCATTGGC‑GGTAATCGGTGTGTT‑CTGTTTCTCTCATTCTGACGTATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCa                                                     <  2:464278/236‑1 (MQ=255)
                                                                                                                                                                                      gttcccccTTACCGTCCCGGGTATATCATGTCCCTTGGC‑GGTAAGCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCa                                                     <  1:861656/234‑1 (MQ=255)
                                                                                                                                                                                      gttCGCCCTAACCCTACCGGCGATATCATGGCCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTGCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCa                                                     <  2:108374/237‑1 (MQ=255)
                                                                                                                                                                                      ggtcccttttgctgtgccGGGGATTTCATTGCCATTGGT‑GCTAGTGGGTGGGTG‑CTGTTGCTCTCATTCTGGCGGATTAACGACCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTGTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCa                                                     <  1:62046/223‑1 (MQ=255)
                                                                                                                                                                                      g‑tCGCCCTAACCGTACCGGCGATATCATGGCCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCa                                                     <  2:735343/238‑1 (MQ=255)
                                                                                                                                                                                      caccccccTACCCGTTCCGGCGATTTCATGCCCATTGGC‑GGTAATCGGCTGGTG‑CTTTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCa                                                     <  1:180878/234‑1 (MQ=255)
                                                                                                                                                                                        tCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCCGGGG‑CTGTTTCTCTCATTTTGTCGGATTAAAGAGCTGGAGAA‑CGGGGTGGAT‑TATTTTTTCCCATTGCCTCAT‑GGCGACAATTTTTTAAAACGTTTTGTCCGGTTTTGAGTTGGACAAAGCAGTGCACAGCCGGGGGCAGAT‑CGACTGCT‑TG‑CCCCCGAACAAATggggaagcccggggcctaaacaccaa                                                    >  1:977298/1‑217 (MQ=255)
                                                                                                                                                                                            ccgaaCCGGTCCGGCGATTTCAGGTACATTGGG‑GGTTAGCGGCTGGGT‑CTTTTTCTCTCATTCTTGCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATAGTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGtt                                                <  2:252688/235‑1 (MQ=255)
                                                                                                                                                                                                     aCCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCTTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGt                                      <  2:664873/239‑1 (MQ=255)
                                                                                                                                                                                                      ccGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGGG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTTCCTACT‑GGCGACATTTTCTCAAAACGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAA‑AATGGGGCCTGCCGCTGCCTTTAACGCCAAGTTCTaaa                                           >  2:180878/1‑233 (MQ=255)
                                                                                                                                                                                                         gCGATATCATGACCATTTGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTTGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtcatc                                  <  2:856265/239‑1 (MQ=255)
                                                                                                                                                                                                                cATGACCATTGGC‑GGTAATCGGCTGGGG‑CTGTTTCTCTCATTCTTTCGGATTAACGATCTGGATAC‑CGCGTTGTAT‑CCTATTTTTCCATTGCCTACT‑GGCGGCATTTTTTCAAACCGTATTGTCTGGTTTTAAGATGATACAATCAGTGCCCCG‑CTGGTGCAGATGCCCTTTCT‑TGCCCCCAGAA‑‑AATGGGGGATGCCCCCGCCCTTAACGCCAAGATCTAAACCCGCCACCAAAGCCg                           >  1:55960/1‑239 (MQ=255)
                                                                                                                                                                                                                    aCCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑cgcg                                                                                                                                                                                                          <  1:752065/66‑1 (MQ=255)
                                                                                                                                                                                                                    aCCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑cgcg                                                                                                                                                                                                          >  2:752065/1‑66 (MQ=255)
                                                                                                                                                                                                                     ttaggtccgggTTATGGGGTGGGG‑TTGTTTCTCTTTTTTGGTCGGATTAACGATCTTGATCC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCa                       <  1:389492/230‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATTTTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTgg       <  2:662269/239‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCTGGTG‑CTGTTGCTCTCATTCTGTCGGATTAACGATCTGGAGAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCGACT‑GGCGACAGTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTgg       <  2:489979/239‑1 (MQ=255)
                                                                                                                                                                                                                                          gTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCgcttgct                                                                                               >  2:188609/1‑149 (MQ=255)
                                                                                                                                                                                                                                          gTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCgcttgct                                                                                               <  1:188609/149‑1 (MQ=255)
                                                                                                                                                                                                                                          gTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCa                                                                        >  2:984048/1‑169 (MQ=255)
                                                                                                                                                                                                                                          gTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCa                                                                        <  1:984048/169‑1 (MQ=255)
                                                                                                                                                                                                                                          gTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGccctgcctcttt  >  2:389492/1‑229 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
TATCACTACCTTGCTGTCAATGACCCAGCCGCTCAAACTGCGTGGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTG‑TCGCCCTAACCGTACCGGCGATATCATGACCATTGGC‑GGTAATCGGCTGGTG‑CTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAC‑CGCGTTGAAT‑CATATTTTCCCATTGCCTACT‑GGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAG‑CTGGTGCAGATGCGCTTGCT‑TG‑CCCCAGAAC‑AATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACG  >  NC_000913/3697506‑3697976

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