Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 3,845,705 +G intergenic (‑72/+26) PP_3394 ← / ← PP_3395 putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/Transcriptional regulator, LysR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154513,845,7031.G92.6% 74.1 / ‑3.0 30intergenic (‑70/+28)PP_3394/PP_3395putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/Transcriptional regulator, LysR family
Reads supporting (aligned to +/- strand):  ref base . (2/0);  major base G (12/13);  minor base A (2/0);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AGCACGGTTACTTGCAGGTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑C‑GG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑GGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCA  >  AE015451/3845404‑3845990
                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                   
aGCACGGTTACTTGCAGGTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTcg‑c‑g                                                                                                                                                                                                                                                                                                    >  1:132303/1‑301 (MQ=255)
                     gggtagggtgggggtggggCCACCACTGCGCGGGCAGCGGCCATCTGCGGCAGCAAGCGTGCCGGAACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCAc                                                                                                                                                                                                                                                                               <  2:620187/283‑1 (MQ=255)
                                     gAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCcgcg                                                                                                                                                                                                                                                               <  1:78691/301‑1 (MQ=255)
                                         caccacTTCGCGGGCATCGGCCATCTGCGCCAGCAAGCGGTCCGGCACGAACTACGTGCCTTCCATGTGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCt                                                                                                                                                                                                                                                           <  2:461046/301‑1 (MQ=255)
                                              ctgcgcgGGCATCGGCCATCTGCGGCAGCAAGCGGGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAg                                                                                                                                                                                                                                                      <  1:659004/298‑1 (MQ=255)
                                              cTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGAATTCAGTGACTGCTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGAGTTGCGCCAGACTAATAGCTGGATGCCGGGCCGAAAAGACGTTTTTCG‑CAGT‑CTGAGCCATCGTAAGACATGATGGATCAGccagagtccttccgg                                                                                                                                                                                                                                                      >  2:386505/1‑288 (MQ=255)
                                                   cgGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTCTTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CTGG‑CTGATCCATCGCCAGCCACGATGGacacccccgacccccgcccgtcccg                                                                                                                                                                                                                                                 >  2:298482/1‑277 (MQ=255)
                                                   cgGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCC‑GG‑CTGATCCATCGTCAGCCACGATGGAGCACCCCGCGCCCTGCCAGGCCgg                                                                                                                                                                                                                                                 >  1:12568/1‑301 (MQ=255)
                                                       caccGGCCATCTGCGGAAGCAAGCGGGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCGCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGCCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCa                                                                                                                                                                                                                                             <  2:109599/298‑1 (MQ=255)
                                                            gCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgcag                                                                                                                                                                                                                                        >  1:52207/1‑301 (MQ=255)
                                                                    cGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCCGAAAAGACGTGTTTCG‑CAGG‑CTGATCCATCGACAGCCACGATGGGTCAGCCCGCGCCCTGCCAGGCCGGTCCAccgagcaaacggt                                                                                                                                                                                                                                >  1:114435/1‑289 (MQ=255)
                                                                        agcaAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtggag                                                                                                                                                                                                                            >  1:611292/1‑299 (MQ=255)
                                                                             gcttGCCGGCCCGAACGACGCGACTTCCATATGGCGCACGCCGGGGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCa                                                                                                                                                                                                                       <  1:139785/298‑1 (MQ=255)
                                                                             gCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCa                                                                                                                                                                                                                       <  1:86392/301‑1 (MQ=255)
                                                                                                                 cACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgc                                                                                                                                                                                   >  1:207268/1‑301 (MQ=255)
                                                                                                                           ccGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCCACAGGTGGGCAAGCAGGCGCCGAGCCGACAGCCAGAGCGTGTGCCCGACGGGCATGCAGaa                                                                                                                                                                         >  1:424292/1‑300 (MQ=255)
                                                                                                                                                         gtgtctgTCGGCATGATCACTTGCAGCCTCTGCAAGCCATCGCGCAGGCCTACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCgtcagt                                                                                                                                          <  1:263036/297‑1 (MQ=255)
                                                                                                                                                                     aTGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGGCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGAc                                                                                                                              >  1:18831/1‑301 (MQ=255)
                                                                                                                                                                                      ttGCAGGCCATCGCGCAGGCCGACTTCATTGATTGTTATGGTTGTCATCGCGGTTTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAAcg                                                                                                            <  2:207184/300‑1 (MQ=255)
                                                                                                                                                                                              cttagcgcaggccagtttCATTGTGTGTTATTGTTTGCATCGCGGTCTCTCTGTCTGAAGCGTTGCGACAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAg                                                                                                    <  2:132246/285‑1 (MQ=255)
                                                                                                                                                                                                                 attgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGGCAAGCAGGCGCCGAGCCGACAGCGAGAGCGGGTGCCGGTAGCGCATGCAGATGACGAATTCACCTAG‑AGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCACGTCCAACAACGCCGCAtcacgggagccggctgcttgcc                                                                                 >  2:159868/1‑280 (MQ=255)
                                                                                                                                                                                                                              gttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGCTAGCCACGATGGATCAGCCCCCGCCCTGCCAGGCCGGTGCACGCCGCAACAGGCGGTCAAGCAGGCGCCGCGCCGACAGCCAGAGCGTGTCCCGGTCGCGCCTGCAGACGACGAATTCACCAAG‑CGCCCAGCCGTCAGTCAGCGCAATGACC‑CGCAGTTCGAACAACGCCCCATAGCC‑GGCCACGCCCCATTGCGCGAACATCGCCag                                                                   >  1:332625/1‑301 (MQ=255)
                                                                                                                                                                                                                                        cgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCTTCATCCACGATGGATCAGCCCGCGCCCTGCCCGTCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCACCGAGCCGACCGCGAAACCGTGTGCCGGTCGCGCATCCAGATGACGAATTCACCAAG‑CCCCCAGTCGTCAGTCGGCGCAAT‑ACCACGCAGGCCCAACAACGCCCCATACCG‑GGACACGCCCTCTTCCGGGAACGCCGCCaggcccagcccc                                                         >  2:95456/1‑294 (MQ=255)
                                                                                                                                                                                                                                               tgtgtgtgAAGAGTTGCGCCAGACTAAGGGCTGTGTGAAGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATTGTCAGCCACGATGGAGCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑GGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCa                                                  <  2:167408/295‑1 (MQ=255)
                                                                                                                                                                                                                                               tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGGAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑GGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCa                                                  <  1:77343/301‑1 (MQ=255)
                                                                                                                                                                                                                                                        gccgttcggccagcatagagccgggaagccGGGCCGAAGAGACGTGTTTCG‑CGGG‑CTGTTCCAGCGTCGTCCACGGTGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑GGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGt                                         <  2:18830/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTTATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑GGACACGGCCTCTGTCGGGAACACCGCCAGGCCAAGGCCAAACTGCACAAGACGGTATGCcgcc                                 >  2:567663/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                   aGACTACGCGCTGGATGCCGTGCAGAACAAACGTGTTTCC‑C‑GGCCTGATCCCTCGTCAGCCCCGATGGCTCAGCCCACGCCCTGCCAGGCCGGTGCCCGCAGCAACAGGTCGGGAAGCAGGCACCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCGTGCAGATGCCGGCTTTA‑CAAGCCGCCCACTCGTCAGTCAGCGCGATCACC‑CGCCGCGCGCACAACGCCGTAT‑CCGCGGACACGGCCTCTTGCGTGCACACGGCCACGCCAACGCCAAACTGCACCAGGCgggttgctccgtcc                              >  2:205351/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                      cTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑GGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAAATGCACCAGACGGTATGCCGCGTCGaaa                           >  1:664205/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                             aaaGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑GGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCACACGTTATGCCGCGTCGAAAGACGACACGTACGAACGAAAACg    >  2:1400/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               aGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGCCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑CGACACGGCCTCTTGCGGCAACACCGCCAGGCCAAGGCCAAACCGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACTACAACGCa  >  2:124327/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                   
AGCACGGTTACTTGCAGGTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑C‑GG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAG‑CGCCCAGTCGTCAGTCAGCGCAATGACC‑CGCAGGTCGAACAACGCCGCATAGCG‑GGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCA  >  AE015451/3845404‑3845990

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑