Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 201,652 G→C 11.8% V2386V (GTG→GTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451201,6520GC11.8% 151.3 / 9.3 50V2386V (GTG→GTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (13/31);  new base C (3/3);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

CAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCG  >  AE015451/201431‑201869
                                                                                                                                                                                                                             |                                                                                                                                                                                                                         
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)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                         
CAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCG  >  AE015451/201431‑201869

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