Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,499,498 1 bp→CG 100% 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.C87.8% 145.6 / 16.5 49intergenic (+74/‑90)trpS/zapETryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base C (11/32);  total (11/38)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*AE0154511,499,4980TG89.8% 146.4 / 9.6 49intergenic (+75/‑89)trpS/zapETryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (11/33);  minor base C (0/5);  total (11/38)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑TT‑T‑CGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCAC  >  AE015451/1499276‑1499700
                                                                                                                                                                                                                                |||                                                                                                                                                                                                          
cGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑CGCGGGCAcg‑c                                                                                                                                                                                                   <  2:970582/233‑1 (MQ=255)
                       gCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                                                                                           <  1:578441/235‑1 (MQ=38)
                              cTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                                                                                                   <  1:1001451/236‑1 (MQ=255)
                              cTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑CGCGGGCAcg‑cccgctc                                                                                                                                                                                             <  2:736231/209‑6 (MQ=255)
                               ttCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                                                                                     <  2:707038/233‑1 (MQ=255)
                               tgcGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑CGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                                                                                     <  1:549797/231‑1 (MQ=255)
                                       gtgGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                                                                                  <  1:391689/228‑1 (MQ=38)
                                                               cGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                                                                        >  2:944593/1‑214 (MQ=255)
                                                                   gggCCGCCGTTCGACAACGCGGAAGGCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCCTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCAcg‑c                                                                                                                                                                                                   <  2:652694/166‑1 (MQ=37)
                                                                    ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                                                                                             <  1:747164/236‑1 (MQ=255)
                                                                     gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                                                                                                                   <  1:669774/213‑1 (MQ=38)
                                                                              cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                                                                                   <  1:358276/236‑1 (MQ=38)
                                                                                          aGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                                                                                        <  1:719448/235‑1 (MQ=38)
                                                                                            cgcGTGATGCGGCCATGCAGAGCAGGAACGAAGCCGTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGCTC‑G‑CGGGCACGCCCGCTACCACAGACt                                                                                                                                                                                     <  2:185847/154‑1 (MQ=21)
                                                                                            cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                                                                                                                            >  2:473274/1‑213 (MQ=255)
                                                                                               gggATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCcaccgttt                                                                                                                                                                                     <  1:35866/151‑6 (MQ=21)
                                                                                                 gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                                                                                                                            >  1:144530/1‑208 (MQ=38)
                                                                                                 gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                                                                          >  2:707030/1‑226 (MQ=38)
                                                                                                  aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt                                                                                               >  1:1041161/1‑236 (MQ=255)
                                                                                                    gCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                                                                              <  1:707033/235‑1 (MQ=39)
                                                                                                      ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                                                                                    <  2:854112/195‑1 (MQ=37)
                                                                                                                agCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCTCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                                                                                                                 <  2:796116/156‑1 (MQ=25)
                                                                                                                  ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                                                                                  <  1:384447/201‑1 (MQ=38)
                                                                                                                    tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                                                                                   <  1:438012/198‑1 (MQ=37)
                                                                                                                           aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                                                                                                               <  1:631251/179‑1 (MQ=37)
                                                                                                                                     ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAATGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATg                                                             <  1:1013613/235‑1 (MQ=38)
                                                                                                                                           cagGACTGAGGCTGACAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                                                                                   <  1:521699/175‑1 (MQ=25)
                                                                                                                                                      cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGCTTGCCATTCAGCGCGCCTATCGCCACAGt                                                                                                                                                      >  2:22538/1‑129 (MQ=17)
                                                                                                                                                      cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTGTTTTGGTGCCTCGCTAACGAAGCATGAATcccc                                                                                                          >  1:115181/1‑172 (MQ=25)
                                                                                                                                                        gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCCCGGCTTTTCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACg                                                                                                                                                  >  2:634753/1‑131 (MQ=21)
                                                                                                                                                                   ggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCAcc                               <  1:697088/235‑1 (MQ=255)
                                                                                                                                                                   ggAGCGGGCGGGCCCGCGAAACCGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCt                                                                                    <  1:628628/182‑1 (MQ=37)
                                                                                                                                                                          gCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                                                                                                                      >  2:695265/1‑109 (MQ=18)
                                                                                                                                                                            gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                                                                              >  2:1053043/1‑163 (MQ=38)
                                                                                                                                                                            gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCACGACTCCCTTAGAACGCTATc                                                                                              <  1:1053046/163‑1 (MQ=37)
                                                                                                                                                                             tGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCTGGTG‑TT‑CGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCt                     <  1:825586/235‑1 (MQ=39)
                                                                                                                                                                               cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCTTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGt                   <  1:807779/235‑1 (MQ=39)
                                                                                                                                                                               cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑CGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGt                   <  1:999327/235‑1 (MQ=255)
                                                                                                                                                                                ccGCGAATCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAAc                                                                               <  2:913069/174‑1 (MQ=38)
                                                                                                                                                                                ccGCGAAACAGGCGCCGCGGTGCATTGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAg                          <  1:450233/227‑1 (MQ=39)
                                                                                                                                                                                 cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATg                                                             <  1:614814/191‑1 (MQ=38)
                                                                                                                                                                                         aGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcgc                                                         <  1:950264/187‑1 (MQ=39)
                                                                                                                                                                                           gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAAc                                                                               <  2:543211/163‑1 (MQ=38)
                                                                                                                                                                                               cgcgGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCAc  >  1:1034513/1‑236 (MQ=255)
                                                                                                                                                                                                 cgGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCa   <  2:1041158/233‑1 (MQ=255)
                                                                                                                                                                                                          gcacaggCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTg                                      <  2:734835/184‑1 (MQ=255)
                                                                                                                                                                                                            accGGCTTTGCCGGTG‑TTCG‑CGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTcc                            <  2:683715/195‑1 (MQ=255)
                                                                                                                                                                                                             ccGGCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                                                                                      <  2:268635/121‑1 (MQ=38)
                                                                                                                                                                                                             caggCTTTGCCGGTG‑TTCG‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTg                                      <  1:101259/184‑1 (MQ=255)
                                                                                                                                                                                                                                 gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCatg                               <  2:767567/174‑3 (MQ=255)
                                                                                                                                                                                                                                |||                                                                                                                                                                                                          
CGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑TT‑T‑CGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCAC  >  AE015451/1499276‑1499700

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