Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 215,383 A→G 6.9% E6963E (GAA→GAG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451215,3830AG 6.9% 261.5 / 8.5 100E6963E (GAA→GAGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (28/64);  new base G (2/5);  total (30/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

CAAGTTCGAACAGCTTACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCACC  >  AE015451/215152‑215614
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
cAAGTTCGAACAGCTTACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                                     <  2:873731/236‑1 (MQ=31)
  aGTTCGACCAGCTGACACTCGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGCGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                                     <  1:585160/234‑1 (MQ=9)
                  gCTGGACAAGCCGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                                <  1:93990/223‑1 (MQ=14)
                                 ggTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCCGCGCGCCATCCACCCa                                                                                                                                                                                                                     >  1:60450/1‑219 (MQ=17)
                                   tCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCa                                                                                                                                                                                                                             >  2:61109/1‑209 (MQ=31)
                                             ggTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGcc                                                                                                                                                                                                                              >  1:411908/1‑198 (MQ=31)
                                                    gacgaGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCt                                                                                                                                                                                   <  1:784876/234‑1 (MQ=37)
                                                         gCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                                <  2:710245/184‑1 (MQ=32)
                                                         gCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGcc                                                                                                                                                                               <  2:1077627/233‑1 (MQ=39)
                                                          ccAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGa                                                                                                                                                                            <  1:486151/235‑1 (MQ=37)
                                                           cAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATc                                                                                                                                                                          <  1:434791/236‑1 (MQ=38)
                                                                  gCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCt                                                                                                                                                                                                                 <  2:891911/190‑1 (MQ=34)
                                                                            gAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAc                                                                                                                                                            <  2:971172/233‑1 (MQ=37)
                                                                             aaGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTc                                                                                                                                                         <  1:92252/235‑1 (MQ=37)
                                                                              aGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                       <  1:1006956/236‑1 (MQ=39)
                                                                                   aCTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCGACGACCTGGTCGTGACCCTGAGCAACAACGCCCACGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCGCCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGAt                                                                                                                                                                           >  2:698539/1‑211 (MQ=255)
                                                                                    cTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtgt                                                                                                                                                   >  2:1078758/1‑234 (MQ=39)
                                                                                              gTCCCCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACgg                                                                                                                                       <  1:139766/236‑1 (MQ=37)
                                                                                                      cACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                       <  1:550256/212‑1 (MQ=38)
                                                                                                       aCCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                              <  1:926485/236‑1 (MQ=38)
                                                                                                        ccgccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                              <  2:779786/235‑1 (MQ=38)
                                                                                                          gccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTc                                                                                                                                                         <  2:1075079/206‑1 (MQ=37)
                                                                                                          gccgACCAGACCTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCg                                                                                                                                                                                                       <  2:427574/160‑1 (MQ=16)
                                                                                                            cgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                              <  1:897460/231‑1 (MQ=35)
                                                                                                            cgACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAc                                                                                                                                                            <  2:530200/201‑1 (MQ=17)
                                                                                                                 aGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAGGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGCCCGCtggg                                                                                                                                                  >  1:372372/1‑204 (MQ=25)
                                                                                                                      ccGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACt                                                                                                               <  1:646556/235‑1 (MQ=25)
                                                                                                                       cGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCt                                                                                                                                                                                                                 <  2:835590/137‑1 (MQ=34)
                                                                                                                        ggTGGCTGAGTACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAGGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCTGCGCGCCATATACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                              <  1:698540/235‑1 (MQ=16)
                                                                                                                        ggTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                              <  1:1078759/235‑1 (MQ=37)
                                                                                                                        ggTGGCCGAGAACGTCAAACCGACCTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCGGGCCAGATCAGCCTGGGCATCGCGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                              <  1:490021/235‑1 (MQ=255)
                                                                                                                                 gaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGGCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                              <  1:184612/226‑1 (MQ=21)
                                                                                                                                  aaCGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCgg                                                                                                    <  1:61109/235‑1 (MQ=33)
                                                                                                                                   aCGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGGGGATGTAGACGGTCGCACCTTCGAAAACCTGCaa                                                                                                                 >  2:292252/1‑221 (MQ=21)
                                                                                                                                   aCGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAACGCGATGACGTCTACAaggacg                                                                                                                                                                                     >  1:362349/1‑153 (MQ=25)
                                                                                                                                      tCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGcc                                                                                                                                                                               <  2:220843/156‑1 (MQ=35)
                                                                                                                                      gcAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAGCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACAccc                                                                                                                                                                                                                      <  1:817505/116‑1 (MQ=25)
                                                                                                                                        caGCCGACCGTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGGGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGAGGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                       <  2:23842/177‑1 (MQ=255)
                                                                                                                                         aGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCa                                                                                                                                                                                                                                     >  2:675794/1‑99 (MQ=18)
                                                                                                                                          gCCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtg                                                                                                                                                    <  2:793788/179‑1 (MQ=19)
                                                                                                                                           ccGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                        <  2:478963/142‑1 (MQ=34)
                                                                                                                                            cGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGCCCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                       <  2:984105/174‑1 (MQ=18)
                                                                                                                                            cGACCTACACCGGGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                              <  2:136222/199‑1 (MQ=18)
                                                                                                                                                 ttCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGGTGCGCAAGGCGATGACGTCTAGAACGACGCTGACCAGATAAGCCTGGGCATCa                                                                                                                                                             <  1:816646/163‑1 (MQ=255)
                                                                                                                                                 ttCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGtt                                                                                            >  2:228258/1‑228 (MQ=32)
                                                                                                                                                   cACCGTCCACGTCAACCAGCCGCTGGCGCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                                                                                                                                                             <  2:891602/129‑1 (MQ=25)
                                                                                                                                                        ggCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCa                                                                                                                                                                              <  2:319319/138‑1 (MQ=21)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACAcc                                                                                                                                                                                                                       <  1:692216/95‑1 (MQ=18)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACAcc                                                                                                                                                                                                                       >  2:692215/1‑95 (MQ=18)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                        <  2:320293/158‑1 (MQ=33)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCaa                                                                                         >  2:734600/1‑221 (MQ=33)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCaa                                                                                         <  1:734601/221‑1 (MQ=33)
                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATg                                                                                                                                              >  2:664358/1‑163 (MQ=33)
                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGACACTGTGCGGTGCGGCt                                                                                                  >  1:43877/1‑207 (MQ=17)
                                                                                                                                                                 aCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCAt                                                                                                                                                               <  2:1039283/145‑1 (MQ=34)
                                                                                                                                                                          tGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGc                                                            >  1:936355/1‑235 (MQ=33)
                                                                                                                                                                            gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGCCGGTCGCACGTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAagtagt                                                               >  2:715982/1‑230 (MQ=18)
                                                                                                                                                                            gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACCTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATac                                                                              <  2:1079073/215‑1 (MQ=17)
                                                                                                                                                                              ccACGAACTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACATCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGAc                                                                                    <  1:469776/207‑1 (MQ=17)
                                                                                                                                                                                  caCCGGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCGCGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                              <  2:260963/160‑1 (MQ=18)
                                                                                                                                                                                   acgtGGTCGTGGCCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATc                                                                                                                                                                          <  1:990400/113‑1 (MQ=12)
                                                                                                                                                                                     cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGCCCGCTGTGGATGTAGACGGTCGCACCTTCGACAACCt                                                                                                                     >  2:277761/1‑167 (MQ=18)
                                                                                                                                                                                      tGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACCCCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAagtagt                                                               <  1:57652/220‑1 (MQ=17)
                                                                                                                                                                                      gggTCGTGACCCTGAGCAACAACGCCCAGGTCACCAGCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACt                                                                                                               <  2:386988/171‑1 (MQ=18)
                                                                                                                                                                                        gTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACcg                                                 <  2:1089453/232‑1 (MQ=32)
                                                                                                                                                                                         tCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                        <  2:880428/96‑1 (MQ=18)
                                                                                                                                                                                           gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCg                                                                                              >  2:476531/1‑184 (MQ=33)
                                                                                                                                                                                               cccTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAAGTCGGCGGTGGATGTAGACGGTGGTACGTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcg                                                                         <  1:545972/201‑1 (MQ=11)
                                                                                                                                                                                                ccTGAGCAACAACGCCCAGGTCACCATCACGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGtc                                                                                                                                            >  2:82120/1‑132 (MQ=21)
                                                                                                                                                                                                 cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTTTACAacg                                                                                                                                                                                        >  2:255292/1‑88 (MQ=18)
                                                                                                                                                                                                 cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                        <  2:1045639/88‑1 (MQ=35)
                                                                                                                                                                                                 cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGATGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                        >  1:1045640/1‑88 (MQ=18)
                                                                                                                                                                                                      caacaaCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                        <  2:162867/115‑1 (MQ=33)
                                                                                                                                                                                                      caacaaCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                        >  1:162867/1‑115 (MQ=33)
                                                                                                                                                                                                       aacaacGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCg                               <  1:1030070/235‑1 (MQ=32)
                                                                                                                                                                                                         caacGCTCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCGGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcgacga                                                                     <  1:58260/195‑1 (MQ=11)
                                                                                                                                                                                                          cacGCCCCGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGAc                                                                                    <  2:154868/178‑1 (MQ=17)
                                                                                                                                                                                                          aacGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGc                                                                                                      <  1:914589/161‑1 (MQ=33)
                                                                                                                                                                                                            cGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGc                                                                                                      >  1:220690/1‑159 (MQ=32)
                                                                                                                                                                                                            cGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAgg                           <  1:461157/234‑1 (MQ=31)
                                                                                                                                                                                                            cGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCCGCTGTGGATGGCGACGGCCGCACCTTCGAGAACCTGCAACTGGGCGGTGCGGCTACCGTTCa                                                                                          >  2:243046/1‑171 (MQ=255)
                                                                                                                                                                                                            cGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCCGCTGTAGATGtcaacg                                                                                                                                        >  2:545971/1‑120 (MQ=255)
                                                                                                                                                                                                              cccAGGTCCCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGc                       <  1:192023/236‑1 (MQ=14)
                                                                                                                                                                                                                 aGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAAGTCCGCTGTAGATGTAGACGGTGGTACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTgg                        >  2:1018141/1‑232 (MQ=9)
                                                                                                                                                                                                                     cACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGAc                                             <  1:698441/207‑1 (MQ=31)
                                                                                                                                                                                                                      aCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGtaa                                  >  1:83499/1‑215 (MQ=14)
                                                                                                                                                                                                                       ccATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCAc               <  1:765184/235‑1 (MQ=31)
                                                                                                                                                                                                                        cATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCAc               <  2:501824/234‑1 (MQ=31)
                                                                                                                                                                                                                         aTCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                              <  2:365470/122‑1 (MQ=17)
                                                                                                                                                                                                                         aTCAAGGACGGCGAAACCAGCGCGCCAGACACCCACGCTGCGCAAGGCGATTACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCg                                                                                                                            <  2:99788/124‑1 (MQ=12)
                                                                                                                                                                                                                           cAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                        <  1:552139/46‑1 (MQ=38)
                                                                                                                                                                                                                             aGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTaca          <  1:168492/234‑1 (MQ=14)
                                                                                                                                                                                                                              ggCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGaa                                                                    <  2:268843/175‑1 (MQ=31)
                                                                                                                                                                                                                               gctggCGATACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGAGGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGAc                                                   <  2:235804/188‑1 (MQ=9)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCgg                                                                                                    <  1:449610/141‑1 (MQ=32)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCgg                                                                                                    >  2:449609/1‑141 (MQ=32)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCCCCGaa                             >  2:510129/1‑212 (MQ=14)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCgag                    >  2:46834/1‑221 (MQ=31)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATc     <  1:228259/236‑1 (MQ=31)
                                                                                                                                                                                                                                    cGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCAcc  <  1:757898/235‑1 (MQ=11)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
CAAGTTCGAACAGCTTACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCACC  >  AE015451/215152‑215614

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