Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,120 T→C 7.8% D1542D (GAT→GAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451199,1200TC 7.8% 208.1 / 3.1 73D1542D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (4/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCAC  >  AE015451/198897‑199350
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      
gCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGACAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACCCCGTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTgt                                                                                                                                                                                                                             >  1:872913/1‑235 (MQ=255)
   tCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGGCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTgtcgtc                                                                                                                                                                                                                         >  1:210256/1‑236 (MQ=255)
        tCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGc                                                                                                                                                                                                                    >  1:904156/1‑236 (MQ=255)
            cggTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGcc                                                                                                                                                                                                                <  1:565750/236‑1 (MQ=255)
            cggTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGCTGGGCGGTGCCGCGTCAGTGGAAGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGc                                                                                                                                                                                                                 >  1:1274303/1‑235 (MQ=11)
             ggTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCa                                                                                                                                                                                                               >  1:1073470/1‑236 (MQ=255)
                ggACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACAGCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCaa                                                                                                                                                                                                              <  2:869386/234‑1 (MQ=255)
                           gTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGt                                                                                                                                                                                                 <  1:436480/236‑1 (MQ=255)
                                 cGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACTAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACGGc                                                                                                                                                                                                          >  2:1185601/1‑221 (MQ=255)
                                   ttCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGc                                                                                                                                                                                            <  2:1073465/233‑1 (MQ=255)
                                   ttCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGc                                                                                                                                                                                            >  2:896895/1‑233 (MQ=255)
                                   ttCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCa                                                                                                                                                                                          >  1:288735/1‑235 (MQ=255)
                                   ttCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCAc                                                                                                                                                                                         >  1:355778/1‑236 (MQ=255)
                                      gagaACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGcc                                                                                                                                                                                           >  1:60391/1‑231 (MQ=255)
                                               gAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTc                                                                                                                                                                                               <  2:586680/218‑1 (MQ=255)
                                                 gTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGt                                                                                                                                                                                      <  2:553175/225‑1 (MQ=255)
                                                 gTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                                                                                                                             >  2:1264720/1‑234 (MQ=255)
                                                 gTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                                                                                                                             <  1:321473/234‑1 (MQ=255)
                                                 gTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCg                                                                                                                                                                            >  1:622292/1‑235 (MQ=255)
                                                 gTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCg                                                                                                                                                                            <  1:490661/235‑1 (MQ=255)
                                                          tGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCt                                                                                                                                                                  <  1:1246759/236‑1 (MQ=255)
                                                            ccGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCt                                                                                                                                                                  >  2:1115953/1‑234 (MQ=255)
                                                             cgcgAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGAGGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCAc                                                                                                                                                                                         >  1:337776/1‑210 (MQ=255)
                                                               cgAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAAc                                                                                                                                                                      <  1:971906/227‑1 (MQ=255)
                                                               cgAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTAt                                                                                                                                                                >  2:1072684/1‑233 (MQ=255)
                                                                   aGTGGATGTCACCGACCCCCTCGACGAGGTGGTCGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACCGTCGTCGTGGCa                                                                                                                                                                                                                   <  1:1012047/178‑1 (MQ=255)
                                                                       gATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAgccgcc                                                                                                                                                                               >  1:495352/1‑210 (MQ=255)
                                                                          gTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCg                                                                                                                                                            >  2:695913/1‑226 (MQ=255)
                                                                                     ccgcgacgaGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGca                                                                                                                                                                                                         <  1:357926/167‑1 (MQ=255)
                                                                                            gaGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCAt                                                                                                                                >  1:984795/1‑236 (MQ=255)
                                                                                               gtggtgGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                                                                                                                             <  2:288731/188‑1 (MQ=40)
                                                                                                tggtggCCAAACTGACTGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAgccg                                                                                                                                                                                 <  1:1126773/183‑1 (MQ=18)
                                                                                                      ccAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCggtg                                                                                                                                           >  2:293918/1‑215 (MQ=39)
                                                                                                      ccAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCAt                                                                                                                                >  2:590382/1‑226 (MQ=40)
                                                                                                           cTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACGCACGGCACTACTCTCGTCGTGACAGCCAACAGCACCACCAGt                                                                                                                                                                                                 >  2:634617/1‑156 (MQ=14)
                                                                                                              ccTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTCCCCCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCg                                                                                                                                                                            <  2:354168/173‑1 (MQ=255)
                                                                                                                   cACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATc                                                                                                                               <  2:807594/214‑1 (MQ=38)
                                                                                                                   cACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACTACTGTCGTCGTGGCAGc                                                                                                                                                                                                                 >  1:1154801/1‑132 (MQ=255)
                                                                                                                   cACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGc                                                                                                                                                                                                                 <  2:1154796/132‑1 (MQ=12)
                                                                                                                       ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGc                                                                                                                                                                                                          >  1:605505/1‑135 (MQ=34)
                                                                                                                       ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGc                                                                                                                                                                                                          <  2:605500/135‑1 (MQ=35)
                                                                                                                         gTCGGTTACCGAAGGCGGTGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAgcggcg                                                                                                               <  1:546481/224‑1 (MQ=11)
                                                                                                                          ccGGTCACCGAAGCCGGCGAGATCACCTCCACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACGGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCacca                                                                                                                                                      <  2:397967/183‑1 (MQ=9)
                                                                                                                            ggTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGt                                                                                                <  1:1156055/236‑1 (MQ=34)
                                                                                                                              tCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAgcgg                                                                                                                 <  2:725041/217‑1 (MQ=33)
                                                                                                                                aCCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGAc                                                                                                                            >  1:981833/1‑204 (MQ=33)
                                                                                                                                aCCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGAc                                                                                                                            <  2:981828/204‑1 (MQ=17)
                                                                                                                                             agaTCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGCCTATGTCGGCACCAACCAGCCGGCGGTCAATGCCATCGACGCGGTCAGCGGCGCAGa                                                                                                           >  1:179105/1‑208 (MQ=14)
                                                                                                                                               aTCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTgg                                                                              <  1:36671/235‑1 (MQ=31)
                                                                                                                                                tCACCTCCACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAaacctgaacct                                                                                <  2:451009/233‑1 (MQ=14)
                                                                                                                                                   ccTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAg                                                                         <  1:590387/236‑1 (MQ=31)
                                                                                                                                                        acCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                                                                    <  1:1026273/236‑1 (MQ=31)
                                                                                                                                                         ccATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCTATGGCACTACTGTCGTCGTTGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGACGCAGATGCATGGAAGTTCGAAAAGCTGAACCTGGACAAGAcccc                                                                    <  2:678396/235‑1 (MQ=9)
                                                                                                                                                           aTCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGc                                                                                                                                                >  2:1028883/1‑157 (MQ=31)
                                                                                                                                                           aTCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGc                                                                                                                                                <  1:1028888/157‑1 (MQ=31)
                                                                                                                                                           aTCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGt                                                                                                                       >  2:1056977/1‑182 (MQ=31)
                                                                                                                                                            tCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTc                                                                <  1:1264725/236‑1 (MQ=31)
                                                                                                                                                               cGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCa                                                                                                             >  2:448180/1‑188 (MQ=31)
                                                                                                                                                               cGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAg                                                              >  1:228320/1‑235 (MQ=11)
                                                                                                                                                                gCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGaa                                                                                                  >  2:577881/1‑198 (MQ=11)
                                                                                                                                                                                        gATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCa                                    <  1:1115958/236‑1 (MQ=11)
                                                                                                                                                                                           cAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTc                                 >  1:407694/1‑236 (MQ=32)
                                                                                                                                                                                               aaCCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAg                              >  1:1033157/1‑235 (MQ=32)
                                                                                                                                                                                                     tCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCaa                          <  2:1152219/233‑1 (MQ=33)
                                                                                                                                                                                                       ggAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCa                                    <  1:695918/221‑1 (MQ=11)
                                                                                                                                                                                                         aGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCa                           <  2:1104451/228‑1 (MQ=33)
                                                                                                                                                                                                         aGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAg                     >  2:443577/1‑234 (MQ=32)
                                                                                                                                                                                                               aCTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGa                 >  1:410941/1‑232 (MQ=32)
                                                                                                                                                                                                                 ttCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGt           >  1:1273833/1‑236 (MQ=33)
                                                                                                                                                                                                                   cAAGATGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTc          <  1:1037884/235‑1 (MQ=18)
                                                                                                                                                                                                                       gCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCCGGCACTCCAGGCAACGAAGGCGACATCGTCAAg       >  1:418290/1‑235 (MQ=18)
                                                                                                                                                                                                                       gCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAg       >  1:785937/1‑235 (MQ=33)
                                                                                                                                                                                                                            ccGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCAc  >  1:557116/1‑235 (MQ=33)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      
GCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCAC  >  AE015451/198897‑199350

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