Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 200,326 C→G 5.3% A1944A (GCC→GCG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451200,3260CG 5.3% 295.9 / 3.0 87A1944A (GCC→GCGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (35/47);  new base G (2/3);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCA  >  AE015451/200102‑200557
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                       
cTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCg                                                                                                                                                                                                                                >  2:787857/1‑234 (MQ=255)
           tGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCg                                                                                                                                                                                                                   >  1:464766/1‑236 (MQ=255)
             cTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                                                                                                                                                 <  1:601393/236‑1 (MQ=255)
              tCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTg                                                                                                                                                                                                                <  1:715507/236‑1 (MQ=255)
                aGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTg                                                                                                                                                                                                                >  1:549061/1‑234 (MQ=255)
                 ggTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTg                                                                                                                                                                                                                <  2:558327/233‑1 (MQ=255)
                      ccATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCa                                                                                                                                                                                                          <  1:108828/234‑1 (MQ=255)
                        aTCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                                                                                                                                                         <  2:703080/233‑1 (MQ=255)
                        aTCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTa                                                                                                                                                                                                      >  1:51290/1‑236 (MQ=255)
                         tCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTac                                                                                                                                                                                                     <  1:813403/236‑1 (MQ=255)
                               cTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCAt                                                                                                                                                                                                >  2:658657/1‑235 (MQ=255)
                               cTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATc                                                                                                                                                                                               <  1:793036/236‑1 (MQ=255)
                                       acaggcaggcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                                                                                                                                                 <  1:360023/202‑1 (MQ=255)
                                             gcgcCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTg                                                                                                                                                                                                                <  2:1254176/205‑1 (MQ=255)
                                               gcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCt                                                                                                                                                                                                                                   >  2:872521/1‑184 (MQ=25)
                                               gcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCt                                                                                                                                                                                                                                   <  1:872526/184‑1 (MQ=255)
                                                     cacCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaa                                                                                                                                                                                    >  2:407484/1‑225 (MQ=40)
                                                     cacCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                                         <  1:653022/236‑1 (MQ=255)
                                                     cacCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                                         <  1:888600/236‑1 (MQ=40)
                                                     cacCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                                         <  1:396458/236‑1 (MQ=40)
                                                      acCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGAc                                                                                                                                                                                                                           <  1:258053/185‑1 (MQ=39)
                                                         cACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAg                                                                                                                                                                                                                      <  1:449921/187‑1 (MQ=255)
                                                             acgcgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGAcc                                                                                                                                                                                                                          >  2:360018/3‑179 (MQ=255)
                                                             aTGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                                                                                                                                                 <  1:568054/188‑1 (MQ=255)
                                                              tGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGt                                                                                                                                                                                                                 >  1:459376/1‑187 (MQ=14)
                                                              tGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGt                                                                                                                                                                                                                 <  2:459371/187‑1 (MQ=14)
                                                                   aaGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGcc                                                                                                                                                                                                                                         <  2:335797/158‑1 (MQ=255)
                                                                           gACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGa                                                                                                                                                   <  1:319498/236‑1 (MQ=14)
                                                                            aCGGCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                                                                                                                                                 <  1:870026/173‑1 (MQ=255)
                                                                             cGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGAGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCGAAGCTGACTGCCACCCCTTCGGTGACCg                                                                                                                                                                                                                         >  1:244046/1‑164 (MQ=255)
                                                                              gTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGa                                                                                                                                                   <  2:482943/233‑1 (MQ=40)
                                                                               tCTATCAGGACGCTGGCCAGATCAGCCTCGGCATCACCTCGGCAGTAGACGTCGCTGGTGCCGCGTTCGAGAACCTTGAGCTGGGCGGCTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCAc                                                                                                                                                                                                         <  2:1138321/178‑1 (MQ=14)
                                                                                   tCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACt                                                                                                                                           <  1:1278760/236‑1 (MQ=255)
                                                                                    cAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACggcc                                                                                                                                                                         >  2:896740/1‑203 (MQ=255)
                                                                                    cAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGAGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGCAGGCGGTGAGATCACCTACACCAt                                                                                                                                                                                                >  2:188802/1‑182 (MQ=255)
                                                                                     aGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                                         <  2:1020439/204‑1 (MQ=40)
                                                                                           cTGGCCAGATCCGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGGGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCa                                                                                                                                                                                                          >  2:745595/1‑165 (MQ=255)
                                                                                            tGGCCAGATCAGCTTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTg                                                                                                                                   <  1:12962/235‑1 (MQ=18)
                                                                                            tGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCAccc                                                                                                                                                                                           >  1:1138326/1‑179 (MQ=255)
                                                                                             ggCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGAACAACAACCACTCGGAGCTGTACTTCaaac                                                                                                                                     >  2:1127052/1‑230 (MQ=255)
                                                                                              gCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGAt                                                                                                                                                                  >  1:903616/1‑202 (MQ=255)
                                                                                              gCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGAcc                                                                                                                                <  1:1172006/236‑1 (MQ=255)
                                                                                                    tCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCg                                                                                                                                                                    >  2:1137094/1‑194 (MQ=40)
                                                                                                           gggCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTg                                                                                                                                                                       <  2:633946/184‑1 (MQ=40)
                                                                                                           gggCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTg                                                                                                                                                                       >  1:633951/1‑184 (MQ=255)
                                                                                                                tCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaac                                                                                                                                                                                   >  2:51711/1‑167 (MQ=40)
                                                                                                                tCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGcaccac                                                                                                                     >  2:482576/1‑229 (MQ=255)
                                                                                                                tCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACCGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGAcc                                                                                                                                                                                      >  2:638531/1‑164 (MQ=255)
                                                                                                                tCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGCGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACt                                                                                                                                                       >  1:163129/1‑195 (MQ=255)
                                                                                                                    cTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTg                                                                                                           <  1:1241173/235‑1 (MQ=255)
                                                                                                                     tCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTc                                                                                                                                                      >  2:453969/1‑191 (MQ=40)
                                                                                                                         cAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGCGGCAGc                                                                                                      >  1:318983/1‑235 (MQ=39)
                                                                                                                          aGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCgg                                                                                                                                                    >  2:884901/1‑188 (MQ=40)
                                                                                                                            tcgACGCCACTGGTGCGGCGGTCGAGAACCTTGAGCTGGGCGGTGCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTaca                                                                                                                                                                                                    <  2:141645/136‑1 (MQ=255)
                                                                                                                              gACGCCACTGGTGCTGCGATCGAGAACCTTGAGTTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                                         <  2:710047/163‑1 (MQ=25)
                                                                                                                                      tGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCaa                                                                                                   <  2:475530/225‑1 (MQ=37)
                                                                                                                                      tGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGCACCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCTTCAGTAACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGACcaacaaacc                                                                                                                                                         >  1:343925/1‑168 (MQ=11)
                                                                                                                                              cGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                          <  2:244097/162‑1 (MQ=37)
                                                                                                                                              cGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGcc                                                                                >  1:332416/1‑236 (MQ=40)
                                                                                                                                                ttCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGa                                                                                                                                                                                                                            <  1:670735/94‑1 (MQ=38)
                                                                                                                                                ttCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGcc                                                                                >  2:1028915/1‑234 (MQ=40)
                                                                                                                                                ttCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTAGCCAAGCTGACTGCCACCCCTTCGGTGa                                                                                                                                                                                                                            >  2:670730/1‑94 (MQ=255)
                                                                                                                                                  cGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTg                                                                            <  1:407489/236‑1 (MQ=40)
                                                                                                                                                  cGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTg                                                                            <  1:325764/236‑1 (MQ=40)
                                                                                                                                                   gagaACCTTGAGCTGGGCGGTCCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGc                                                                           <  1:1137007/236‑1 (MQ=38)
                                                                                                                                                      aaCCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAAc                                                                        >  1:715063/1‑236 (MQ=255)
                                                                                                                                                         cTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAAcc                                                                       <  2:910442/234‑1 (MQ=40)
                                                                                                                                                             aGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                 <  1:658662/236‑1 (MQ=255)
                                                                                                                                                              gCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCAc                                                                                                                                                        <  2:135861/148‑1 (MQ=37)
                                                                                                                                                              gCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCAc                                                                                                                                                        >  1:135864/1‑148 (MQ=37)
                                                                                                                                                              gCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                 <  2:876184/235‑1 (MQ=255)
                                                                                                                                                              gCTGGGCGGTGCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                 <  2:1027667/235‑1 (MQ=39)
                                                                                                                                                                     ggTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGt                                                         >  1:1121965/1‑236 (MQ=40)
                                                                                                                                                                     ggTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCAACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGt                                                         >  1:750129/1‑236 (MQ=255)
                                                                                                                                                                        tCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCa                                                                                                        <  1:231035/186‑1 (MQ=9)
                                                                                                                                                                          cGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCAcc                                                                                                                                                                                            >  2:561310/1‑100 (MQ=31)
                                                                                                                                                                          cGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCAcc                                                                                                                                                                                            <  1:561315/100‑1 (MQ=31)
                                                                                                                                                                               cGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                          <  2:774591/129‑1 (MQ=33)
                                                                                                                                                                               cGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                          >  1:774596/1‑129 (MQ=33)
                                                                                                                                                                                   tCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCa                                           <  1:913359/236‑1 (MQ=38)
                                                                                                                                                                                    cAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGtt                                                                                     >  2:962230/1‑193 (MQ=17)
                                                                                                                                                                                    cAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCa                                           <  1:1192802/235‑1 (MQ=37)
                                                                                                                                                                                        taaCCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCgg                                                <  2:1166880/224‑1 (MQ=39)
                                                                                                                                                                                            cGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                          <  2:324092/116‑1 (MQ=32)
                                                                                                                                                                                                  ccTCGACGAAGTGGTGGCCAAGCTGGCTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCAGCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGt                                                                                      >  2:657213/1‑178 (MQ=25)
                                                                                                                                                                                                         gaAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCa                      <  1:962235/235‑1 (MQ=38)
                                                                                                                                                                                                                  gCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGt            <  1:787862/236‑1 (MQ=40)
                                                                                                                                                                                                                          gACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcggcg    <  1:482581/236‑1 (MQ=40)
                                                                                                                                                                                                                               ccACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCACCGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCa  >  2:1210037/1‑233 (MQ=37)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                       
CTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCA  >  AE015451/200102‑200557

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