Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 1,499,498 1 bp→CG intergenic (+75/‑89) trpS → / → zapE Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,499,4971.C96.8% 108.3 / ‑5.0 31intergenic (+74/‑90)trpS/zapETryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (17/13);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE0154511,499,4980TG100.0% 108.9 / NA 31intergenic (+75/‑89)trpS/zapETryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/13);  total (18/13)

CGCGCAGCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCC‑GGTGTT‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCTGGGGTGGGGTGGGGCGAGG  >  AE015451/1499215‑1499798
                                                                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                                                                            
cacgccgctcttccgatctcAAGGGTGCTGGCGCCGTGAGCCAGCAGCTGCAGCAAGGTGACGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGATCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTcc                                                                                                                                                                                                                                                                                               <  1:436634/282‑1 (MQ=255)
              gACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGc                                                                                                                                                                                                                                                                                 <  4:204674/301‑1 (MQ=255)
                   cAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTccagatcggaagagcacacgg                                                                                                                                                                                                                                                                            >  3:436636/1‑282 (MQ=255)
                    aaGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTaa                                                                                                                                                                                                                                                                           >  3:111992/1‑301 (MQ=255)
                               ggagttcagacgtgtgctcttccgatctGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                                                                                                                                                                                                <  4:121215/274‑1 (MQ=255)
                                     gAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                                                                                                                                                                          >  3:226597/1‑301 (MQ=255)
                                                           ggCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttgg                                                                                                                                                                                                                                    >  3:62135/1‑300 (MQ=255)
                                                           ggCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACCGGCTTTGCCGGGTGTC‑GCGGGCACGCCCGCTTCCACCGACTGATTGCCATTACGCGGGCCTaacagatcagaagagagtcgtagaagggaag                                                                                                                                                                                                                                    >  2:121190/1‑271 (MQ=255)
                                                              gAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                                                                                                                                                                                                                 <  1:62133/301‑1 (MQ=255)
                                                                       gTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACCGTGACGGCCCGGTTTTTTTTTCCGCGCTAACGaa                                                                                                                                                                                                                        >  2:459341/1‑301 (MQ=255)
                                                                                  aGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                                                                                                                                                                             <  4:459470/301‑1 (MQ=255)
                                                                                                                       cccggtccagggCCACAGTGCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGACGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct                                                                                                                                                                        <  2:531411/292‑1 (MQ=255)
                                                                                                                           ccgatcggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGTCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt                                                                                                                                                                    <  1:415192/295‑1 (MQ=255)
                                                                                                                               cGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGAAGATCTGAAACGTCCGGActagatcgga                                                                                                                                                                >  4:531561/1‑293 (MQ=255)
                                                                                                                                  gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGTCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtcttagatcgg                                                                                                                                                             >  3:415194/1‑294 (MQ=255)
                                                                                                                                                           cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACGGCCCGTTTTTTTTTGCCCCGCTAACGAATCATGACTCCCTTTAAACGCTATCAATCAGATCTGAAACGTCCCGACTTCTTCCATTACGCGGCGCAAGAAACTGACGTGCGTc                                                                                                                                    >  1:114820/1‑301 (MQ=255)
                                                                                                                                                               aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAAGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAACCGTCCCGACTTCTTCCATGAAGCGGCACAGGAACCTGCGGTGCGTccccc                                                                                                                                >  1:173481/1‑297 (MQ=255)
                                                                                                                                                                ggCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTc                                                                                                                               <  1:226599/300‑1 (MQ=255)
                                                                                                                                                                  gggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCa                                                                                                                             <  1:111989/300‑1 (MQ=255)
                                                                                                                                                                                       aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGt                                                                                                        <  3:114823/301‑1 (MQ=255)
                                                                                                                                                                                                          ggACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAAGCGGGGCCGGAAACTGAGGTGCGTCACCTACCGCGCCTGTCCGCAGACCTGGGGCACGCGCAGaacaaaaca                                                                                     >  1:365866/1‑297 (MQ=255)
                                                                                                                                                                                                                   cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAACCTGCGGTGCGTCCCCCCCAGCGCCTGTAAGACGACCTGGCGCACGCGCCGCACAACAAGCCGGGccag                                                                            >  2:694019/1‑298 (MQ=255)
                                                                                                                                                                                                                        aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGGGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGGTCgg                                                                       >  4:57049/1‑301 (MQ=255)
                                                                                                                                                                                                                                 gAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGtt                                                              <  3:173483/301‑1 (MQ=255)
                                                                                                                                                                                                                                    cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCAGGTTCGGCAAGCTGTTCgg                                                           >  4:448945/1‑301 (MQ=255)
                                                                                                                                                                                                                                                  aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCACGAAGGACCAGa                                             >  4:166176/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         acggtgtgctcttccgatctgCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGt                                      <  3:597384/281‑1 (MQ=255)
                                                                                                                                                                                                                                                           ccgcGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGCTTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGCTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGAACTGGTGCAAGCGCAGAACAACAAGCCTGGCAAGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGCCTCCGGtaa                                    >  1:163419/1‑299 (MQ=255)
                                                                                                                                                                                                                                                           ccgcGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCa                                    <  3:365867/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGtagatcggaagagcgtcgtgg                  >  1:597382/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                ttGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGGTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACCTCggggggg               >  2:781754/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                   gatcttgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGGTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCtggggtggg             <  2:167486/290‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ccggtgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCtggggtgggg            >  3:613424/11‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        tgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCtggggtgggagacg        >  4:167528/7‑296 (MQ=255)
                                                                                                                                                                                                                                                                                         gttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTACCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCtgggggggggtgggg       >  2:78664/6‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          ttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCtggggtggggtggggg      >  3:333741/5‑300 (MQ=255)
                                                                                                                                                                                                                                                                                              cGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAACCTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCTGGGGTGGGGTGGGGCGAgg  <  4:782034/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                                                                            
CGCGCAGCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCC‑GGTGTT‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCTGGGGTGGGGTGGGGCGAGG  >  AE015451/1499215‑1499798

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