Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,751 G→C 6.0% T419T (ACG→ACC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,7510GC 6.0% 288.8 / 4.5 94T419T (ACG→ACCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (55/33);  new base C (3/3);  total (58/36)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATC  >  AE015451/195526‑195982
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
ttCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTg                                                                                                                                                                                                                               >  1:39820/1‑236 (MQ=16)
 tCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                                                              <  1:996755/236‑1 (MQ=255)
                                 aggcggTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTcc                                                                                                                                                                                              >  1:795760/2‑236 (MQ=255)
                                         aTTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                                                                         <  1:991877/201‑3 (MQ=2)
                                           tgcGAGCAACTGACCCCGACTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTAcc                                                                                                                                                                                                                                   <  1:661016/187‑1 (MQ=255)
                                                         ccggAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                                                                         <  1:124073/182‑3 (MQ=2)
                                                            gaacccGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                                                              >  1:631477/5‑177 (MQ=255)
                                                            gAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                                                              <  2:631475/177‑1 (MQ=255)
                                                                 cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                                                         <  1:543885/177‑1 (MQ=40)
                                                                       cGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                                                                <  2:599781/212‑1 (MQ=25)
                                                                        ggCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                                                                <  2:1238793/163‑1 (MQ=12)
                                                                                                 gTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGcc                                                                                                                                >  1:491548/1‑234 (MQ=255)
                                                                                                   cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                                                         >  1:983427/6‑143 (MQ=255)
                                                                                                   cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                                                         <  2:983425/138‑1 (MQ=255)
                                                                                                   cGACACCACTACCGCGACCCTGACGACGAGCCCGTCGGTCACCGAAGGAGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCa                                                                                                                               >  1:839929/1‑233 (MQ=255)
                                                                                                               cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACgg                                                                                                                                                               >  1:1253799/1‑189 (MQ=37)
                                                                                                               cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTc                                                                                                                 >  1:1409269/1‑235 (MQ=255)
                                                                                                                      cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAACCGACTAcctc                                                                                                                                                                                                                                 >  2:493501/1‑114 (MQ=255)
                                                                                                                      cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGa                                                                                                                                                                              >  1:55301/1‑167 (MQ=38)
                                                                                                                      cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGc                                                                                                          >  1:787387/1‑235 (MQ=255)
                                                                                                                            gCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                                                                         >  1:389806/1‑150 (MQ=39)
                                                                                                                                 cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCTGTACCCTAAACACATCTGCCCACCCACCAGTAACCGTGACACTGTCCAACCCCCAAACGATTACCGTAGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACggtt                                                                                                                                                             >  1:705607/1‑171 (MQ=11)
                                                                                                                                   cGTCGGTGACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                                                         <  2:797562/111‑1 (MQ=255)
                                                                                                                                      cGGTCACCGAAGGGGGCGTGATCACTTACCCCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                                                  <  2:313529/147‑1 (MQ=255)
                                                                                                                                        gtgaCCGAAGGTGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt                                                                                                                                                                           <  2:175697/149‑1 (MQ=16)
                                                                                                                                        gtgaCCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACa                                                                                                            <  1:64729/212‑1 (MQ=9)
                                                                                                                                         tCACCGAAGGTGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGt                                                                                                                                                                                                      <  2:772713/124‑1 (MQ=255)
                                                                                                                                         tCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAGCGATCACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCa                                                                                                                                                    <  2:378212/174‑1 (MQ=12)
                                                                                                                                                        tGATCACTTACACCGTGATCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                                                         <  1:1479188/90‑1 (MQ=255)
                                                                                                                                                        tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGGCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGAcc                                                                                                                                                        >  1:464821/1‑155 (MQ=255)
                                                                                                                                                        tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGCCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTc                                                                                                                                                            >  1:223860/1‑151 (MQ=255)
                                                                                                                                                         gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATc                                                                       >  1:739644/1‑235 (MQ=255)
                                                                                                                                                          aTCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGac                                                                                                              >  1:1139372/1‑194 (MQ=37)
                                                                                                                                                                     cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCg                                                                                                                                                                                                                          <  2:1041621/76‑1 (MQ=255)
                                                                                                                                                                     cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGTCCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCg                                                                                            <  1:220796/202‑1 (MQ=255)
                                                                                                                                                                        gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGa                                                                                                               >  2:487983/1‑180 (MQ=9)
                                                                                                                                                                        gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                                                                <  1:265213/67‑1 (MQ=255)
                                                                                                                                                                        gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                                                                >  2:265211/1‑67 (MQ=255)
                                                                                                                                                                          cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgagacc                                                                                                                                   >  2:40758/1‑156 (MQ=38)
                                                                                                                                                                            cTGAGCAATCCTGCCCAGACGCCGGTAACAGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCCGGGCAGCGTCGATTTCAAGACCCCGGCCAATGACGTCTACAACAACGGCTCg                                                                                                                                                           >  2:723497/1‑132 (MQ=11)
                                                                                                                                                                            cTGAGCAACGCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccac                                                   <  1:1406306/236‑1 (MQ=255)
                                                                                                                                                                            cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAACCATTAc                                                                                                                                                                                                                                    >  2:1082816/1‑59 (MQ=11)
                                                                                                                                                                            cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                                                            >  1:1007246/1‑99 (MQ=255)
                                                                                                                                                                            cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGaatc                                                                                              >  2:488177/1‑191 (MQ=2)
                                                                                                                                                                            cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCCGCGTCCCCAt                                                                                                                                          >  2:465546/1‑149 (MQ=255)
                                                                                                                                                                            cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCGCGACCGTCAGCGTCACCATCgagaac                                                                                                                                   >  2:420081/1‑154 (MQ=255)
                                                                                                                                                                               aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac                                                <  1:487985/236‑1 (MQ=255)
                                                                                                                                                                               aGCAACCCTGCCCAGACGCCGGTCACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                                                            <  2:1444158/96‑1 (MQ=12)
                                                                                                                                                                               aGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGt                                                                                                                                                                                                      <  2:1422964/86‑1 (MQ=255)
                                                                                                                                                                               aGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTgg                                                                                                                           <  2:510574/161‑1 (MQ=40)
                                                                                                                                                                                gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                                                         <  2:334991/66‑1 (MQ=255)
                                                                                                                                                                                gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCCTc                                                                                                                                         >  2:391955/1‑146 (MQ=255)
                                                                                                                                                                                gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac                                                >  1:1404849/1‑235 (MQ=255)
                                                                                                                                                                                      cTGCCAAGACGCCGGTAACCGTGATCCTGACCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGACGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGt                                                            <  1:385177/217‑1 (MQ=255)
                                                                                                                                                                                       tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                                                                >  2:1362123/1‑52 (MQ=255)
                                                                                                                                                                                       tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                                                              >  1:1444568/1‑54 (MQ=255)
                                                                                                                                                                                       tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                                                              <  2:1444565/54‑1 (MQ=255)
                                                                                                                                                                                       tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAcc                                                                             <  2:674289/199‑1 (MQ=255)
                                                                                                                                                                                       tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAcc                                                                             >  1:674291/1‑199 (MQ=255)
                                                                                                                                                                                       tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCAACTTCGaa                                                                                                              >  1:415499/1‑166 (MQ=25)
                                                                                                                                                                                       tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGGCTACAACAACGGCTCGAACGTCAGCGTCACCATCga                                                                                                                                       >  1:5260/1‑141 (MQ=255)
                                                                                                                                                                                       tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTccc                                                                                                                                                                                             >  1:654066/1‑86 (MQ=255)
                                                                                                                                                                                        gCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGTAACTTCGac                                                                                                              >  2:1017523/1‑164 (MQ=255)
                                                                                                                                                                                            aGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTg                                                                                                                                                                                                                               <  1:1132635/48‑1 (MQ=255)
                                                                                                                                                                                            aGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                                                            >  1:193682/1‑83 (MQ=38)
                                                                                                                                                                                              ccGCCGGTAACCGTGCCCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAc                                                                              <  1:234641/190‑1 (MQ=255)
                                                                                                                                                                                               cGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGac                                                        <  1:1225827/212‑1 (MQ=255)
                                                                                                                                                                                                 ccGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTCAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCAcaa                                                                         >  2:1013454/1‑191 (MQ=255)
                                                                                                                                                                                                 ccGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGGAATTTc                                                                                                                 >  1:202174/1‑153 (MQ=255)
                                                                                                                                                                                                 ccGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCgg                                                                                                                                                                                                                                 >  2:326101/1‑40 (MQ=255)
                                                                                                                                                                                                    agacccGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCg                                                                                      <  1:142014/173‑1 (MQ=255)
                                                                                                                                                                                                      aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTATAACCCAGGGCAGCGTCGATTTCCAGTCCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTggcggc                                                                                                                       >  2:927946/1‑142 (MQ=255)
                                                                                                                                                                                                      aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                                                                         >  1:641557/1‑76 (MQ=40)
                                                                                                                                                                                                      aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACgga                                                                                                                                                              >  1:294919/1‑102 (MQ=25)
                                                                                                                                                                                                      aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAACAACTGACCCCGAACCCGACGCCGGCCCAGACCACGATCAACGACt                                                                >  2:1173711/1‑197 (MQ=255)
                                                                                                                                                                                                      aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCCACaacaac                                                                                                                                                                 >  2:1365617/1‑100 (MQ=255)
                                                                                                                                                                                                            gAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCt                                                                                                                                                             >  1:771458/1‑98 (MQ=35)
                                                                                                                                                                                                            gAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCg                                                              <  1:456095/193‑1 (MQ=18)
                                                                                                                                                                                                               ccTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCggcc                    >  1:980320/1‑230 (MQ=255)
                                                                                                                                                                                                               ccTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                 >  1:1456977/1‑235 (MQ=255)
                                                                                                                                                                                                               ccTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGGCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGAGTCACCATCGAGAACGCCACCGGCGGCa                                                                                                                      >  2:234640/1‑134 (MQ=11)
                                                                                                                                                                                                                 tGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCg                                                                                            >  2:728696/1‑158 (MQ=255)
                                                                                                                                                                                                                 tGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATAACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCg                                                                                                                                        >  1:288462/1‑114 (MQ=17)
                                                                                                                                                                                                                   ataagagacagaaCGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAgg            <  1:730350/225‑1 (MQ=255)
                                                                                                                                                                                                                      aaCGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAAccc                                                                                             >  2:1382044/1‑152 (MQ=40)
                                                                                                                                                                                                                      aaCGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAACTGAccc                                                                                                    >  2:713340/1‑145 (MQ=25)
                                                                                                                                                                                                                      aaCGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGGCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCCCGACCGTCAGCGTCACCACCga                                                                                                                                       >  2:255404/1‑110 (MQ=255)
                                                                                                                                                                                                                       aCGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc        <  1:269954/236‑1 (MQ=255)
                                                                                                                                                                                                                        cGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACCggcggct                                                                                                                      >  2:142013/1‑124 (MQ=17)
                                                                                                                                                                                                                             aaaCGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc  >  1:1032255/1‑236 (MQ=255)
                                                                                                                                                                                                                              aaCGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCAACTTCGAGCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc  <  1:1216112/235‑1 (MQ=255)
                                                                                                                                                                                                                                cGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGCAggcgg         >  1:1215160/1‑226 (MQ=255)
                                                                                                                                                                                                                                cGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACcg                                              >  2:663872/1‑189 (MQ=255)
                                                                                                                                                                                                                                cGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAc                                                                              >  1:22608/1‑157 (MQ=25)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
TTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATC  >  AE015451/195526‑195982

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