Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 206,311 G→C 5.2% V3939V (GTG→GTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451206,3110GC 5.2% 290.3 / 6.1 94V3939V (GTG→GTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (38/50);  new base C (2/3);  total (41/53)
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.65e-01

GCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGT  >  AE015451/206079‑206538
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   
gCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGAc                                                                                                                                                                                                                                   >  1:1127911/1‑235 (MQ=255)
           cAGGTCACGATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCTCAAGGCGGTGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCATCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggcta                                                                                                                                                                                                                       <  1:794098/236‑3 (MQ=255)
                 aCCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCAc                                                                                                                                                                                                                 <  1:1254189/236‑1 (MQ=255)
                           ccGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGCCGGTGACCGAAGGCGGCGAGATCACCTACACCAt                                                                                                                                                                                                        >  1:377599/1‑235 (MQ=255)
                               cGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACg                                                                                                                                                                                                    <  1:1028890/235‑1 (MQ=255)
                                gagaCCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGgcggcg                                                                                                                                                                                                                        <  2:471116/214‑1 (MQ=255)
                                      aGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGAc                                                                                                                                                                                               >  2:1187687/1‑233 (MQ=255)
                                           gcCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcg                                                                                                                                                                                                                           >  2:643408/1‑200 (MQ=255)
                                              gtACACCCACGCTGCGCAGGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACg                                                                                                                                                                                    <  1:977084/235‑1 (MQ=255)
                                                  acCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACg                                                                                                                                                                                                    >  2:996838/1‑216 (MQ=255)
                                                    ccACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCCGc                                                                                                                                                                              <  1:632103/236‑1 (MQ=255)
                                                      aCGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCg                                                                                                                                                                            <  1:625886/236‑1 (MQ=255)
                                                         cTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATc                                                                                                                                                                         >  1:580545/1‑236 (MQ=255)
                                                          tGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcg                                                                                                                                                                                                                           <  1:861546/185‑1 (MQ=255)
                                                          tGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcg                                                                                                                                                                                                                           >  2:861541/1‑185 (MQ=255)
                                                          tGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                       <  1:30177/205‑1 (MQ=255)
                                                          tGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATca                                                                                                                                                                        >  1:669156/1‑236 (MQ=255)
                                                            cgcAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGTGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGACGGCGGCGCGATCACCTACACCATCACGCTGACCAACACAgcc                                                                                                                                                                                     >  1:427746/1‑219 (MQ=255)
                                                              cAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                              >  2:1228192/1‑226 (MQ=255)
                                                               aaGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGCTGTAGACGGTCGCACCCTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcgg                                                                                                                                                                                                                          >  1:1177475/1‑181 (MQ=255)
                                                                aGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTc                                                                                                                                                                                 <  1:16823/221‑1 (MQ=255)
                                                                     aTGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGgcggcg                                                                                                                                                                                                                        <  2:240123/177‑1 (MQ=255)
                                                                       gACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATc                                                                                                                                                                         <  1:354052/222‑1 (MQ=255)
                                                                          gTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAACTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCAt                                                                                                                                                         <  1:699789/235‑1 (MQ=255)
                                                                                       cTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTg                                                                                                                                           >  1:87685/1‑236 (MQ=255)
                                                                                           ccAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                              >  1:45231/1‑197 (MQ=255)
                                                                                                tCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                  <  1:858429/236‑1 (MQ=255)
                                                                                                  aGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGCTCACCTACACCATCACGCTGCCcaa                                                                                                                                                                                            >  2:284466/1‑176 (MQ=255)
                                                                                                  aGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGAcc                                                                                                                                                    >  2:747933/1‑216 (MQ=255)
                                                                                                     cTGTGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                             <  1:31780/236‑1 (MQ=255)
                                                                                                       gggtaTCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCg                                                                                                                           >  1:809278/5‑236 (MQ=255)
                                                                                                       gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAg                                                                                                                                                                                                                              <  1:739557/137‑1 (MQ=255)
                                                                                                       gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCga                                                                                                                                                                                                                       >  2:583112/1‑144 (MQ=255)
                                                                                                       gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACcaaca                                                                                                                                                                                          >  2:757100/1‑173 (MQ=255)
                                                                                                       gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                             >  1:956163/1‑234 (MQ=255)
                                                                                                       gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGCCGGTGACCGAAGGCGGCGAGATCAccg                                                                                                                                                                                                               >  2:176599/1‑151 (MQ=255)
                                                                                                                gTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACACCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACACAGACGGTCTGCCGATCAACAACCACAg                                                                                                                                                              >  1:97149/1‑192 (MQ=255)
                                                                                                                     ctgTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAACGACGGTCTGCCAATCAACAACCACAGTGCAt                                                                                                                                                         >  2:883754/1‑192 (MQ=255)
                                                                                                                             gTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACcaaca                                                                                                                                                                                          >  1:540197/1‑151 (MQ=255)
                                                                                                                             gTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                                       <  2:1127906/234‑1 (MQ=255)
                                                                                                                              tAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGt                                                                                                    <  1:1228197/236‑1 (MQ=40)
                                                                                                                                 aCGGTCGGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                         <  1:1001032/116‑1 (MQ=255)
                                                                                                                                 aCGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                              >  1:1088556/1‑223 (MQ=40)
                                                                                                                                 aCGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                              <  2:1088551/223‑1 (MQ=39)
                                                                                                                                  cGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGt                                                                                                                                                             >  2:751886/1‑175 (MQ=255)
                                                                                                                                    gtggCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCgag                                                                                                                                                                                                                      <  1:1240882/113‑1 (MQ=38)
                                                                                                                                     tgcgaCCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                       <  2:1070845/126‑1 (MQ=25)
                                                                                                                                     tgcgaCCATCGAGAACCTGCAACTGGGCGGCAACACTTCGGCTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCATCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGt                                                                                               <  2:720951/230‑1 (MQ=9)
                                                                                                                                     gcGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACAGCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCCTTGACCTTCACGCGGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                              <  2:739639/218‑1 (MQ=255)
                                                                                                                                       gCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTgg                                                                                             >  2:1006792/1‑234 (MQ=16)
                                                                                                                                         aCCTGCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTAc                                                                                                  <  1:1072364/227‑1 (MQ=25)
                                                                                                                                         aCCGTCGAGAACCTGCAGCTGGGCGGCAACGCTTCGGTGCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGAt                                                                                                                                                                          <  2:427741/155‑1 (MQ=255)
                                                                                                                                          ccTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                         <  2:1174238/107‑1 (MQ=38)
                                                                                                                                          ccTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCg                                                                                                                                                                            <  2:753053/152‑1 (MQ=38)
                                                                                                                                            ctCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGc                                                                                      <  1:132024/235‑1 (MQ=37)
                                                                                                                                                 gaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGtcc                                                                                                  >  1:509134/1‑217 (MQ=21)
                                                                                                                                                  aaCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                  <  1:290261/186‑1 (MQ=35)
                                                                                                                                                   aCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCa                                                                                                         <  2:540192/210‑1 (MQ=34)
                                                                                                                                                          ccTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCGGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                              <  2:509129/197‑1 (MQ=21)
                                                                                                                                                          aCTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTggg                                                                                            >  1:607339/1‑216 (MQ=35)
                                                                                                                                                               gcggcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                                                                                                <  1:1223278/191‑1 (MQ=34)
                                                                                                                                                                 ggcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                              <  2:581058/191‑1 (MQ=34)
                                                                                                                                                                 ggcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAAcc                                                                                                                                                                  >  1:727725/1‑139 (MQ=37)
                                                                                                                                                                  gcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                  <  2:1198112/138‑1 (MQ=17)
                                                                                                                                                                  gcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGGCGGTGACCGAAGGCGGCGAGATCAcc                                                                                                                                                                                                                >  2:670748/1‑92 (MQ=255)
                                                                                                                                                                      cGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCgg                                                                     <  1:1091970/227‑1 (MQ=39)
                                                                                                                                                                         ttCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGcccc                                                                       <  1:121089/222‑1 (MQ=39)
                                                                                                                                                                           cGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                  <  1:511634/161‑1 (MQ=14)
                                                                                                                                                                             gTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                              <  2:45230/115‑1 (MQ=39)
                                                                                                                                                                              ttCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGCCCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGt                                                                                                    <  1:836503/188‑1 (MQ=25)
                                                                                                                                                                                 aaGTCACCGACACCACCGACGAAGTCGTGGCCAAGCTGACCGCCACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGa                                                                                                                                                      <  2:361150/135‑1 (MQ=18)
                                                                                                                                                                                           ccccaccGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGAACAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCg                                                                                                             <  2:760096/164‑1 (MQ=21)
                                                                                                                                                                                           caccacCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggc                                                                                                                                                                                                                            >  2:451628/1‑55 (MQ=35)
                                                                                                                                                                                           caccacCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGAcc                                          <  2:669151/233‑1 (MQ=40)
                                                                                                                                                                                                   acgaAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGAc                                                                                                                                                                                     <  1:1249418/86‑1 (MQ=38)
                                                                                                                                                                                                   acgaAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGAt                                <  1:751891/235‑1 (MQ=40)
                                                                                                                                                                                                        gtggtgGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                              <  2:957260/88‑1 (MQ=38)
                                                                                                                                                                                                        gtggtgGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGa                           <  2:1193147/235‑1 (MQ=25)
                                                                                                                                                                                                             ggCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                         <  2:858442/40‑1 (MQ=255)
                                                                                                                                                                                                              gCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGGGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGTTACCGTGGGTACTGCTACCGtca                                                                              >  2:157062/1‑176 (MQ=255)
                                                                                                                                                                                                                 aagttGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCg                                                                      >  1:660303/5‑183 (MQ=255)
                                                                                                                                                                                                                 aaGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCg                                                                      <  2:660298/183‑1 (MQ=255)
                                                                                                                                                                                                                 aaGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGt                 <  1:1187692/236‑1 (MQ=40)
                                                                                                                                                                                                                 aaGCTGACCGCGACCCCGTCGGGGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGggggg                                                                                           >  2:318568/1‑161 (MQ=25)
                                                                                                                                                                                                                     tGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCGACgcc                                                       >  1:241654/1‑192 (MQ=21)
                                                                                                                                                                                                                      gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                  <  1:685485/86‑1 (MQ=38)
                                                                                                                                                                                                                      gACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCg                                                                      <  1:863812/178‑1 (MQ=14)
                                                                                                                                                                                                                         cgcgACCCCGTCGGTGACCGACGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCAc                                                  <  2:205959/195‑1 (MQ=39)
                                                                                                                                                                                                                         cgcgACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACAGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAgg           >  2:693894/1‑234 (MQ=9)
                                                                                                                                                                                                                           cgACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCg       >  1:1077322/1‑236 (MQ=39)
                                                                                                                                                                                                                            gACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCg       >  1:1186313/1‑235 (MQ=38)
                                                                                                                                                                                                                                ccttCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                                                                                                   >  2:1260890/4‑75 (MQ=255)
                                                                                                                                                                                                                                ccGTCGGTGACGGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGc                                                                          >  1:1206259/1‑164 (MQ=39)
                                                                                                                                                                                                                                ccGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGt  >  1:608170/1‑236 (MQ=38)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   
GCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGT  >  AE015451/206079‑206538

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