Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,360 T→C 9.4% A1622A (GCT→GCC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451199,3600TC 9.4% 266.0 / 6.5 90A1622A (GCT→GCCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (37/44);  new base C (4/5);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTG  >  AE015451/199130‑199580
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                            
gtcgtcGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCa                                                                                                                                                                                                                         <  2:1003174/236‑1 (MQ=35)
   gtcgtGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGa                                                                                                                                                                                                                       >  1:644374/1‑235 (MQ=35)
                       caccGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                                                                                   <  1:554466/235‑1 (MQ=35)
                            gTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAAc                                                                                                                                                                                                      >  1:50252/1‑227 (MQ=17)
                                     cTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCaa                                                                                                                                                                                                 <  1:220497/223‑1 (MQ=35)
                                         aGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAAc                                                                                                                                                                                                      >  1:332731/1‑214 (MQ=35)
                                         aGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGGCAAGc                                                                                                                                                                                               >  2:634900/1‑221 (MQ=37)
                                               cccGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCa                                                                                                                                                                          <  1:469433/236‑1 (MQ=255)
                                                  ggACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGc                                                                                                                                                                                                             <  1:1061256/198‑1 (MQ=18)
                                                   gACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAg                                                                                                                                                                                                                        >  2:1035148/1‑186 (MQ=38)
                                                   gACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTg                                                                                                                                                                                 >  1:1059248/1‑225 (MQ=40)
                                                   gACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAg                                                                                                                                                                                                                        >  2:854159/1‑186 (MQ=18)
                                                     cAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGCCAACCAGc                                                                                                                                                                    >  1:1069871/1‑236 (MQ=255)
                                                       aCGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGcc                                                                                                                                                                   <  1:50006/235‑1 (MQ=255)
                                                           cTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTg                                                                                                                                                               <  1:504641/235‑1 (MQ=255)
                                                           cTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTg                                                                                                                                                               >  1:1087521/1‑235 (MQ=255)
                                                                tCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                                                                                   <  1:250418/194‑1 (MQ=37)
                                                                 cGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                                                                                   <  2:405274/193‑1 (MQ=37)
                                                                  ggCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACg                                                                                                                                                        <  1:712006/235‑1 (MQ=255)
                                                                  ggCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACg                                                                                                                                                        <  1:621608/235‑1 (MQ=255)
                                                                     accaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGc                                                                                                                                                                                               >  1:818334/1‑193 (MQ=21)
                                                                     accaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGAc                                                                                                                                                                                                                           >  2:684465/1‑165 (MQ=38)
                                                                     accaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCt                                                                                                                                                    >  1:925773/1‑236 (MQ=255)
                                                                     accaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCCGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTc                                                                                                                                                                                                                   >  2:170352/1‑173 (MQ=38)
                                                                      ccaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCAc                                                                                                                                                                              <  1:519137/209‑1 (MQ=255)
                                                                      ccaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCAc                                                                                                                                                                              >  2:519134/1‑209 (MQ=40)
                                                                      ccaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCaa                                                                                                                                                                         >  2:970983/1‑214 (MQ=255)
                                                                       caaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCg                                                                                                                                                                                   >  2:381326/1‑203 (MQ=40)
                                                                       caaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                                                                  >  1:424544/1‑236 (MQ=255)
                                                                        aaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGt                                                                                                                                                 <  1:731633/236‑1 (MQ=255)
                                                                        aaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTc                                                                                                                                                                                                   <  1:854162/186‑1 (MQ=18)
                                                                           caGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGGCg                                                                                                                                               <  1:1074385/235‑1 (MQ=39)
                                                                            aGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCGTCACCGTGc                                                                                                                                                                                >  2:294159/1‑201 (MQ=255)
                                                                                ggtgggCAATGCCAGCGACGCGGTCAGCGGCGCAGATGCAGGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTc                                                                                                                                                                                                   <  2:50251/178‑1 (MQ=255)
                                                                                     tCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACg                                                                                                                                                        >  2:478605/1‑216 (MQ=39)
                                                                                     tCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACg                                                                                                                                                        <  1:478608/216‑1 (MQ=40)
                                                                                          gCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGCCATCGTCAAGGTCACCATCACGGCTGACCAGACTTCgg                                                                                                                                                                                                                 >  2:305767/1‑154 (MQ=25)
                                                                                          gCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGca                                                                                                                                   <  2:1087518/232‑1 (MQ=38)
                                                                                           ccATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                                                                                                                                             >  2:1037700/1‑221 (MQ=21)
                                                                                                 aCGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCt                                                                                                                                                    >  2:330241/1‑208 (MQ=38)
                                                                                                 aCGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCt                                                                                                                                                    >  2:1011824/1‑208 (MQ=38)
                                                                                                    cgGTCAGCGGCGCAGATGCATGGAAGTTCTAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCATGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                                                                                                                                             <  2:1028075/212‑1 (MQ=9)
                                                                                                           cggcgCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                                                                              <  1:1027599/236‑1 (MQ=255)
                                                                                                           cggcgCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                                                                              <  1:381329/236‑1 (MQ=40)
                                                                                                                    tGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGcc                                                                                                                                                                   <  2:1066872/174‑1 (MQ=40)
                                                                                                                      cATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaacaatg                                                                                                                             >  2:381415/1‑208 (MQ=18)
                                                                                                                           aaGTTCGAAAACCTGAACCTGGCCAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCa                                                                                               <  2:504462/235‑1 (MQ=255)
                                                                                                                                cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGc                                                                                                                                                                    >  2:179656/1‑161 (MQ=37)
                                                                                                                                   aaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaaca                                                                                                                                <  2:468268/194‑1 (MQ=35)
                                                                                                                                   aaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaaca                                                                                                                                >  1:468271/1‑194 (MQ=35)
                                                                                                                                     acctgaacctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGt                                                                                      <  2:925770/234‑1 (MQ=255)
                                                                                                                                                    aGACCCCGGTCACCACCGAGGTCACCGACGAGCCAGGAACTCCAGGCAACGAAGGTGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATc                                                                                                                                                <  2:577252/161‑1 (MQ=255)
                                                                                                                                                        cccGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAg                                                                    >  1:40451/1‑233 (MQ=17)
                                                                                                                                                        cccGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGCCACCATCACGGCTGACCAg                                                                                                                                                                                                                        >  2:79395/1‑85 (MQ=18)
                                                                                                                                                          cGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGCCGCGCAAGGCGAt                                                               <  1:661305/236‑1 (MQ=255)
                                                                                                                                                          cGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCg                                                                 <  2:424541/234‑1 (MQ=255)
                                                                                                                                                           ggTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                                                                                                                                             <  1:485289/157‑1 (MQ=35)
                                                                                                                                                             tCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCa                                                                                                                                                                          <  2:969897/126‑1 (MQ=255)
                                                                                                                                                                    cGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGCCAAGGTCACCATCACGGCTGACCAGACTGCGGTAGCCGAGCACGCCAAGCCGACCATCACCGTg                                                                                                                                                                                 >  2:1046225/1‑112 (MQ=14)
                                                                                                                                                                       ggTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGCCATCGTCACGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAgcg                                                                                            <  1:114837/194‑1 (MQ=255)
                                                                                                                                                                        ggcACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGCCATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAgg                                                  <  1:298282/233‑1 (MQ=255)
                                                                                                                                                                        gTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAgg                                                  <  1:768889/235‑1 (MQ=255)
                                                                                                                                                                         tCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTgg                                                                                                       <  2:405753/181‑1 (MQ=40)
                                                                                                                                                                           aCCGACGTGCCAAGCAGTCCAGGGGACGAAGGGGAGTTCGTAAAGGACAGCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGa                                                                                                                                                                                            <  1:876585/94‑1 (MQ=14)
                                                                                                                                                                              gacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTgg                                                                                                       <  1:179658/176‑1 (MQ=39)
                                                                                                                                                                              gacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCa                                                                                               <  1:526332/184‑1 (MQ=255)
                                                                                                                                                                              gacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCa                                                                                               >  2:526329/1‑184 (MQ=39)
                                                                                                                                                                              gacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCt                                             <  1:842321/234‑1 (MQ=255)
                                                                                                                                                                              gacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACACGCCGTTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACAcc                                                                                 >  1:163372/1‑198 (MQ=255)
                                                                                                                                                                                    ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAACGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                                                                              <  2:432584/163‑1 (MQ=39)
                                                                                                                                                                                     cAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                                                                              <  2:818331/162‑1 (MQ=39)
                                                                                                                                                                                         tacgccAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                                                                                                                     >  2:1023679/5‑135 (MQ=255)
                                                                                                                                                                                         cacgccAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                                                                                                                     <  1:1023682/131‑1 (MQ=255)
                                                                                                                                                                                            tCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACg                                                                                                                                                        <  1:691911/113‑1 (MQ=40)
                                                                                                                                                                                            gccAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACg                                                                                                                                                        >  2:691908/2‑113 (MQ=255)
                                                                                                                                                                                                ggCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaaca                                                                                                                                >  2:1055736/1‑133 (MQ=38)
                                                                                                                                                                                                    gagAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGccc                                                                                                                          <  1:170334/133‑1 (MQ=39)
                                                                                                                                                                                                         ggCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTggcg                                                                                                     >  1:540490/1‑149 (MQ=39)
                                                                                                                                                                                                         ggCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAgg                                                                   <  2:451727/185‑1 (MQ=40)
                                                                                                                                                                                                         ggCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAgg                                                                   >  1:451730/1‑185 (MQ=255)
                                                                                                                                                                                                         ggCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGAc                                                >  1:172765/1‑204 (MQ=255)
                                                                                                                                                                                                         ggCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCacg                                                                             >  2:779194/1‑175 (MQ=255)
                                                                                                                                                                                                              cATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCt                                                                                                                                                                                         <  2:388808/62‑1 (MQ=255)
                                                                                                                                                                                                              cATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCt                                                                                                                                                                                         >  1:388811/1‑62 (MQ=255)
                                                                                                                                                                                                              cATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTg            <  1:693682/235‑1 (MQ=255)
                                                                                                                                                                                                                  gTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGca                                                                                                                                   <  2:565253/112‑1 (MQ=255)
                                                                                                                                                                                                                  gTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGca                                                                                                                                   >  1:565256/1‑112 (MQ=255)
                                                                                                                                                                                                                        gTCACCATCACGGCTGCCCAGACTTCGGTAGCCGAGACCGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTg  <  1:287391/235‑1 (MQ=255)
                                                                                                                                                                                                                              aTCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACg                                                           >  2:439659/1‑172 (MQ=16)
                                                                                                                                                                                                                                 aCGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGc                                                                                                          <  1:55083/122‑1 (MQ=40)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                            
GTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTG  >  AE015451/199130‑199580

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: