Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 201,649 G→A 10.1% S2385S (TCG→TCA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451201,6490GA10.1% 137.6 / 6.0 49S2385S (TCG→TCAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (14/30);  new base A (2/3);  total (16/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCG  >  AE015451/201416‑201869
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                            
cTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCg                                                                                                                                                                                                                              >  1:549199/1‑234 (MQ=255)
               cAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGgcggcg                                                                                                                                                                                                              <  1:259943/235‑1 (MQ=255)
                aGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCgag                                                                                                                                                                                                            <  1:85331/236‑1 (MQ=255)
                aGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCAACCCCTTCGGGGACCGAAGgcggcg                                                                                                                                                                                                              >  1:240188/1‑234 (MQ=255)
                   tCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCg                                                                                                                                                                                                                       <  2:352400/222‑1 (MQ=255)
                     aCCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGa                                                                                                                                                                                                                          >  1:577411/1‑217 (MQ=255)
                               ctggcgAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                             <  1:544509/231‑1 (MQ=255)
                                               gcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCGTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGCGACCGAAGGCGGCGAGATCACCTACACCATCACGCt                                                                                                                                                                                        >  1:346279/1‑225 (MQ=255)
                                                        ccACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCt                                                                                                                                                                      <  2:408873/234‑1 (MQ=255)
                                                              cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGCTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCGACACCATCACGCTGa                                                                                                                                                                                      >  1:206540/1‑212 (MQ=255)
                                                                tgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTc                                                                                                                                                                       <  2:598244/224‑1 (MQ=255)
                                                                   aaGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACAcc                                                                                                                                                                                                >  2:622639/1‑197 (MQ=255)
                                                                     ggCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                             <  2:514340/198‑1 (MQ=255)
                                                                      gCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCAc                                                                                                                                                      <  1:246099/236‑1 (MQ=255)
                                                                       ccaTGCCGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTg                                                                                                                                                                     <  2:531403/218‑1 (MQ=17)
                                                                       cGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACt                                                                                                                                                     <  1:238296/236‑1 (MQ=255)
                                                                         aTGACGTCGATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCg                                                                                                                                                   <  1:40186/236‑1 (MQ=255)
                                                                                  aTCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                             <  2:705871/185‑1 (MQ=255)
                                                                                                  gATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGAcc                                                                                                                              <  2:680940/232‑1 (MQ=255)
                                                                                                                            tGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTAGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGCCCAACAAAGACGGGCTGCCGATCAATAAcc                                                                                                                                                        <  2:540275/180‑1 (MQ=255)
                                                                                                                              taCGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTg                                                                                                                                                                     <  2:2026/164‑1 (MQ=255)
                                                                                                                               aCGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTac                                                                                                                                                                                                   >  1:44360/1‑134 (MQ=255)
                                                                                                                                            caccttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCg                                                                                                                                                   <  2:565837/165‑1 (MQ=255)
                                                                                                                                              cGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGgcggcg                                                                                                                                                                                                              <  1:711121/108‑1 (MQ=255)
                                                                                                                                                ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGc                                                                                                                                                                                         >  2:103282/1‑127 (MQ=255)
                                                                                                                                                ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGc                                                                                                                                                                                         <  1:103282/127‑1 (MQ=255)
                                                                                                                                                ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTCAcc                                                                                                                                                                                                                        >  2:611512/1‑96 (MQ=18)
                                                                                                                                                     gaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTg                                                                          <  2:648541/233‑1 (MQ=255)
                                                                                                                                                           gCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCCTCGTGGCAGCCAACAGCACCACCGGTTCgg                                                                                >  1:40415/1‑221 (MQ=255)
                                                                                                                                                                   gcggcAACGCTTCGGTTCAAGTGACCGACACCCGCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAAt                                                                                                                                                            <  2:16091/137‑1 (MQ=255)
                                                                                                                                                                     ggcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCaa                                                                                                                                     >  1:105781/1‑158 (MQ=255)
                                                                                                                                                                     ggcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCaa                                                                                                                                     <  2:105781/158‑1 (MQ=255)
                                                                                                                                                                                   tCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCa                                         <  1:55473/236‑1 (MQ=255)
                                                                                                                                                                                    cAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTg                                                                                                         <  1:88683/171‑1 (MQ=255)
                                                                                                                                                                                               cacCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCggt                               >  2:638796/1‑234 (MQ=11)
                                                                                                                                                                                               cacCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                >  2:436778/1‑201 (MQ=255)
                                                                                                                                                                                               cacCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                <  1:436779/201‑1 (MQ=255)
                                                                                                                                                                                                acCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCggtgg                             <  1:597446/235‑1 (MQ=255)
                                                                                                                                                                                                           aGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGc                      >  2:472886/1‑231 (MQ=255)
                                                                                                                                                                                                               gTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                >  1:460625/1‑185 (MQ=255)
                                                                                                                                                                                                               gTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGAc               <  2:309364/234‑1 (MQ=255)
                                                                                                                                                                                                               gTGGCGAAGCTGACCGCCACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                <  2:460624/185‑1 (MQ=37)
                                                                                                                                                                                                                 ggCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCaa                                                                                                 >  2:474796/1‑150 (MQ=255)
                                                                                                                                                                                                                 ggCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCaa                                                                                                 <  1:474797/150‑1 (MQ=255)
                                                                                                                                                                                                                      aGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTg                                                                          <  2:549198/168‑1 (MQ=39)
                                                                                                                                                                                                                         tGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcggc   <  1:311806/236‑1 (MQ=40)
                                                                                                                                                                                                                            ccGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcggcg  <  2:722727/234‑1 (MQ=255)
                                                                                                                                                                                                                                   cccTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACaa                                                          >  2:315129/1‑171 (MQ=34)
                                                                                                                                                                                                                                        cGGTGACCGAAGGCGGCGAGATCACCTACACCTTCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGAc               <  2:397683/209‑1 (MQ=25)
                                                                                                                                                                                                                                        cGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCATGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCg                                                               <  1:552707/161‑1 (MQ=14)
                                                                                                                                                                                                                                         ggTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAAcc                                                                                                                                                        <  1:653856/71‑1 (MQ=34)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                            
CTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCG  >  AE015451/201416‑201869

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