Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 201,859 G→C 7.1% A2455A (GCG→GCC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451201,8590GC 7.1% 249.5 / 4.9 84A2455A (GCG→GCCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (43/34);  new base C (2/4);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAC  >  AE015451/201634‑202084
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
gACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcgg                                                                                                                                                                                                                           <  2:675248/234‑1 (MQ=40)
                          gcggcgAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGc                                                                                                                                                                                                                  >  1:774467/1‑217 (MQ=32)
                           cggcgAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacct                                                                                                                                                                                                >  2:862994/1‑234 (MQ=33)
                            tgcgAGATCACCTGCACCATCACGCTGACCAACAAAGACGGTCGGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctg                                                                                                                                                                                               <  2:473459/233‑1 (MQ=21)
                                        tACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAg                                                                                                                                                                                   <  2:594024/234‑1 (MQ=34)
                                                  cGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAg                                                                                                                                                                        >  1:553206/1‑235 (MQ=14)
                                                   gCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACAGTCGTCGTGGCAGCCAATAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATTCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGc                                                                                                                                                                       <  1:572488/235‑1 (MQ=21)
                                                          aacaaAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGGCAGCGGCGCAg                                                                                                                                                                                                                      >  1:54172/1‑181 (MQ=18)
                                                                 aCGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGAc                                                                                                                                                        <  1:501579/236‑1 (MQ=255)
                                                                    gTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAg                                                                                                                                                     >  1:76464/1‑236 (MQ=255)
                                                                    gTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                      >  2:639114/1‑219 (MQ=25)
                                                                        gCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGtt                                                                                                                                                                                                         <  1:932799/180‑1 (MQ=35)
                                                                               aaTAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTGGACAAGACCCCGGTCAGc                                                                                                                                                                       <  2:819778/207‑1 (MQ=21)
                                                                                          ggAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcc                                                                                                                                                                                <  1:255999/187‑1 (MQ=35)
                                                                                          ggAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcc                                                                                                                                                                                >  2:255996/1‑187 (MQ=35)
                                                                                          ggAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGTTCAGCACTACCGTCACTGACTAg                                                                                                                                                     >  1:233472/1‑214 (MQ=255)
                                                                                          ggAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGCCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAg                                                                                                                               >  1:671993/1‑236 (MQ=255)
                                                                                           gAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTggc                                                                                                                                                                                                             >  1:299843/1‑156 (MQ=25)
                                                                                              cTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcg                                                                                                                                                                                >  2:221080/1‑182 (MQ=25)
                                                                                                  aCTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGa                                                                                                                                                         >  2:872832/1‑202 (MQ=255)
                                                                                                  aCTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATc                                                                                                                       <  1:872115/236‑1 (MQ=255)
                                                                                                   cTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATc                                                                                                                       <  1:639772/235‑1 (MQ=255)
                                                                                                         gCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTa                                                                                                                                                                   >  2:612265/1‑185 (MQ=37)
                                                                                                         gCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTc                                                                                                                    >  1:449487/1‑232 (MQ=255)
                                                                                                         gCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAg                                                                                                                 >  1:832335/1‑235 (MQ=255)
                                                                                                         gCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAg                                                                                                                 >  1:576835/1‑235 (MQ=255)
                                                                                                         gCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                      >  2:921269/1‑182 (MQ=25)
                                                                                                            gACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTc                                                                                                                    >  1:558343/1‑229 (MQ=255)
                                                                                                            gACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGt                                                                                                               >  1:30266/1‑234 (MQ=255)
                                                                                                              ccGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCCGGCAACGAAGGCGACATCGTCAAGGGCAc                                                                                                            >  1:84538/1‑235 (MQ=255)
                                                                                                                gATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCg                                                                                                                      >  1:999695/1‑223 (MQ=255)
                                                                                                                gATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCAATCCAGGCAACGAAGGCGACATCg                                                                                                                      <  2:999692/223‑1 (MQ=255)
                                                                                                                    gccccactGTCGTCGTGGCAGCCAACAGCACCACCGGCTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGcc                                                                                                                                                   <  2:726548/185‑1 (MQ=255)
                                                                                                                               gtcgtGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTc                                                                                                                    <  2:30265/210‑1 (MQ=255)
                                                                                                                                tcgtGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACa                                                                                                                                                                                     >  2:439829/1‑144 (MQ=37)
                                                                                                                                tcgtGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACt                                                                                         >  1:999820/1‑236 (MQ=255)
                                                                                                                                     gCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                      >  1:846449/1‑154 (MQ=38)
                                                                                                                                     gCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTa                                                                                                                                                                   >  2:857207/1‑157 (MQ=38)
                                                                                                                                          cAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGcc                                                                                >  1:1036345/1‑235 (MQ=255)
                                                                                                                                            aCAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGt                                                                                                                                                                                                                                 >  2:1081049/1‑88 (MQ=40)
                                                                                                                                            aCAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGt                                                                                                                                                                                                                                 <  1:1081052/88‑1 (MQ=40)
                                                                                                                                            aCAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCGTCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCg                                                                                                                            >  2:406479/1‑189 (MQ=255)
                                                                                                                                            aCAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGATCAATGCCGTCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCg                                                                                                                            <  1:406482/189‑1 (MQ=255)
                                                                                                                                                  ccaccGGCTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTc                                                                       <  1:72650/236‑1 (MQ=25)
                                                                                                                                                   caccGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAgg                                                                                                                              >  2:443025/1‑180 (MQ=255)
                                                                                                                                                   caccGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCa                                                                      >  1:37252/1‑236 (MQ=255)
                                                                                                                                                      cGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                      >  1:595490/1‑137 (MQ=39)
                                                                                                                                                       ggTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAAGTCACCATCACTGCCGACCAGACCTCGGTGGCCGAGAACGTCaaacc                                                                  <  1:862997/236‑4 (MQ=255)
                                                                                                                                                        gTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctgaatct                                                                                                                                                                                          <  1:864184/115‑4 (MQ=0)
                                                                                                                                                          tCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGa                                                                                                                                 <  1:786389/170‑1 (MQ=255)
                                                                                                                                                           cGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGAc                                                                                                                                                        >  1:375657/1‑146 (MQ=255)
                                                                                                                                                           cGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGAc                                                                                                                                                        <  2:375654/146‑1 (MQ=255)
                                                                                                                                                              ccACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                      >  1:957397/1‑129 (MQ=25)
                                                                                                                                                                aCTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCAc                                                                                                                                                                     >  2:147391/1‑128 (MQ=38)
                                                                                                                                                                     aGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACGACCGTCa                                                                                                                                                             >  2:265386/1‑131 (MQ=255)
                                                                                                                                                                      gccgcaCCGGACAACGTCGATGTCGGCACCAACGCGCCGGTGGTCAAAGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTcc                                                                                                                                          <  1:245660/149‑1 (MQ=255)
                                                                                                                                                                               gACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCCTCAcc                                                        >  2:550797/1‑222 (MQ=255)
                                                                                                                                                                                 cAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACGTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGc                                        >  1:1028207/1‑236 (MQ=255)
                                                                                                                                                                                 cAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGc                                        >  1:759866/1‑236 (MQ=255)
                                                                                                                                                                                  agcGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCt                                                             <  1:138192/212‑1 (MQ=39)
                                                                                                                                                                                   aCGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGcc                                       <  1:413726/235‑1 (MQ=255)
                                                                                                                                                                                   aCGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGcc                                       <  1:691563/235‑1 (MQ=255)
                                                                                                                                                                                       cTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAgg                                                                                                                              <  2:430526/144‑1 (MQ=37)
                                                                                                                                                                                       cTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCACCGTCACCGACGAGCCAGGCACTCCAGGCAACGAAgg                                                                                                                              >  1:430529/1‑144 (MQ=255)
                                                                                                                                                                                        tATGTCGGCACCAACGCACCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCCCGGCCGACCAGACTTCGGTAGCCga                                                                              >  2:1083087/1‑191 (MQ=255)
                                                                                                                                                                                         aTGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGTCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCTGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACGTCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGc                                 <  1:1026553/235‑1 (MQ=255)
                                                                                                                                                                                          tGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACGGTGACTGACGAGCCAGGCACTCCAGGCAACGAAGGCg                                                                                                                            >  1:1042750/1‑143 (MQ=255)
                                                                                                                                                                                          tGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCg                                                                                                                            <  2:1042747/143‑1 (MQ=38)
                                                                                                                                                                                                 aCCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCAGCACggg                                                                                                   >  2:413723/1‑160 (MQ=25)
                                                                                                                                                                                                   cAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTgg                      >  1:550153/1‑236 (MQ=255)
                                                                                                                                                                                                   cAACGCGCCGGTGGGCAATGCCATCGCCGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTCCCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGc                                 <  1:175162/225‑1 (MQ=255)
                                                                                                                                                                                                     aCGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAgg                                                                                                                              <  2:50686/130‑1 (MQ=34)
                                                                                                                                                                                                      cgcgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGc                                                                                                                             <  2:261572/130‑1 (MQ=33)
                                                                                                                                                                                                        cgcCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGt                                                      <  2:476602/199‑1 (MQ=21)
                                                                                                                                                                                                         gcCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTg                 <  1:957255/235‑1 (MQ=255)
                                                                                                                                                                                                          ccGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCt                                                             <  1:940576/190‑1 (MQ=38)
                                                                                                                                                                                                              tggtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTc                                                                       <  2:710236/176‑1 (MQ=33)
                                                                                                                                                                                                              tggtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCg                   >  1:1062248/1‑228 (MQ=255)
                                                                                                                                                                                                              tggtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCg                   <  2:1062245/228‑1 (MQ=255)
                                                                                                                                                                                                                        cATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCaacaac  <  1:1068438/235‑1 (MQ=38)
                                                                                                                                                                                                                        cATCGACGCAGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac  <  1:628823/235‑1 (MQ=255)
                                                                                                                                                                                                                         aTCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCCCTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGcc                                       <  1:265389/197‑1 (MQ=255)
                                                                                                                                                                                                                           cGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAAc                                            >  1:228902/1‑190 (MQ=39)
                                                                                                                                                                                                                            gACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTg                                   <  2:224041/198‑1 (MQ=37)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
GACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAC  >  AE015451/201634‑202084

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: