Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 213,601 G→C 5.4% V6369V (GTG→GTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451213,6010GC 5.4% 312.1 / 4.4 90V6369V (GTG→GTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (36/49);  new base C (2/3);  total (38/52)
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.68e-01

ATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGT  >  AE015451/213375‑213820
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
aTGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAgg                                                                                                                                                                                                                     <  1:693784/235‑1 (MQ=255)
      aGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                 >  2:1077649/1‑233 (MQ=255)
            tcatcaAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                         <  1:975407/235‑1 (MQ=255)
              atcaAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                         >  2:938419/1‑233 (MQ=255)
                caAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACATCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGAcc                                                                                                                                                                                                                          >  1:1045612/1‑214 (MQ=255)
                           aaaCCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGcc                                                                                                                                                                                          >  1:99344/1‑234 (MQ=255)
                           aaaCCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGCTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAggcg                                                                                                                                                                                                                   >  2:623354/1‑210 (MQ=255)
                            aaCCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGa                                                                                                                                                                                        >  1:280766/1‑236 (MQ=255)
                                    cgcCATCCACCCACGCTGCGCAAGGCGAGGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                         <  1:311039/211‑1 (MQ=255)
                                            acCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCt                                                                                                                                                                        <  1:1005284/236‑1 (MQ=255)
                                              ccACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                      <  1:938768/236‑1 (MQ=255)
                                                 cGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAACTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAggcg                                                                                                                                                                                                                   <  1:683531/188‑1 (MQ=255)
                                                    tGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACTGACGCAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACAcc                                                                                                                                                                                                  >  1:35456/1‑202 (MQ=255)
                                                     gcgcAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACcaa                                                                                                                                                                                    <  1:248317/215‑1 (MQ=255)
                                                         aaGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCgag                                                                                                                                                                                                              >  2:445522/1‑185 (MQ=255)
                                                         aaGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATc                                                                                                                                                                                               <  1:215800/200‑1 (MQ=255)
                                                         aaGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATcaacaa                                                                                                                                                            <  1:834190/235‑1 (MQ=255)
                                                                  aCGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                   <  1:512655/235‑1 (MQ=255)
                                                                  aCGTCGACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                         <  2:252557/181‑1 (MQ=255)
                                                                   cGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGcct                                                                                                                                                                    >  1:756072/1‑216 (MQ=255)
                                                                   cGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                         >  2:250648/1‑180 (MQ=255)
                                                                   cGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACcaaca                                                                                                                                                                                  >  2:825398/1‑203 (MQ=255)
                                                                        aCAACGATGCTGGCCAGAGCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGgcggcg                                                                                                                                                                                                                <  2:290607/168‑1 (MQ=255)
                                                                                 cTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCCCGCTg                                                                                                                                   >  1:388445/1‑236 (MQ=255)
                                                                                 cTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCt                                                                                                                                    >  1:941886/1‑235 (MQ=255)
                                                                                 cTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTg                                                                                                                                   <  1:247745/236‑1 (MQ=255)
                                                                                  tGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCGGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGa                                                                                                                                  >  1:241096/1‑236 (MQ=255)
                                                                                         aTCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATc                                                                                                                                                                                               <  2:144542/168‑1 (MQ=255)
                                                                                            aGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                          <  2:941883/234‑1 (MQ=255)
                                                                                                     aTTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATcaac                                                                                                                                                              >  1:240770/1‑189 (MQ=255)
                                                                                                      ttAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACcaca                                                                                                                                                       >  1:648460/1‑195 (MQ=255)
                                                                                                      ttAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGtc                                                                                                                 >  2:715986/1‑233 (MQ=255)
                                                                                                        aactCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAGCTGGGCGGCGCCGCTACGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCaacaac                                                                                                                                                           <  1:233973/186‑1 (MQ=1)
                                                                                                          cgcggcggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGGACCGCTACGGTTGATGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAgg                                                                                                                                                                                                                     <  2:184560/127‑1 (MQ=11)
                                                                                                            cGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTg                                                                                                                                               <  1:309862/197‑1 (MQ=255)
                                                                                                             ggCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                                                                                        <  1:1080238/235‑1 (MQ=255)
                                                                                                             ggCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                                                                                        <  1:546484/235‑1 (MQ=255)
                                                                                                                aGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGATCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                         <  2:661250/135‑1 (MQ=255)
                                                                                                                aGTGGACGCCACTAGTGCGACCTTCTAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                         >  1:661253/1‑135 (MQ=255)
                                                                                                                    gACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGCGACGCGTCGGTGCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                          <  1:704490/178‑1 (MQ=11)
                                                                                                                    gACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCaa                                                                                                <  1:734123/236‑1 (MQ=255)
                                                                                                                         cACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                               <  1:497024/232‑1 (MQ=39)
                                                                                                                         cACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCa                                                                                           <  1:1077652/236‑1 (MQ=255)
                                                                                                                          aCTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACTCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCAc                                                                                                                                                                                             >  2:498489/1‑137 (MQ=255)
                                                                                                                           cTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                          <  1:666272/203‑1 (MQ=40)
                                                                                                                            tGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACcac                                                                                                                                                        >  1:24604/1‑172 (MQ=255)
                                                                                                                            tGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                               >  1:408841/1‑133 (MQ=37)
                                                                                                                              gTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGgcggcg                                                                                                                                                                                                                <  1:659007/114‑1 (MQ=39)
                                                                                                                                      ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATAACCGTGCCGGCCAACGGccc                                                                                          >  1:331740/1‑222 (MQ=25)
                                                                                                                                            aaCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                          <  1:250651/186‑1 (MQ=39)
                                                                                                                                                   aaCTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGCCCAACAAAGACGGCCTGCCGATcaac                                                                                                                                                              >  2:57867/1‑143 (MQ=255)
                                                                                                                                                    aCTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACg                                                                                                                                      >  2:310780/1‑166 (MQ=40)
                                                                                                                                                      tGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGc                                                                                                   <  2:922316/199‑1 (MQ=40)
                                                                                                                                                      tGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                              <  1:557543/236‑1 (MQ=255)
                                                                                                                                                      tGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                              <  1:825401/236‑1 (MQ=255)
                                                                                                                                                         gcggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                              <  2:491112/233‑1 (MQ=255)
                                                                                                                                                           ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCg                                                                                                                                                                    >  1:915774/1‑129 (MQ=255)
                                                                                                                                                           ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCg                                                                                                                                                                    <  2:915771/129‑1 (MQ=255)
                                                                                                                                                           ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAACACCGTCATCACCGTGCCGGACAACGGCAcc                                                                                         >  2:243951/1‑204 (MQ=255)
                                                                                                                                                           ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTcc                                                               >  1:925699/1‑230 (MQ=255)
                                                                                                                                                                ggCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                               <  1:484216/192‑1 (MQ=37)
                                                                                                                                                                ggCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                               >  2:484213/2‑193 (MQ=37)
                                                                                                                                                                cGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGt                                                                        >  2:731337/1‑216 (MQ=255)
                                                                                                                                                                 gCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTa                                                   <  1:715989/236‑1 (MQ=255)
                                                                                                                                                                      ggTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACg                                                 <  2:834726/233‑1 (MQ=255)
                                                                                                                                                                      ggTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTc                                               <  2:1022206/235‑1 (MQ=255)
                                                                                                                                                                              tGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAAc                                      <  1:785240/236‑1 (MQ=255)
                                                                                                                                                                                cccGACCCCACCGACGAAGTAGTGGCCAAGTGGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCGACCCCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGt                                                                                   <  2:353838/188‑1 (MQ=255)
                                                                                                                                                                                    acacCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCt                                >  1:1054913/1‑236 (MQ=14)
                                                                                                                                                                                     caccacCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCt                                <  1:730475/235‑1 (MQ=255)
                                                                                                                                                                                      accaccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCgg                                                             <  1:206488/205‑1 (MQ=40)
                                                                                                                                                                                           cgacgaAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTcc                                                               >  2:146564/1‑198 (MQ=18)
                                                                                                                                                                                           cgacgaAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTg                          <  1:498220/235‑1 (MQ=17)
                                                                                                                                                                                           cgacgaAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTc                                               <  1:19282/214‑1 (MQ=40)
                                                                                                                                                                                              cgaAGTCGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                 <  2:19448/49‑1 (MQ=255)
                                                                                                                                                                                              cgaAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACggt                                                                                                                                                                          >  2:333800/1‑87 (MQ=17)
                                                                                                                                                                                                  gtagtGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGAt                  >  1:106652/1‑236 (MQ=40)
                                                                                                                                                                                                    agtGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATc                 <  1:938422/235‑1 (MQ=255)
                                                                                                                                                                                                      gggCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                     <  2:244386/132‑1 (MQ=21)
                                                                                                                                                                                                         ccAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCTTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACg                                     <  1:577403/210‑1 (MQ=255)
                                                                                                                                                                                                                gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                   >  2:1082722/1‑93 (MQ=37)
                                                                                                                                                                                                                gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                   <  1:1082725/93‑1 (MQ=37)
                                                                                                                                                                                                                gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAAc                                                        >  1:898306/1‑184 (MQ=255)
                                                                                                                                                                                                                gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAAc                                                        <  2:898303/184‑1 (MQ=40)
                                                                                                                                                                                                                    gcgACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGt  >  2:935677/1‑234 (MQ=40)
                                                                                                                                                                                                                      gACCCCGACGGTGACCGAAGGCGGCGAGAATACCTACCCCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTa                                                   <  2:280763/183‑1 (MQ=255)
                                                                                                                                                                                                                           cGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACACCGTCTACGTCGGTACCAACGACCCTGTCGCGAGGt                      >  2:346385/1‑207 (MQ=255)
                                                                                                                                                                                                                            gTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCa                                                                                                                                                                                              >  2:70273/1‑38 (MQ=40)
                                                                                                                                                                                                                             tCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGc                                                                                                                                                                                           <  2:747360/40‑1 (MQ=40)
                                                                                                                                                                                                                              cGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGc                                                                                                                                                                                           <  1:197535/39‑1 (MQ=40)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
ATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGT  >  AE015451/213375‑213820

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