Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 200,113 T→C 8.9% N1873N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451200,1130TC 8.9% 227.9 / 6.9 74N1873N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base C (3/4);  total (36/38)
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.97e-01

CATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCG  >  AE015451/199879‑200348
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           
cagcGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                                                                                                                                             <  1:371178/232‑1 (MQ=37)
 aTCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATg                                                                                                                                                                                                                                            <  1:502093/235‑1 (MQ=12)
     aCGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTc                                                                                                                                                                                                                                         >  2:653017/1‑234 (MQ=40)
     aCGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                                       >  1:1166885/1‑236 (MQ=255)
     aCGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                                       <  1:165772/236‑1 (MQ=40)
             aGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCAt                                                                                                                                                                                                                               >  1:876189/1‑236 (MQ=255)
               cggcgCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCa                                                                                                                                                                                                                             <  1:276907/236‑1 (MQ=40)
                  cgcAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCa                                                                                                                                                                                                                             >  2:813398/1‑233 (MQ=255)
                     aGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCt                                                                                                                                                                                                                        <  1:1261052/235‑1 (MQ=40)
                               aaGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                                                                                                                             >  1:992994/1‑236 (MQ=255)
                               aaGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                                                                                                                             >  1:1254181/1‑236 (MQ=255)
                                 gTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                           <  1:249674/236‑1 (MQ=255)
                                    cgAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCaa                                                                                                                                                                                                                            >  2:888595/2‑216 (MQ=11)
                                       aaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATg                                                                                                                                                                                                                                            >  2:325759/1‑197 (MQ=35)
                                         acatgaacctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCGCAg                                                                                                                                                                                                                                       >  1:706912/4‑200 (MQ=18)
                                              aacctgGACAAGACCCCGGGCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGACGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgc                                                                                                                                                                                                            <  1:216301/222‑1 (MQ=255)
                                               acctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTgg                                                                                                                                                                                                                      <  1:609235/211‑1 (MQ=38)
                                                 ctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCCCTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCa                                                                                                                                                                                                                             <  1:184624/202‑1 (MQ=38)
                                                          aCCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAgg                                                                                                                                                                                  <  2:898695/236‑1 (MQ=255)
                                                              cGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGAt                                                                                                                                                                              <  1:328856/236‑1 (MQ=255)
                                                              cGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGAt                                                                                                                                                                              <  1:618932/236‑1 (MQ=255)
                                                              cGGCCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGAt                                                                                                                                                                              <  1:532741/236‑1 (MQ=255)
                                                                 tgctgaCCGGGGTTACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGAt                                                                                                                                                                                          >  1:901614/6‑221 (MQ=255)
                                                                 tgctgaCCGGGGTTACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGAt                                                                                                                                                                                          <  2:901609/216‑1 (MQ=255)
                                                                      aCCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCt                                                                                                                                                                       <  1:959556/235‑1 (MQ=39)
                                                                          aGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGATCAACAATGATCAGCTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTa                                                                                                                                                                      <  2:693587/232‑1 (MQ=255)
                                                                             tCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGaa                                                                                                                                                                                                                  <  2:816324/185‑1 (MQ=255)
                                                                             tCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGaa                                                                                                                                                                                                                  >  1:816329/1‑185 (MQ=255)
                                                                             tCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGCAAc                                                                                                                                                                                                                >  2:386337/1‑187 (MQ=255)
                                                                             tCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAgg                                                                                                                                                                 <  2:992989/234‑1 (MQ=255)
                                                                                 ccacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCa                                                                                                                                                                                                              <  1:92249/183‑1 (MQ=255)
                                                                                  gacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                                       <  2:906899/159‑1 (MQ=40)
                                                                                  gacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCg                                                                                                                                                                                                      >  1:146576/1‑192 (MQ=40)
                                                                                       gCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCCCCATCCCGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCa                                                                                                                                                                                                                                        <  2:277695/153‑1 (MQ=255)
                                                                                             cACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCa                                                                                                                                                  <  2:488295/233‑1 (MQ=255)
                                                                                               cTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGAt                                                                                                                                                                                          <  2:968335/191‑1 (MQ=40)
                                                                                                  cAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCa                                                                                                                                                  <  1:1030714/228‑1 (MQ=255)
                                                                                                   aGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTgggt                                                                                                                                          <  1:596493/235‑2 (MQ=255)
                                                                                                     gCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAACCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCACGGCTGCCGAAACCAGCGCGCCg                                                                                                                                                                                                      >  2:126106/1‑173 (MQ=255)
                                                                                                             gtctaCATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGTCCTTCACCGTGCACGTCAACCACCCGCTGGCCAACGACCTGGTCGTAACCCTGAGCAACAATACTCAGGTCACAATCAAGGCt                                                                                                                                                                                                                        >  2:1163422/5‑147 (MQ=11)
                                                                                                               cGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                                                                                                                                       >  2:994693/1‑130 (MQ=25)
                                                                                                               cGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                                                                                                                                       <  1:994698/130‑1 (MQ=25)
                                                                                                                           ggTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACACCCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                                                                                                                             >  1:646836/1‑144 (MQ=11)
                                                                                                                               aCCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGtg                                                                                                              >  1:475535/1‑235 (MQ=255)
                                                                                                                                  aTCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGa                                                                                                                         >  2:1196779/1‑221 (MQ=255)
                                                                                                                                  aTCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGtgct                                                                                                            >  2:793031/1‑234 (MQ=255)
                                                                                                                                  aTCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGggctgc                                                                                                          >  1:217142/1‑236 (MQ=255)
                                                                                                                                  aTCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAACCCAGCGCGCCGCACACCCACGCGGCGCAAGGc                                                                                                                                                                                 >  2:376461/1‑165 (MQ=255)
                                                                                                                                     aCGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTg                                                                                                                                             <  1:920047/198‑1 (MQ=40)
                                                                                                                                       ggCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGtt                                                                                                       <  1:1178396/234‑1 (MQ=255)
                                                                                                                                         ccGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACgcgc                                                                                                                                                                                      >  2:1192530/1‑153 (MQ=14)
                                                                                                                                         ccGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAg                                                                                                                                                                                   >  2:1261851/1‑156 (MQ=14)
                                                                                                                                            aCCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAAc                                                                                                >  1:558332/1‑236 (MQ=255)
                                                                                                                                              cAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAAc                                                                                                >  2:601388/1‑234 (MQ=255)
                                                                                                                                              cAGACTTCGGTAGCCGAGGACGTCAAGCCGACCTTCACCATGCACGTCAACCAGCCGCTGACCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGtgctg                                                                                                           >  2:715502/1‑223 (MQ=255)
                                                                                                                                                     cGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgc                                                                                                                                                                                                            <  2:343625/119‑1 (MQ=32)
                                                                                                                                                           ccGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCgg                                                                                  >  1:482948/1‑235 (MQ=255)
                                                                                                                                                               gaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                                                                                                                             >  1:665844/1‑108 (MQ=33)
                                                                                                                                                               gaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                                                                                                                             <  2:665839/108‑1 (MQ=33)
                                                                                                                                                                       aaCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGtt                                                                     <  1:791807/236‑1 (MQ=255)
                                                                                                                                                                                         aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                           >  2:621165/1‑84 (MQ=35)
                                                                                                                                                                                         aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                           <  1:621170/84‑1 (MQ=33)
                                                                                                                                                                                         aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCa                                                                                                                                                                                              >  1:1044426/1‑97 (MQ=33)
                                                                                                                                                                                         aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCa                                                                                                                                                                                              <  2:1044421/97‑1 (MQ=33)
                                                                                                                                                                                         aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCt                                                    >  1:879614/1‑235 (MQ=255)
                                                                                                                                                                                          cGTCAACCAGCCGCTGGCCTACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCa                                                                                                                                                       <  1:219698/135‑1 (MQ=255)
                                                                                                                                                                                                        tGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGcc                                                                                                                                                        <  2:988729/120‑1 (MQ=12)
                                                                                                                                                                                                          gCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGcc                                    >  2:574557/1‑234 (MQ=255)
                                                                                                                                                                                                                 aCCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCAc                                                             <  1:1031850/202‑1 (MQ=9)
                                                                                                                                                                                                                       tCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGcc                                    <  2:732951/221‑1 (MQ=255)
                                                                                                                                                                                                                          tGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCgaga                                                                                                  <  2:754838/156‑1 (MQ=255)
                                                                                                                                                                                                                             cccTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGa                                                                                                                         <  1:992975/130‑1 (MQ=255)
                                                                                                                                                                                                                               cTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGAc                                                                                                                        >  2:755221/1‑129 (MQ=255)
                                                                                                                                                                                                                               cTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCg               >  2:787857/1‑234 (MQ=255)
                                                                                                                                                                                                                                    caacaaTGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCGGCTTCGGt                                                                      >  1:264900/1‑174 (MQ=255)
                                                                                                                                                                                                                                          tGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCg  >  1:464766/1‑236 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           
CATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCG  >  AE015451/199879‑200348

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