Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 210,319 C→G 13.6% T5275T (ACC→ACG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451210,3190CG13.6% 257.6 / 20.7 86T5275T (ACC→ACGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base G (7/5);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTC  >  AE015451/210102‑210553
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
acacCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGc                                                                                                                                                                                                                          >  1:649047/1‑236 (MQ=255)
        gCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGAt                                                                                                                                                                                                                  <  1:684971/236‑1 (MQ=255)
            cgcAGGGTGATGACGTTTACACCGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                      <  1:150271/228‑1 (MQ=255)
                 ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATcaacaa                                                                                                                                                                                                            <  2:576087/233‑1 (MQ=255)
                 ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACca                                                                                                                                                                                                         >  1:842352/1‑236 (MQ=255)
                       gaCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCCCCGACGAAGTGGTGGCCAAGCTGACTGCGCCCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                   <  1:255523/236‑1 (MQ=255)
                         cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACcaacca                                                                                                                                                                                                                                 >  1:135253/1‑202 (MQ=255)
                         cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCAt                                                                                                                                                                                                 >  1:1116399/1‑236 (MQ=255)
                             tACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGAc                                                                                                                                                                                             <  1:420431/236‑1 (MQ=255)
                             tACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGAc                                                                                                                                                                                             <  1:188571/236‑1 (MQ=255)
                              aCAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGc                                                                                                                                                                                                                          >  1:308963/1‑206 (MQ=255)
                              aCAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGCACGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCAATCAACAACCACAGTGCATTGAc                                                                                                                                                                                             >  1:859223/1‑235 (MQ=255)
                                          gCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAg                                                                                                                                                                                                                               >  2:316890/1‑189 (MQ=255)
                                           ccAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGc                                                                                                                                                                                >  1:857512/1‑235 (MQ=255)
                                               aTCCGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCATCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATc                                                                                                                                                                                                                 <  2:512753/198‑1 (MQ=255)
                                                     cTGGGCATCCCCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                                     <  1:235406/236‑1 (MQ=255)
                                                       gggCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGCTCAACAACCACAGTg                                                                                                                                                                                                    >  2:732683/1‑203 (MQ=255)
                                                            tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACcaacaa                                                                                                                                                                                                                                 >  1:702749/1‑169 (MQ=255)
                                                              cccTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                      <  1:180097/177‑1 (MQ=255)
                                                                cTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCt                                                                                                                                                                                           >  1:1107835/1‑203 (MQ=255)
                                                                 tCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGCCCAACACAGACGGcc                                                                                                                                                                                                                         >  2:61765/1‑172 (MQ=255)
                                                                     cggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCGCAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCg                                                                                                                                                     >  1:67380/1‑236 (MQ=255)
                                                                     cggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCg                                                                                                                                                                   >  1:413045/1‑222 (MQ=255)
                                                                     cggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                                                                      >  1:630626/1‑235 (MQ=255)
                                                                      ggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaa                                                                                                                                                                        <  1:1100662/216‑1 (MQ=255)
                                                                      ggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaa                                                                                                                                                                        >  2:1100661/1‑216 (MQ=255)
                                                                           aCGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                                                                               >  2:577195/1‑236 (MQ=12)
                                                                           aCGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAACACCGTCATCACCGTGCCGGCCaa                                                                                                                                                >  1:625448/1‑235 (MQ=255)
                                                                             gCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACg                                                                                                                                              >  1:553691/1‑235 (MQ=255)
                                                                                 cTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                         <  1:888523/236‑1 (MQ=255)
                                                                                   ggtgacaCCTTCGAGAACCTGGAGCTGGGCGGTGCCGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCaacaac                                                                                                                                                                                                           <  1:800763/162‑1 (MQ=255)
                                                                                   ggtcgcaCCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACg                                                                                                                                                                            <  2:1111745/193‑1 (MQ=255)
                                                                                    gTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                         <  2:553690/233‑1 (MQ=255)
                                                                                       cGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGc                                                                                                                                   <  1:96979/236‑1 (MQ=255)
                                                                                       cGACCGTCGAGAGCCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCg                                                                                                                                                     <  2:1056260/218‑1 (MQ=255)
                                                                                         aCCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCt                                                                                                                                                                                                                        <  2:1047027/149‑1 (MQ=255)
                                                                                          ccTTCGAGAACCTGGAGCTGGGCGGCGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGAcc                                                                                                                                                                                            >  1:619209/1‑176 (MQ=9)
                                                                                           ctccGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                                          <  2:152964/158‑1 (MQ=255)
                                                                                            ttCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGTCCTGCCGATCAACAACCACAGTGCATTGACCTTCa                                                                                                                                                                                        >  1:568907/1‑178 (MQ=255)
                                                                                                  aaCCTGGAGCTGGGCGGTGCGGCTTCGTTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGcc                                                                                                                                                                                                                         <  1:77859/139‑1 (MQ=255)
                                                                                                  aaCCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGt                                                                                                                        >  1:761249/1‑236 (MQ=255)
                                                                                                  aaCCTGGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCActg                                                                                                                         <  1:270075/235‑3 (MQ=255)
                                                                                                       cgAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGG‑‑CTGCCACCGTCACTGcc                                                                                                                 >  1:301002/2‑236 (MQ=11)
                                                                                                        gAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTg                                                                                                                   <  1:122983/235‑1 (MQ=255)
                                                                                                          gCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACggtc                                                                                                                                                                                                                         <  1:411778/131‑3 (MQ=37)
                                                                                                          gCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACggtc                                                                                                                                                                                                                         >  2:411777/1‑129 (MQ=37)
                                                                                                           cTGGGCGGTGCGGCTTCGGTTCAAGTGACCGCCACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAAc                                                                                                                                               <  2:54956/204‑1 (MQ=255)
                                                                                                           cTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGcccc                                                                                                               <  1:316891/236‑1 (MQ=255)
                                                                                                            ggggCGGAGCGGCTTCGGTTCAAGTGACCGACACCACCGGCGACGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCt                                                                                                                                                                                                                        <  1:334342/129‑1 (MQ=255)
                                                                                                              ggCGGTGCGGCGTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                              <  2:991372/234‑1 (MQ=255)
                                                                                                                   cgcggcgTCGGTCCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACaat                                                                                                        <  2:774870/234‑2 (MQ=11)
                                                                                                                      ggCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTc                                                                                                     >  1:841653/1‑235 (MQ=40)
                                                                                                                      ggCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCAATCAACAACca                                                                                                                                                                                                         >  1:1025332/1‑135 (MQ=18)
                                                                                                                      ggCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCg                                                                                                                                        >  2:843128/1‑200 (MQ=11)
                                                                                                                                     gcccGACACCCCCGACGAAGTGGTGGCCAACCTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGt                                                                                                                                                                                                     <  2:280762/122‑1 (MQ=11)
                                                                                                                                      aCCGACACCCTCGACGAAGTGGTGGCCAAACTGACTGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGc                                                                                                                                            <  2:393274/180‑1 (MQ=255)
                                                                                                                                      aCCGACACCACCGACGAAGTGGTGGCCAAGCTGACGGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTg                                                                                                                   <  2:103552/205‑1 (MQ=255)
                                                                                                                                         tacacCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCAACACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                          <  2:126394/146‑1 (MQ=37)
                                                                                                                                           caccacCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCgc                                                                                                                                                                              >  2:846587/1‑140 (MQ=37)
                                                                                                                                           caccacCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCCTCACCGTGCCGGCCAACGGCCACGTGGGCACt                                                                                                                                >  2:104221/1‑187 (MQ=255)
                                                                                                                                           caccacCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTg                                                                               >  1:457993/1‑236 (MQ=255)
                                                                                                                                                    gagacaggtggCCAAACTGACCGCCACCCCGTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGTAAGACCGTCATCACCGTGCCGGCCa                                                                                                                                                 <  1:600086/155‑1 (MQ=11)
                                                                                                                                                       aGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAATAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGcccc                                                                                                               <  1:942254/192‑1 (MQ=255)
                                                                                                                                                           gtggCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGa                                                                         >  2:722902/1‑226 (MQ=255)
                                                                                                                                                           gtggCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATc                                                                 <  2:955157/234‑1 (MQ=255)
                                                                                                                                                                 aaGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCCCCGGACAACGCCCACgc                                                                                                >  2:133618/1‑196 (MQ=255)
                                                                                                                                                                       accgCGACCCCGTCCGTCACCGAAGGCGGCGAGATCACTTACACCAGCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCa                                                                                                                                                                         <  2:198520/115‑1 (MQ=1)
                                                                                                                                                                         tGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCgg                                                   >  1:572579/1‑234 (MQ=255)
                                                                                                                                                                              ccccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                          <  2:315595/110‑1 (MQ=35)
                                                                                                                                                                                ccgccgGTCACCGAAGGCGGCGAGCTCACCTACACCATCACCCTGACCCACAAAGACGGTCTGCCGCTcc                                                                                                                                                                                                                >  2:308927/1‑69 (MQ=11)
                                                                                                                                                                                ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                          <  2:480835/108‑1 (MQ=38)
                                                                                                                                                                                ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                          >  1:480836/1‑108 (MQ=38)
                                                                                                                                                                                ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACg                                                                                                 >  2:301124/1‑181 (MQ=37)
                                                                                                                                                                                ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTc                                                                                               <  2:311362/183‑1 (MQ=37)
                                                                                                                                                                                ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTc                                                                                               >  1:311363/1‑183 (MQ=37)
                                                                                                                                                                                 cGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCg                                                                            >  2:883263/1‑201 (MQ=255)
                                                                                                                                                                                  gacGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTTCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACTGTGCCGCCCAACGGCACTGTGGGCACt                                                                                                                                <  2:695373/146‑1 (MQ=255)
                                                                                                                                                                                    cGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGt                                      >  1:573339/1‑236 (MQ=255)
                                                                                                                                                                                      gTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGc                                                                                                                                                   >  2:1081637/1‑125 (MQ=38)
                                                                                                                                                                                      gTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGt                                      >  2:1010054/1‑234 (MQ=255)
                                                                                                                                                                                      gTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCAACCGt                                                         >  2:606148/1‑215 (MQ=255)
                                                                                                                                                                                        cACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTa                                     >  2:618228/1‑233 (MQ=255)
                                                                                                                                                                                            gAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTg                                                                                                                                                        <  1:183659/114‑1 (MQ=17)
                                                                                                                                                                                                   gccagaTCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGa                                                                   <  2:946892/189‑1 (MQ=255)
                                                                                                                                                                                                      agaTCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                          <  1:536701/86‑1 (MQ=33)
                                                                                                                                                                                                          cACCTACACCATCACCCTGACCAACAAAGCCGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCGCGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCa                                                                                                                                  <  2:48923/122‑1 (MQ=255)
                                                                                                                                                                                                                        cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGCGCGACGTCAAGACCGTCATCACCGTGCCGGCCAACgtcc                                                                                                                                           >  2:68496/1‑96 (MQ=255)
                                                                                                                                                                                                                        cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGa                                                                                                            >  2:70885/1‑130 (MQ=38)
                                                                                                                                                                                                                        cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCg                                                          >  2:367460/1‑180 (MQ=255)
                                                                                                                                                                                                                        cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCg     >  2:1054960/1‑233 (MQ=255)
                                                                                                                                                                                                                        cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTc  >  1:191321/1‑236 (MQ=255)
                                                                                                                                                                                                                        cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCg     >  2:228358/1‑233 (MQ=255)
                                                                                                                                                                                                                        cccTGACCAACAAAGACGGCCTGCCGATCAACAACAACAGCGCATTGACCTTCACGCTGAGCGACGGCCAGACCGTCATCACCGTGCCGGCCCACGGCACCGTGGGCa                                                                                                                                  >  2:727255/1‑108 (MQ=11)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTC  >  AE015451/210102‑210553

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