Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 212,155 T→C 12.1% G5887G (GGT→GGC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451212,1550TC12.1% 131.6 / 7.9 49G5887G (GGT→GGCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (21/22);  new base C (2/4);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACC  >  AE015451/211928‑212386
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
gAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTc                                                                                                                                                                                                                                 >  1:646659/1‑236 (MQ=32)
       cccTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGt                                                                                                                                                                                                                           <  2:393637/234‑1 (MQ=18)
                                cGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCgaga                                                                                                                                                                                                 <  1:619974/236‑1 (MQ=33)
                                                                            ttCACGTCACCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGacac                                                                                                                                                      <  1:81218/235‑1 (MQ=255)
                                                                             tCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACCCCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTcg                                                                                                                                                                                                                                <  2:130126/160‑1 (MQ=18)
                                                                                    aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGcggcg                                                                                                                                                                       >  1:209950/1‑206 (MQ=25)
                                                                                        aGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGaa                                                                                                                                            <  2:392074/233‑1 (MQ=255)
                                                                                                     cGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCgg                                                                                                                                                                            >  2:533986/1‑188 (MQ=37)
                                                                                                     cGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGa                                                                                                                            <  1:251125/236‑1 (MQ=255)
                                                                                                           ggTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCg                                                                                                                       <  2:515495/235‑1 (MQ=255)
                                                                                                                          caacaaTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCg                                                                                                                 <  1:582514/226‑1 (MQ=255)
                                                                                                                          caacaaTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGAcc                                                                                                        >  1:654713/1‑235 (MQ=255)
                                                                                                                           aacaaTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCCCCGGTGCGACCTTCGCGAACCTGGAGCTGGGCg                                                                                                                                                                                  >  1:520789/1‑160 (MQ=255)
                                                                                                                               aTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAgg                                                                                                   >  1:423135/1‑235 (MQ=255)
                                                                                                                                   ccAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                               <  1:515087/235‑1 (MQ=255)
                                                                                                                                       gTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCTAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTg                                                                                              <  1:365601/232‑1 (MQ=255)
                                                                                                                                        tcgccATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGGCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATc                                                                                         <  1:139653/233‑1 (MQ=255)
                                                                                                                                        tccccATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATc                                                                                         <  1:153053/233‑1 (MQ=255)
                                                                                                                                          aCCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                       >  1:656284/1‑236 (MQ=255)
                                                                                                                                               cAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTcggcggcg                                                                                                                                                                                                                          >  1:67225/1‑98 (MQ=37)
                                                                                                                                                    ccGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGGGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTc                                                                                                                                                                                                                                 >  2:230694/1‑88 (MQ=18)
                                                                                                                                                     tggcgagaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                       <  1:293048/218‑1 (MQ=255)
                                                                                                                                                           aaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTAccc                                                                         <  2:656283/233‑1 (MQ=255)
                                                                                                                                                            aCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                               <  1:480029/210‑1 (MQ=255)
                                                                                                                                                                                  cTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGAccc                                                                                                                   >  1:708042/1‑168 (MQ=9)
                                                                                                                                                                                  cTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCgat                                                                                                                                                                                                  >  2:49731/1‑88 (MQ=18)
                                                                                                                                                                                   tGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                       >  2:609417/1‑195 (MQ=21)
                                                                                                                                                                                   tGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGCCGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGc                                                                                                                                                                                   >  1:75998/1‑103 (MQ=25)
                                                                                                                                                                                   tGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCt                                                                        >  1:478109/1‑210 (MQ=255)
                                                                                                                                                                                   tGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACg                                                          <  1:259841/224‑1 (MQ=255)
                                                                                                                                                                                   tGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACg                                                          >  2:259840/1‑224 (MQ=255)
                                                                                                                                                                                   tGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGGGAGATCACTTACACCATTACCCTGACCCACAAAGACGGTCTGCCGAGTa                                              >  1:168988/1‑236 (MQ=255)
                                                                                                                                                                                   tGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGGGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGCCGGtgga                                                                                                                                                                  >  2:241133/1‑117 (MQ=1)
                                                                                                                                                                                   tGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                       <  1:609418/195‑1 (MQ=255)
                                                                                                                                                                                     cgcAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                               >  1:328348/1‑185 (MQ=255)
                                                                                                                                                                                     cgcAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                               <  2:328347/185‑1 (MQ=255)
                                                                                                                                                                                        aGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCGTCGGTcc                                                                                                                                                                   >  1:645071/1‑114 (MQ=255)
                                                                                                                                                                                                      tACAACGACGGTGGCCAGATCAGCCTGGGCATCAACGCTGCGGTGGACGCCACTAGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                               <  2:708041/168‑1 (MQ=11)
                                                                                                                                                                                                       aCAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCa                                                                                                                                  <  2:369872/132‑1 (MQ=255)
                                                                                                                                                                                                       aCAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCa                                                                                                                                  >  1:369873/1‑132 (MQ=255)
                                                                                                                                                                                                        cAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACg                                                                                                                                                                                                                      >  1:705367/1‑47 (MQ=35)
                                                                                                                                                                                                        cAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACg                                                                                                                                                                                                                      <  2:705366/47‑1 (MQ=35)
                                                                                                                                                                                                            gacgCTGGCCAGATCAGCCTGGGCATCAGCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                       <  1:86389/167‑1 (MQ=255)
                                                                                                                                                                                                                cTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGc                                 <  2:643536/220‑1 (MQ=255)
                                                                                                                                                                                                                          cAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCATTCGGTGACCGAAGGCGGTGAGATCACTTaca                                                                                  <  1:394915/161‑1 (MQ=255)
                                                                                                                                                                                                                           aGCCTGGGTATCAACTCGGCGGTGGGCGCCACTGGTGCGCCCTTCGAGAACCTCGAGCTGGGCGGTGCGGCGGCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTga                                                                                             <  2:165237/149‑1 (MQ=255)
                                                                                                                                                                                                                             ccTGGGTATCAACGCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGacac                                                                                                                                                      <  2:181291/90‑1 (MQ=255)
                                                                                                                                                                                                                                gggTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCAGTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAAcc  >  1:615230/1‑235 (MQ=255)
                                                                                                                                                                                                                                gggTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGGGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGAc                                                                     >  1:200976/1‑168 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
GAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACC  >  AE015451/211928‑212386

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