Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 206,821 T→C 9.6% G4109G (GGT→GGC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451206,8210TC 9.6% 183.6 / 2.6 68G4109G (GGT→GGCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (4/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACC  >  AE015451/206594‑207053
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
gacgaGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGt                                                                                                                                                                                                                                  >  1:948963/1‑236 (MQ=33)
gacgaGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGt                                                                                                                                                                                                                                  >  1:93264/1‑236 (MQ=32)
gacgaGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGAc                                                                                                                                                                                                                                    >  2:915312/1‑234 (MQ=18)
                  ggCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCa                                                                                                                                                                                                                <  1:696792/236‑1 (MQ=32)
                             cggcAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGc                                                                                                                                                                                                                      <  2:880954/219‑1 (MQ=31)
                              cgcagcaACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGAc                                                                                                                                                                                                                                    <  1:595283/203‑1 (MQ=14)
                                       gAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACt                                                                                                                                                                                            >  1:946893/1‑235 (MQ=31)
                                             gACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                                                                                     <  1:280842/236‑1 (MQ=31)
                                              aCCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                                                                                     <  1:651750/235‑1 (MQ=31)
                                                                tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCCTCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAac                                                                                                                                                                                                                           >  1:53353/1‑179 (MQ=14)
                                                                tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGAcc                                                                                                                                                                  >  1:1067767/1‑236 (MQ=31)
                                                                tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGAc                                                                                                                                                                   >  1:976345/1‑235 (MQ=31)
                                                                tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTCCGCAGGGTGATGACGTCTa                                                                                                                                                                                                                               >  1:1087713/1‑175 (MQ=14)
                                                                tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCa                                                                                                                                                                                                                >  1:538586/1‑190 (MQ=9)
                                                                     gCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGTTGCGACCTTCg                                                                                                                                                              >  1:492389/1‑235 (MQ=14)
                                                                           cAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGACGCTGGCCAGAt                                                                                                                                                                                                             >  1:966239/1‑182 (MQ=14)
                                                                                ttCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTg                                                                                                                                                      <  2:673070/232‑1 (MQ=9)
                                                                                gtCGGGGGCCGAGAACGTCAACCCGACGTTCCCCGTACACGTCACCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTa                                                                                                                                                                                                                               <  2:180244/158‑1 (MQ=9)
                                                                                              aCGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCcgcg                                                                                                                                    <  1:915313/236‑1 (MQ=35)
                                                                                                gTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGaa                                                                                                                                  <  1:1128596/236‑1 (MQ=35)
                                                                                                    aaCCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAgg                                                                                                                                <  1:714864/234‑1 (MQ=37)
                                                                                                                                gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCg                                                                                                  >  2:594420/1‑236 (MQ=37)
                                                                                                                                   gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCg                                                                                                  >  2:1115057/1‑233 (MQ=37)
                                                                                                                                       gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCCACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGACCCTGGCGCTGGGCGGTGCCGCGCAGGTggg                                                                                                                            >  2:86189/1‑202 (MQ=255)
                                                                                                                                       gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCaa                                                                                           <  1:841042/236‑1 (MQ=35)
                                                                                                                                       gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCa                                                                                            >  1:1082907/1‑235 (MQ=35)
                                                                                                                                       gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGAGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCaa                                                                                           <  1:898103/236‑1 (MQ=25)
                                                                                                                                          cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCaa                                                                                           <  2:948962/233‑1 (MQ=35)
                                                                                                                                            cGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCt                                                                                                                                                 <  1:33742/177‑1 (MQ=33)
                                                                                                                                            cGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGtt                                                                                        <  2:976344/234‑1 (MQ=35)
                                                                                                                                                     ccagACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTgg                                                                                                                             <  1:313124/185‑1 (MQ=37)
                                                                                                                                                      gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCCTCCAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCAcc                                                                                                                     >  2:494161/1‑195 (MQ=18)
                                                                                                                                                            cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGATGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTg                                                                                                                              >  2:627242/1‑180 (MQ=39)
                                                                                                                                                                       cGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAg                                                                                                     >  2:425770/1‑194 (MQ=255)
                                                                                                                                                                       cGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAg                                                                                                     <  1:425771/194‑1 (MQ=255)
                                                                                                                                                                          ccAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTg                                                                                                                                                >  2:855718/1‑148 (MQ=34)
                                                                                                                                                                          ccAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc                                                         >  1:883410/1‑235 (MQ=35)
                                                                                                                                                                               tCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATc                                                   <  1:1115058/236‑1 (MQ=34)
                                                                                                                                                                                 aCCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAAGTCGGCGGTGGACGCCACCGGTGCGACCTTCGAGAACCTGCCACTGGGCGGTGCcg                                                                                                                                      >  1:549203/1‑151 (MQ=255)
                                                                                                                                                                                 aCCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCaa                                                 >  2:782575/1‑235 (MQ=21)
                                                                                                                                                                                 aCCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                 >  1:83780/1‑236 (MQ=33)
                                                                                                                                                                                  ccATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                 <  1:1002981/235‑1 (MQ=34)
                                                                                                                                                                                      cAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGaa                                                                                                                                  >  2:34826/1‑150 (MQ=37)
                                                                                                                                                                                      cAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGa                                                                >  2:959938/1‑216 (MQ=34)
                                                                                                                                                                                      cAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTaca                                            >  1:495676/1‑236 (MQ=35)
                                                                                                                                                                                         ggCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTaca                                            <  2:1067766/233‑1 (MQ=37)
                                                                                                                                                                                            tGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGCCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc                                                         <  1:447065/217‑1 (MQ=21)
                                                                                                                                                                                               cGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGc                                   <  1:1036415/236‑1 (MQ=33)
                                                                                                                                                                                                gagcCCAGCGCGCCGTACCCCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                 <  1:904356/221‑1 (MQ=18)
                                                                                                                                                                                                gagaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGTGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCAccc                                                                             >  1:235805/1‑193 (MQ=21)
                                                                                                                                                                                                gagaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTGCGAGAACCTGGCGCGGGGCGGt                                                                                                                                          >  2:8717/1‑132 (MQ=21)
                                                                                                                                                                                                 agaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAgg                                                                                                                                >  2:767449/1‑141 (MQ=35)
                                                                                                                                                                                                    ccAGCGCGCCATACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCgag                                                      <  1:24673/212‑1 (MQ=17)
                                                                                                                                                                                                      aGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTg                                 <  2:495675/231‑1 (MQ=32)
                                                                                                                                                                                                           gcCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGGAGGTGGATGTCACCGACACCACCGACGAAGTCgcg                                                                                                >  1:125659/1‑161 (MQ=18)
                                                                                                                                                                                                              gTACACCCACGCTGCGCAGGGTGATGACGTGTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACg                    <  1:61159/236‑1 (MQ=17)
                                                                                                                                                                                                                      aCGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCGCTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                 <  2:70539/199‑1 (MQ=18)
                                                                                                                                                                                                                      aCGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATc                                       <  2:693014/209‑1 (MQ=33)
                                                                                                                                                                                                                      aCGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATc                                       >  1:693015/1‑209 (MQ=32)
                                                                                                                                                                                                                      aCGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGACGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCa            <  1:1070114/236‑1 (MQ=12)
                                                                                                                                                                                                                            cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCcg                                                                                                                                      >  2:1113825/1‑108 (MQ=38)
                                                                                                                                                                                                                            cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAg                                                                                                     >  2:591787/1‑141 (MQ=34)
                                                                                                                                                                                                                            cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAg                       <  1:662998/219‑1 (MQ=32)
                                                                                                                                                                                                                             gcAGGGCGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGCTGCGACCTCCGAGAACCTGGCGCt                                                                                                                                                 <  2:633223/96‑1 (MQ=255)
                                                                                                                                                                                                                                gggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGgcggcg                                                        <  2:293449/182‑1 (MQ=33)
                                                                                                                                                                                                                                gggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCaacaac   <  1:761624/235‑1 (MQ=33)
                                                                                                                                                                                                                                 ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAAcc  <  1:594421/235‑1 (MQ=33)
                                                                                                                                                                                                                                 ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAAcc  <  1:476596/235‑1 (MQ=32)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
GACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACC  >  AE015451/206594‑207053

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