Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 215,413 A→G 8.0% Q6973Q (CAA→CAG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451215,4130AG 8.0% 226.9 / 5.3 83Q6973Q (CAA→CAGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (32/43);  new base G (4/3);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATC  >  AE015451/215182‑215611
                                                                                                                                                                                                                                       |                                                                                                                                                                                                      
gCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                     <  1:380662/235‑1 (MQ=34)
 ccGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                     <  2:632302/234‑1 (MQ=35)
  cGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  <  1:568673/236‑1 (MQ=37)
  cGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGAAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  <  1:734135/236‑1 (MQ=37)
                 tGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  >  1:951449/1‑221 (MQ=35)
                     cgacgaGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                                                                                                                                                          <  1:390659/225‑1 (MQ=37)
                         gaGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGcc                                                                                                                                                                            <  1:1257968/235‑1 (MQ=39)
                             cAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTc                                                                                                                                                                                            >  2:658783/1‑215 (MQ=18)
                               ggTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATc                                                                                                                                                                       <  1:707332/234‑1 (MQ=37)
                                  aCGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTAc                                                                                                                                                                                         >  2:518573/1‑213 (MQ=38)
                                    gCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAg                                                                                                                                                                                                       <  1:136458/197‑1 (MQ=34)
                                     ccGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTg                                                                                                                                                                 <  2:1023206/234‑1 (MQ=37)
                                        ggCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCg                                                                                                                                                                                                    >  2:501814/1‑196 (MQ=25)
                                             cGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTCTACaa                                                                                                                                                                                       >  1:10731/1‑204 (MQ=255)
                                             cGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                     <  1:608252/206‑1 (MQ=21)
                                                    gACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                    <  2:621167/232‑1 (MQ=38)
                                                     aCTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGc                                                                                                                                                   <  2:1254976/232‑1 (MQ=38)
                                                                gTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATc                                                                                                                                                           <  2:1132016/213‑1 (MQ=38)
                                                                   aCCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                    <  1:537699/217‑1 (MQ=38)
                                                                           cgccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTc                                                                                                                          <  1:1097033/235‑1 (MQ=37)
                                                                            gccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAGGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAc                                                                                                                                                         >  1:374097/1‑203 (MQ=25)
                                                                                  cAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAccc                                                                                                                   >  1:484822/1‑235 (MQ=37)
                                                                                  cAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCGGGGCATCACGTCCg                                                                                                                                                    >  2:418207/1‑202 (MQ=25)
                                                                                   aGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATc                                                                                                                                                                       <  1:377486/182‑1 (MQ=37)
                                                                                   aGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATc                                                                                                                                                                       >  2:377481/1‑182 (MQ=37)
                                                                                       cccGGTGGCCGCGACCGCCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTg                                                                                                                                                                 <  2:171007/182‑1 (MQ=11)
                                                                                             ggCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGCTGTAGACGGt                                                                                                                                   >  1:1277383/1‑208 (MQ=25)
                                                                                               ccGAGAACGTCAAGCCGACCTTCACAGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATc                                                                                                                                                                       <  2:1115511/170‑1 (MQ=37)
                                                                                                       gTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAc                                                                                                                                                         <  2:903245/176‑1 (MQ=35)
                                                                                                               gACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                           <  2:935679/198‑1 (MQ=34)
                                                                                                                    tCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAc                                                                                                                                                         >  2:3501/1‑163 (MQ=33)
                                                                                                                      aCCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCCCGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGAcc                                                                                <  1:1100520/234‑1 (MQ=17)
                                                                                                                          tGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGc                                                                                                   <  2:942565/211‑1 (MQ=33)
                                                                                                                              cGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                           <  1:206773/199‑1 (MQ=33)
                                                                                                                              cGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                           <  2:804201/199‑1 (MQ=33)
                                                                                                                                  aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCACGGCGACGACGTCTACAACGACGCTGGCCCGACCAGCCTGGGc                                                                                                                                                              >  2:138233/1‑144 (MQ=255)
                                                                                                                                  aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTGTACAACGACGCTGGCCAGAACAGCCTgg                                                                                                                                                                >  2:18000/1‑142 (MQ=21)
                                                                                                                                     cAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTg                                                                                                                                                                 <  2:426084/138‑1 (MQ=34)
                                                                                                                                       gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAagta                                                              >  1:265729/1‑235 (MQ=255)
                                                                                                                                       gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGAcc                                                                                >  2:539787/1‑217 (MQ=33)
                                                                                                                                       gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAagtag                                                             >  1:536067/1‑236 (MQ=33)
                                                                                                                                         cGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAagtagt                                                            <  1:518578/235‑1 (MQ=32)
                                                                                                                                          gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTa                                                                                              >  1:977233/1‑200 (MQ=34)
                                                                                                                                            tGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGcaa                                                                                                                      <  2:694965/174‑1 (MQ=17)
                                                                                                                                              gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCaa                                                      <  1:1137761/236‑1 (MQ=33)
                                                                                                                                               cccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAg                                                     >  1:913749/1‑236 (MQ=32)
                                                                                                                                                ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGc                                                    >  1:993006/1‑236 (MQ=33)
                                                                                                                                                     aCCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGa                                                 >  2:706515/1‑234 (MQ=33)
                                                                                                                                                      ccTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACcg                                              >  1:640569/1‑236 (MQ=33)
                                                                                                                                                           tCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                                                                                                                                                          <  1:749580/91‑1 (MQ=34)
                                                                                                                                                           tCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                                                                                                                                                          >  2:749575/1‑91 (MQ=34)
                                                                                                                                                             gggACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                           <  2:388468/166‑1 (MQ=14)
                                                                                                                                                             gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGaa                                                                                                                       >  1:388473/1‑156 (MQ=11)
                                                                                                                                                              tGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGAccc                                        <  2:265725/234‑1 (MQ=255)
                                                                                                                                                                  ccTGAGCAACAACGCGCAGGTCACCATCAAGGCCCGCGAAACCAGCGCGCCATACACCCCCGCTGCGCAAGGc                                                                                                                                                                                                     >  2:1147398/1‑73 (MQ=14)
                                                                                                                                                                   cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATc                                                                                                                                                           >  1:716821/1‑114 (MQ=33)
                                                                                                                                                                   cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATc                                                                                                                                                           <  2:716816/114‑1 (MQ=33)
                                                                                                                                                                   cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGCTa                                                                            >  2:183189/1‑193 (MQ=18)
                                                                                                                                                                   cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcgacg                                                                   >  2:429308/1‑202 (MQ=32)
                                                                                                                                                                   cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGt                                 <  1:796334/236‑1 (MQ=33)
                                                                                                                                                                     aggCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGcc                                                                                                                                                                            <  2:213694/93‑1 (MQ=16)
                                                                                                                                                                       gCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGaa                                                                 >  1:1047956/1‑200 (MQ=33)
                                                                                                                                                                        caacaaCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCg                                                                               >  1:904709/1‑185 (MQ=33)
                                                                                                                                                                         aacaacGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCg                            <  1:732017/235‑1 (MQ=32)
                                                                                                                                                                                cccAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTTCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGc                    <  1:116027/236‑1 (MQ=14)
                                                                                                                                                                                 ccAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGc                    <  1:731432/235‑1 (MQ=31)
                                                                                                                                                                                  cAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCAc                              <  1:191345/224‑1 (MQ=31)
                                                                                                                                                                                   aGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGcc                                                        >  2:1178211/1‑197 (MQ=14)
                                                                                                                                                                                      tCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATc              <  1:539792/236‑1 (MQ=31)
                                                                                                                                                                                        aCCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTct                                                                                                                                                     <  2:822204/99‑2 (MQ=9)
                                                                                                                                                                                        aCCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTct                                                                                                                                                     >  1:822209/1‑98 (MQ=9)
                                                                                                                                                                                        aCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATgca                                                                                                                                         >  1:221091/1‑109 (MQ=17)
                                                                                                                                                                                        aCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATc              <  2:556684/234‑1 (MQ=31)
                                                                                                                                                                                        aCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCAc            <  2:640564/236‑1 (MQ=31)
                                                                                                                                                                                        aCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTAACGTTCAAGTGACCGATaccac                                                                        >  2:355131/1‑176 (MQ=14)
                                                                                                                                                                                         ccATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTCCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCCCCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCAc            <  1:703689/235‑1 (MQ=14)
                                                                                                                                                                                             cAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTTCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAag                                                                >  2:116024/1‑179 (MQ=14)
                                                                                                                                                                                             cAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATac                                                                           >  2:1059752/1‑168 (MQ=31)
                                                                                                                                                                                             cAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGc                    <  2:653120/223‑1 (MQ=31)
                                                                                                                                                                                             cAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGc                    >  1:653125/1‑223 (MQ=31)
                                                                                                                                                                                              aaGGCCGGCGAAACCAGCGCGCCATACTCCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTaca       <  1:706244/235‑1 (MQ=14)
                                                                                                                                                                                                  ctggCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGACGCTGGCCAGATCCGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAGCTGGGCGGTGCCGCTACCGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATc  <  1:983208/234‑1 (MQ=11)
                                                                                                                                                                                                   cGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAAGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCga                  >  2:1097839/1‑219 (MQ=14)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                      
GCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATC  >  AE015451/215182‑215611

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