Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 202,064 G→A 7.8% V2524I (GTC→ATC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451202,0640GA 7.8% 198.4 / 5.4 64V2524I (GTC→ATC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (2/3);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCG  >  AE015451/201861‑202297
                                                                                                                                                                                                           |                                                                                                                                                                                                                                         
tCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAcc                                                                                                                                                                                                           <  1:407904/236‑1 (MQ=40)
    cggcgCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCa                                                                                                                                                                                                       <  1:153130/236‑1 (MQ=40)
     ggcgCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCTGGTCACCATCa                                                                                                                                                                                                       <  1:757232/235‑1 (MQ=255)
     ggcgCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCa                                                                                                                                                                                                       <  1:482488/235‑1 (MQ=255)
     ggcgCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGCCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCa                                                                                                                                                                                                       <  1:215140/235‑1 (MQ=255)
      gcgcAGATGCGTGGCAGTTCGAAACCCTGAACCTGGACAAGACACCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGca                                                                                                                                                                                                                            <  1:87951/213‑1 (MQ=255)
             tGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGCGCAACAACGCGCAggc                                                                                                                                                                                                               >  2:81394/1‑218 (MQ=255)
                       ttCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGt                                                                                                                                                                                      <  2:758113/234‑1 (MQ=255)
                         cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGccaca                                                                                                                                                                                                                         >  1:189120/1‑197 (MQ=255)
                         cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCa                                                                                                                                                                                        <  2:698734/230‑1 (MQ=255)
                            aaaCCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGAcc                                                                                                                                                                                         >  2:308264/1‑226 (MQ=255)
                            aaaCCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                                                                                                                 >  2:177088/1‑202 (MQ=255)
                                       tgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCacgac                                                                                                                                                                    <  1:698130/236‑1 (MQ=255)
                                        ggACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCacgac                                                                                                                                                                    <  1:782691/235‑1 (MQ=255)
                                           cAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGAcgc                                                                                                                                                                  <  2:422469/234‑1 (MQ=255)
                                                ccccGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCGGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAgg                                                                                                                                                            <  1:527305/235‑1 (MQ=255)
                                                       cAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGGAACGAAGGCGACATCGGCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGaa                                                                                                                                                                                                                                    >  2:665533/1‑155 (MQ=255)
                                                       aaGCACTACCGTCACTGACGAGCCAGGCAATCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTgcg                                                                                                                                                                                   <  2:397568/204‑1 (MQ=255)
                                                          cACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCa                                                                                                                                                                                        <  2:587121/197‑1 (MQ=255)
                                                           aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAg                                                                                                                                                                                                                              <  1:443163/158‑1 (MQ=255)
                                                           aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCacg                                                                                                                                                                      >  2:504146/1‑214 (MQ=255)
                                                           aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGAt                                                                                                                                                        >  1:84531/1‑228 (MQ=255)
                                                           aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAg                                                                                                                                                                                                                              >  2:443161/1‑158 (MQ=255)
                                                            cTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTgaga                                                                                                                                                                                           <  1:769127/192‑1 (MQ=255)
                                                                     cTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCt                                                                                                                                      <  1:177090/236‑1 (MQ=40)
                                                                        aCGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCgg                                                                                                                                                                                                <  2:178614/175‑1 (MQ=39)
                                                                                 gCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGc                                                                                                                          <  1:781246/236‑1 (MQ=255)
                                                                                   aCTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCg                                                                                                                                                          >  2:647435/1‑202 (MQ=40)
                                                                                             aCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCAc                                                                                                               <  1:473387/235‑1 (MQ=255)
                                                                                                         tCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACAATCGAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACgg                                                                                                                                                   >  1:458352/1‑186 (MQ=255)
                                                                                                         ccGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTa                                                                                                                                                <  1:647210/189‑1 (MQ=255)
                                                                                                          cGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                                                                                                         >  1:113067/1‑100 (MQ=255)
                                                                                                           gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGcc                                                                                                                                  <  2:717634/202‑1 (MQ=255)
                                                                                                           gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATg                                                                                                 <  1:81395/235‑1 (MQ=255)
                                                                                                            tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGAt                                                                                                  >  1:288220/1‑233 (MQ=255)
                                                                                                                        tCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTg                                                                                                                                                                                                                                      <  2:506511/89‑1 (MQ=255)
                                                                                                                        tCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTg                                                                                                                                                                                                                                      >  1:506513/1‑89 (MQ=255)
                                                                                                                         cACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCt                                                                                  >  1:428944/1‑236 (MQ=255)
                                                                                                                         cACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCt                                                                                  >  1:672941/1‑236 (MQ=255)
                                                                                                                              ccGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACCACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGCACACCCACGACgcgc                                                                                                                                                                >  1:121738/1‑153 (MQ=255)
                                                                                                                                 aCCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                          <  1:318993/204‑1 (MQ=255)
                                                                                                                                 aCCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                          >  2:318991/1‑204 (MQ=255)
                                                                                                                                 aCCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGt                                                                                         >  2:695723/1‑221 (MQ=255)
                                                                                                                                         tCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGCCCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAAc                                                                   >  1:415345/1‑235 (MQ=255)
                                                                                                                                           ggTAGCCGAGAACGTCAAACCAACCTTCACTGTGCACGTCAACCAGCCGCTGGCTCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATg                                                                                                 <  2:288218/203‑1 (MQ=255)
                                                                                                                                              aGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTggg                                                               >  1:597570/1‑234 (MQ=255)
                                                                                                                                                ccGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAgg                                                                                                                                           <  2:165097/156‑1 (MQ=255)
                                                                                                                                                ccGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAgg                                                                                                                                           >  1:165099/1‑156 (MQ=255)
                                                                                                                                                        gTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGcccgcg                                                                                                                     >  2:166962/1‑166 (MQ=255)
                                                                                                                                                            aGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCt                                                                                                                                      <  2:411054/149‑1 (MQ=255)
                                                                                                                                                                aaCCTTCACCGTGCACGTCAACCAGCCGCTGGCGCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAg                                                                                                                                                             <  2:263/122‑1 (MQ=255)
                                                                                                                                                                  ccTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTcgc                                                                                                                                                                                                            <  2:116733/73‑3 (MQ=21)
                                                                                                                                                                           tGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACg                                                       <  2:676446/213‑1 (MQ=255)
                                                                                                                                                                            gaaCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                           <  2:113595/158‑1 (MQ=255)
                                                                                                                                                                                gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATg                                                                                                 >  2:598261/1‑166 (MQ=255)
                                                                                                                                                                                gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATg                                                                                                 <  1:598263/166‑1 (MQ=255)
                                                                                                                                                                                gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTg                                                                       >  1:600596/1‑192 (MQ=255)
                                                                                                                                                                                gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTcga                           >  1:131361/1‑236 (MQ=255)
                                                                                                                                                                                        gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGt                                                                                         >  2:205074/1‑166 (MQ=255)
                                                                                                                                                                                        gCCGCTGGCTCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCa                                                                                                                            <  1:704211/131‑1 (MQ=38)
                                                                                                                                                                                        gCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCa                                                                                                                            >  2:704209/1‑131 (MQ=17)
                                                                                                                                                                                             tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGc               >  1:119194/1‑235 (MQ=255)
                                                                                                                                                                                                     gACCTGGGCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGAcgc                                                                                                                                                                  <  2:119991/80‑1 (MQ=255)
                                                                                                                                                                                                         tGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACcgcg  <  1:308266/236‑1 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                         
TCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCG  >  AE015451/201861‑202297

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: