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.C97.7% 145.3 / ‑5.7 43intergenic (+74/‑90)trpS/zapETryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (21/21);  minor base A (1/0);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
*AE0154511,499,4980TG100.0% 152.5 / NA 43intergenic (+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 (22/21);  total (22/21)

GCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGGCAC‑G‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCTGGGGTGGGG  >  AE015451/1499221‑1499788
                                                                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                                                                    
gCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAc‑g‑ccc                                                                                                                                                                                                                                                                                          <  2:70467/290‑2 (MQ=255)
           ggCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCa                                                                                                                                                                                                                                                                     <  1:63720/300‑1 (MQ=255)
                  gTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAg                                                                                                                                                                                                                                                                                <  2:143230/282‑1 (MQ=255)
                                   aagtagcTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTg                                                                                                                                                                                                                                             <  1:520606/296‑1 (MQ=255)
                                   aagcaGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                                                                                                                                                                                                              <  1:620229/299‑1 (MQ=255)
                                         cTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGCGTTCGCGGGCCC‑GCCC‑GCCCCCACACACCGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                                                                                                                                                                                      >  1:557827/1‑301 (MQ=255)
                                           gcagaaAGGTGGCAAGGCCGATGTCCGAATCGAAGGCCTGAGCTTCAGCCTGTGGCTCAATGGCGAGCAGGTTACCAACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCAGCCATTCAGAGCCTGAACGAAGCCTTGGAACCACAGCAGGACTGAGGCTGCCAAACGGTGGGGGCGGGCGTACCCGCGAAACAGGCACCACGGTGCATAGCACCGGCTTTGCCAGTGTTCGCGGGCACAC‑CC‑GCTCCCACAGACTGATTACCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                                                                                                                                                                                      <  1:417507/299‑1 (MQ=255)
                                              gcaAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTa                                                                                                                                                                                                                                                        <  2:197025/278‑1 (MQ=255)
                                              gcaAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCCCCGCCGTGCATGGCACCAGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCCCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                                                                                                                                                                                                         >  1:197029/1‑277 (MQ=255)
                                                          ggCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCg                                                                                                                                                                                                                      >  1:754107/1‑300 (MQ=255)
                                                             cGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCt                                                                                                                                                                                                                                                   >  1:192452/1‑268 (MQ=255)
                                                                     gCATCGAAGGCCTGGGCTTCGGCATGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCg                                                                                                                                                                                                                      <  2:192448/289‑1 (MQ=255)
                                                                       aTCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCcaca                                                                                                                                                                                                                                                                                 <  2:299416/228‑1 (MQ=39)
                                                                         cGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCCCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACCGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACGGCCCGTTTTTTGTTGCCTCGc                                                                                                                                                                                                                     >  2:393062/1‑286 (MQ=255)
                                                                          gAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCcac                                                                                                                                                                                                                                                                                  >  1:829946/1‑224 (MQ=38)
                                                                                cTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTACCACAGACTGCTTGcc                                                                                                                                                                                                                                                                      >  1:863053/1‑230 (MQ=38)
                                                                                        cGGCCGGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                                                                                                                                                                                        <  1:760523/300‑1 (MQ=255)
                                                                                                                             ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcg                                                                                                                                                  <  1:830733/301‑1 (MQ=255)
                                                                                                                                  gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTTCATGGCACCGGCTTTGCTGGTTTTCGCGGGCAC‑GCCC‑GCCCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                                                                                                                                                                >  2:779993/1‑202 (MQ=255)
                                                                                                                                  gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                                                                                                                                                                <  1:780019/202‑1 (MQ=38)
                                                                                                                                        cAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCTACTTCTTCCATGAcgc                                                                                                                                                 >  1:530289/1‑291 (MQ=255)
                                                                                                                                            gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCAC‑GCCC‑GCTCCCACAGACAGATTGCCATTAAGCGGGCCGATCg                                                                                                                                                                                                                                                     >  1:382129/1‑187 (MQ=255)
                                                                                                                                                             ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCCCCCACAGACTGACGGCCATTAAGCGGGCCTAACGCCACAGTGAc                                                                                                                                                                                                                                           >  2:837715/1‑180 (MQ=255)
                                                                                                                                                                 aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                                                                                                                                                                                  >  2:770034/1‑233 (MQ=39)
                                                                                                                                                                 aTTCAGAGCCTGAACGAAGCCTTGGAACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATg                                                                                                                                                     <  1:770060/262‑1 (MQ=39)
                                                                                                                                                                         ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCt                                                                                                                      <  2:754081/285‑1 (MQ=255)
                                                                                                                                                                                gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGGGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAg                                                                                                                                          <  2:817530/258‑1 (MQ=39)
                                                                                                                                                                                     cTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGGGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                                                                                                                                                                                                    >  2:617890/1‑195 (MQ=38)
                                                                                                                                                                                        ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATTACTCCCTTAGAACGCt                                                                                                                                                                                         >  2:872950/1‑203 (MQ=38)
                                                                                                                                                                                               cagcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGCGTTCGCGGGCAC‑G‑CCTGCCCCCACAGACTGATTGCCATTAAGCGGGCCTATAGCTACAGTGCAGGCCCGTTTTTTGTTGCCTCGCTAACGACTCATGACTCCCTTAGAACg                                                                                                                                                                                           >  1:408768/1‑194 (MQ=255)
                                                                                                                                                                                                   agGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa                                                                                                                                                       <  2:758558/226‑1 (MQ=39)
                                                                                                                                                                                                   agGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCCAACGAATCATGACTCCCTTAGAACGCTATCCAGCAGATCTGAAAc                                                                                                                                                                       >  1:758584/1‑210 (MQ=38)
                                                                                                                                                                                                   agGACTGAGGCCGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCACACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc                                                                                                                                                                                                                 <  2:257485/168‑1 (MQ=37)
                                                                                                                                                                                                                         gggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCCTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGATTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAATGTCCCGACTTCTTCCATg                                                                                                                                                     >  1:943616/1‑206 (MQ=38)
                                                                                                                                                                                                                                     gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGt                                                                                                                               <  2:382117/216‑1 (MQ=39)
                                                                                                                                                                                                                                                         ggTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTTCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTc                                >  1:907036/1‑291 (MQ=255)
                                                                                                                                                                                                                                                            gCATGGCACCGGCTTTGCCGGGGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTc                                <  2:907010/288‑1 (MQ=255)
                                                                                                                                                                                                                                                             catggcaccGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCaagaag                                          <  2:363195/269‑1 (MQ=255)
                                                                                                                                                                                                                                                             catggcaccGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCaagaag                                          >  1:363207/9‑277 (MQ=255)
                                                                                                                                                                                                                                                              atggcaccGGCTTCGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGCTTGCCAATAAGTGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAAGGTCCCGACTTCTTCCATg                                                                                                                                                     >  2:419341/8‑169 (MQ=37)
                                                                                                                                                                                                                                                                 gcaccGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggc                                                                                                                                              >  2:655028/5‑173 (MQ=39)
                                                                                                                                                                                                                                                                 gcaccGGCTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggc                                                                                                                                              <  1:655053/169‑1 (MQ=39)
                                                                                                                                                                                                                                                                        cTTTGCCGGTGTTCGCGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGCTTTTTGTTGCCTCGCTAACGAATCGTGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGTAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGt                                                                                                           >  2:126029/1‑201 (MQ=255)
                                                                                                                                                                                                                                                                               ggtgttcgcGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTTCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTTtggggggggg  >  1:339931/9‑299 (MQ=255)
                                                                                                                                                                                                                                                                               ggtgttcgcGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGTCCTGTACTTCtggggtggg   >  1:631587/9‑298 (MQ=255)
                                                                                                                                                                                                                                                                                       cGGGCAC‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCGTGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGTAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGAACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACt              <  1:126031/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                                                                    
GCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGGCAC‑G‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCAAGGGCCTGTACTTCTGGGGTGGGG  >  AE015451/1499221‑1499788

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