Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 206,311 G→C 6.6% 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 6.6% 278.7 / 6.6 89V3939V (GTG→GTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (37/45);  new base C (2/4);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCA  >  AE015451/206079‑206543
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
gCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGccc                                                                                                                                                                                                                                       >  1:129149/1‑236 (MQ=255)
          ccAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                              <  1:992901/235‑1 (MQ=255)
          acAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGTCACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                              <  1:940784/234‑1 (MQ=255)
            aGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCTCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAgg                                                                                                                                                                                                                                  <  2:948887/229‑1 (MQ=11)
                             ggCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCAc                                                                                                                                                                                                          <  1:180925/236‑1 (MQ=255)
                                            ccATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                              <  2:681626/201‑1 (MQ=255)
                                            ccATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGAt                                                                                                                                                                                                                         <  2:64601/206‑1 (MQ=255)
                                                 cacCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTc                                                                                                                                                                                      <  1:889553/236‑1 (MQ=255)
                                                  acCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTc                                                                                                                                                                                      <  1:902512/235‑1 (MQ=255)
                                                            cgcAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCt                                                                                                                                                                                     <  2:929581/226‑1 (MQ=255)
                                                                aGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                       <  1:920753/236‑1 (MQ=255)
                                                                  gCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                                                                                                        <  2:921284/233‑1 (MQ=255)
                                                                  gCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                       <  1:1030000/234‑1 (MQ=255)
                                                                         cGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACAcc                                                                                                                                                                                                               >  2:1075613/1‑187 (MQ=255)
                                                                         cGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTacac                                                                                                                                                                                                                >  1:452446/1‑186 (MQ=255)
                                                                               cAACGACACTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGTCAACGCTTCGCTTCAAGTGACCGACACCAGCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCt                                                                                                                                                        <  1:1028637/236‑1 (MQ=255)
                                                                                aaCGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAccc                                                                                                                                                                      >  2:580477/1‑220 (MQ=255)
                                                                                      gCTGGCCAGAGCAGCCTGCGCATCACGTCTGCTGTGGATGTAGCCGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTGCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                            <  2:290752/177‑1 (MQ=255)
                                                                                       cTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTg                                                                                                                                                >  1:821588/1‑236 (MQ=255)
                                                                                        tGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTAACCGAAGGCGGCGAGATCACCTacac                                                                                                                                                                                                                >  1:749233/1‑171 (MQ=255)
                                                                                                tCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                       <  1:1075616/236‑1 (MQ=255)
                                                                                                       gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGAc                                                                                                                                          >  1:406363/1‑226 (MQ=255)
                                                                                                       gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaaaa                                                                                                                                   >  1:227141/1‑233 (MQ=255)
                                                                                                         gCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACg                                                                                                                                                                                         <  2:625573/177‑1 (MQ=255)
                                                                                                         gCATCACGGCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATc                                                                                                                                                                                                            <  2:757731/158‑1 (MQ=255)
                                                                                                            tCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                       <  2:465132/192‑1 (MQ=255)
                                                                                                            tCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGt                                                                                                                                                                  >  2:861290/1‑197 (MQ=255)
                                                                                                                gTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACAACAGTGCATTGACCTTCACGCTGAgag                                                                                                                                            >  1:76252/1‑213 (MQ=255)
                                                                                                                gTCTGCGGTGGATGTAGCCGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCCTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                            <  1:147086/151‑1 (MQ=255)
                                                                                                                 tCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCgc                                                                                                                      >  1:328127/1‑235 (MQ=255)
                                                                                                                  ctgctgTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTc                                                                                                                                                      >  2:908465/1‑203 (MQ=255)
                                                                                                                                  gggtgcgaCCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCt                                                                                                                                                                                     <  2:1035260/149‑1 (MQ=17)
                                                                                                                                  cGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGCGATCACCTACACCATCACGCTGACcaacaa                                                                                                                                                                                              >  2:283065/1‑147 (MQ=255)
                                                                                                                                  cGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATcaacaa                                                                                                                                                                         >  2:700240/1‑168 (MQ=40)
                                                                                                                                   ggTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGcca                                                                                                                                                                                 >  1:93291/1‑158 (MQ=40)
                                                                                                                                   ggTCGCACCTTCGAAAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGt                                                                                                    <  1:499902/236‑1 (MQ=11)
                                                                                                                                    gtgcgaCCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                                                                                                                          <  2:747705/176‑1 (MQ=17)
                                                                                                                                    gTCCCACCTTCGAAAACCTGCAACTGGGCGGCAACGCTGCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCCCGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGt                                                                                                    <  1:269867/235‑1 (MQ=255)
                                                                                                                                     tcccACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACcaa                                                                                                                                                                                                 <  1:705797/138‑1 (MQ=38)
                                                                                                                                       gCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAACg                                                                                                           >  2:869180/1‑225 (MQ=25)
                                                                                                                                        cACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAg                                                                                                                                                                   >  1:557602/1‑168 (MQ=37)
                                                                                                                                        cACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTgg                                                                                                  >  2:968362/1‑233 (MQ=37)
                                                                                                                                        cACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGACTg                                                                                                                                                                                    >  2:1089767/1‑151 (MQ=37)
                                                                                                                                        cACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTCGTGGCCAAGCTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATcaac                                                                                                                                                                           >  1:454574/1‑160 (MQ=2)
                                                                                                                                        cACCTCCGAGAACCTGCAACGGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCa                                                                                                              <  1:306450/221‑1 (MQ=21)
                                                                                                                                           cTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTaca                                                                                                                                                                                                                 >  2:653232/1‑119 (MQ=40)
                                                                                                                                           cTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                                                                                                        <  2:879192/160‑1 (MQ=35)
                                                                                                                                            ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGt                                                                                                    <  2:676809/227‑1 (MQ=18)
                                                                                                                                            ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGc                                                                                           <  1:631020/236‑1 (MQ=39)
                                                                                                                                                 gaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                              <  1:421114/100‑1 (MQ=37)
                                                                                                                                                 gaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                                              >  2:421111/1‑100 (MQ=37)
                                                                                                                                                  aaCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCg                                                                                      <  1:319367/235‑1 (MQ=37)
                                                                                                                                                        cAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGt                                                                                                                               >  2:425032/1‑188 (MQ=35)
                                                                                                                                                         aaCTGGGCGGCACCGCGTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTgg                                                                                                  <  2:325814/216‑1 (MQ=25)
                                                                                                                                                         aaCTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCg                                                                                                                       >  2:769610/1‑195 (MQ=14)
                                                                                                                                                         aaCTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCg                                                                                                                       <  1:769613/195‑1 (MQ=34)
                                                                                                                                                          aCTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTAc                                                                                              >  2:799014/1‑219 (MQ=16)
                                                                                                                                                           cTGGGCGGCAACGCTTCGGTTGAAGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                                   <  2:766208/133‑1 (MQ=255)
                                                                                                                                                            tGGGCGGCAACGCTTCGGTTCAAGTCGCCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                                                                                                     <  2:133161/194‑1 (MQ=21)
                                                                                                                                                              ggcggcGACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                       <  1:930846/174‑1 (MQ=25)
                                                                                                                                                                 ggcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTACACGCTGAGCGACGGCAAAAcc                                                                                                                                 >  1:6471/1‑177 (MQ=25)
                                                                                                                                                                  gcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                            >  1:23047/1‑101 (MQ=40)
                                                                                                                                                                  gcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                            <  2:23046/101‑1 (MQ=40)
                                                                                                                                                                      cgcgTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACca                                                                                                                                                                      >  1:766211/1‑135 (MQ=21)
                                                                                                                                                                      cGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGAc                                                                        >  2:687915/1‑229 (MQ=39)
                                                                                                                                                                            ggTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGTTGACCGCAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCGCAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAcc                                                                                                                        <  1:536968/175‑1 (MQ=255)
                                                                                                                                                                             gTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGCCCTTCACGCTGAGCGACGGCAAAAcc                                                                                                                                 >  2:280960/1‑165 (MQ=25)
                                                                                                                                                                              ttCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGc                                                         <  1:547700/236‑1 (MQ=39)
                                                                                                                                                                                 aactCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCAc                                                                                                                                                                                                                      <  1:194585/73‑1 (MQ=255)
                                                                                                                                                                                  aGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCa                                                     <  1:506537/236‑1 (MQ=255)
                                                                                                                                                                                   gTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCaa                                                    <  1:663068/236‑1 (MQ=255)
                                                                                                                                                                                     cACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGt                                                                                                    >  2:654123/1‑186 (MQ=35)
                                                                                                                                                                                        cGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGAcc                                               <  1:620845/236‑1 (MQ=255)
                                                                                                                                                                                           ccccaccGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCCCCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGCCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGt                                                             <  1:211297/217‑1 (MQ=255)
                                                                                                                                                                                           caccacCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGt                                                                                                <  1:400990/184‑1 (MQ=17)
                                                                                                                                                                                                 cgacgaAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCg                                                                                      >  1:805366/1‑188 (MQ=39)
                                                                                                                                                                                                 cgacgaAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCg                                                                                      <  2:805363/188‑1 (MQ=39)
                                                                                                                                                                                                   acgaAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGcta                                                                                                              >  2:655966/1‑160 (MQ=18)
                                                                                                                                                                                                    cgaAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATc                                                                                                                                                                                                                        >  2:194583/1‑55 (MQ=35)
                                                                                                                                                                                                       aGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGa                                <  1:908468/236‑1 (MQ=40)
                                                                                                                                                                                                           gtggCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGc                                                                                                                                                                                                        <  1:462219/64‑1 (MQ=40)
                                                                                                                                                                                                           gtggCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCg                             >  1:207690/1‑235 (MQ=38)
                                                                                                                                                                                                            tggCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCAc                                                                                                                                                                                                                      <  1:1004400/49‑1 (MQ=40)
                                                                                                                                                                                                            tggCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCAGCGAGATCAc                                                                                                                                                                                                                      >  2:1004397/1‑49 (MQ=255)
                                                                                                                                                                                                                      gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACg                                                              >  1:584223/1‑191 (MQ=38)
                                                                                                                                                                                                                      gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCg                                         >  2:354869/1‑212 (MQ=39)
                                                                                                                                                                                                                                ccGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTAc                                                               >  1:202686/1‑180 (MQ=39)
                                                                                                                                                                                                                                ccGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGtt                     >  1:311298/1‑222 (MQ=37)
                                                                                                                                                                                                                                     ggTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCa  <  1:377811/236‑1 (MQ=40)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
GCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCA  >  AE015451/206079‑206543

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