Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,550 T→C 8.1% L1686L (TTG→CTG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451199,5500TC 8.1% 236.4 / 3.5 78L1686L (TTG→CTG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base C (3/4);  total (39/39)
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

CGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACG  >  AE015451/199327‑199778
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
cGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACt                                                                                                                                                                                                                          >  1:380787/1‑236 (MQ=255)
                      aCCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCCTCAACTCGGCGGTg                                                                                                                                                                                                                  >  1:319272/1‑222 (MQ=255)
                          tgaCGGCGGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTg                                                                                                                                                                                                        <  1:736310/226‑1 (MQ=255)
                          tCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGtgctgc                                                                                                                                                                                                 >  1:527336/1‑235 (MQ=255)
                                           tCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACCcg                                                                                                                                                                                                                        >  1:855740/1‑195 (MQ=255)
                                            cGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCt                                                                                                                                                                               <  1:463922/235‑1 (MQ=255)
                                              gTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACACCGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCAtt                                                                                                                                                                                              <  1:579875/218‑1 (MQ=255)
                                                aGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCg                                                                                                                                                                          <  1:589842/236‑1 (MQ=255)
                                                     agaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGtgct                                                                                                                                                                     >  1:933278/1‑236 (MQ=255)
                                                                 cGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGt                                                                                                                                                                                                                   <  1:267011/178‑1 (MQ=255)
                                                                                aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCt                                                                                                                                           >  1:695334/1‑235 (MQ=255)
                                                                                        cAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgt                                                                                                                                  >  1:838150/1‑236 (MQ=255)
                                                                                          gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTg                                                                                                                                                                                                        <  2:467463/164‑1 (MQ=40)
                                                                                          gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTg                                                                                                                                                                                                        >  1:467464/1‑164 (MQ=39)
                                                                                            cGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtggtg                                                                                                                              <  1:1073136/236‑1 (MQ=255)
                                                                                             gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACATCGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCt                                                                                                                                                                               <  1:480871/186‑1 (MQ=255)
                                                                                                 gCCCCCGACCCGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCaa                                                                                                                         <  1:265653/236‑1 (MQ=255)
                                                                                                   ccACGACCTGGTCGTGACCCTGAGCAACACCGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTcg                                                                                                                                         <  1:667672/218‑1 (MQ=255)
                                                                                                      cGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGt                                                                                                                                                                        >  1:863725/1‑184 (MQ=255)
                                                                                                       tgccTGGTCGTGACCCTGATCAACCACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGCCGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCa                                                                                                                                                     <  2:318102/200‑1 (MQ=255)
                                                                                                       gACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGAc                                                                                                                   <  1:790113/236‑1 (MQ=255)
                                                                                                          cTGGGCGTGACGCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCACCGTACACCCACGACGCGCACGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCg                                                                                                                                                               <  1:209562/189‑1 (MQ=255)
                                                                                                           tGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTAACAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCAc                                                                                                                                                    <  1:39188/199‑1 (MQ=11)
                                                                                                           tGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTg                                                                                                                 <  2:933277/234‑1 (MQ=255)
                                                                                                                     ccTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGAc                                                                                                                                                                                                               >  2:121303/1‑130 (MQ=38)
                                                                                                                      cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGCCGCTGGCCAGATCGGCCTGGGCGTCAACTcggcg                                                                                                                                                                                                                     >  2:235373/1‑123 (MQ=11)
                                                                                                                          gCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGaa                                                                                                                                    >  1:439688/1‑200 (MQ=255)
                                                                                                                           caacaaCGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGAcc                                                                                                >  1:732291/1‑235 (MQ=255)
                                                                                                                           caacaaCGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTc                                                                                                                                          >  1:937952/1‑193 (MQ=255)
                                                                                                                               aacGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGtg                                                                                                                                                                                                     >  2:550983/1‑130 (MQ=38)
                                                                                                                                 cGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCAcc                                                                                                            >  2:1102011/1‑217 (MQ=255)
                                                                                                                                  gCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAgg                                                                                           >  2:550108/1‑233 (MQ=255)
                                                                                                                                    ccAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTg                                                                                      <  1:858841/236‑1 (MQ=255)
                                                                                                                                           aCCAGCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATc                                                                                 <  2:607978/234‑1 (MQ=255)
                                                                                                                                            ccATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGGGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtgg                                                                                                                                <  2:930516/186‑1 (MQ=255)
                                                                                                                                                cAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGGGGACGCCACTggg                                                                                                                                                                                                      >  2:253452/1‑111 (MQ=25)
                                                                                                                                                          gagaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTgaga                                                                                   >  2:486069/1‑217 (MQ=18)
                                                                                                                                                          gagaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTTCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCa                                                                                                                          >  1:453280/1‑178 (MQ=255)
                                                                                                                                                           agaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGtg                                                                                                                                                                                                     >  2:469010/1‑102 (MQ=31)
                                                                                                                                                              cttgcgcgcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGAc                                                                                                                   <  2:1009931/178‑1 (MQ=255)
                                                                                                                                                                aGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGAcc                                                            >  2:679599/1‑234 (MQ=255)
                                                                                                                                                                aGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACca                                                           >  1:965008/1‑235 (MQ=40)
                                                                                                                                                                   gcgcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGa                                                                                              <  1:767805/197‑1 (MQ=255)
                                                                                                                                                                     gcCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGACGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaacaa                                                       <  1:481113/234‑1 (MQ=255)
                                                                                                                                                                     gcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                       <  1:1086596/202‑1 (MQ=255)
                                                                                                                                                                     gcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                       >  2:1086595/1‑202 (MQ=255)
                                                                                                                                                                            ccccACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGGTCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGGGAACCTTGAGCt                                                                                                                                                                               <  2:193222/106‑1 (MQ=255)
                                                                                                                                                                            acCCACGACGCGCAAGGCGATGCCGGCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCGACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAg                                                                                                                        <  2:14077/162‑1 (MQ=255)
                                                                                                                                                                             cccACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTCCAtt                                                                                                                                                                                              <  2:178739/91‑1 (MQ=17)
                                                                                                                                                                               cACGACGCGCACGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAGCGCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCGTCGGTGACCGAAGGCGGTGAGATCACCTa                                                                            <  1:360268/203‑1 (MQ=255)
                                                                                                                                                                                    cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTTCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCAc                                                                   >  2:445703/1‑207 (MQ=37)
                                                                                                                                                                                    cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCg                                          >  2:679270/1‑232 (MQ=255)
                                                                                                                                                                                     gcgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCaa                                     <  1:519791/236‑1 (MQ=39)
                                                                                                                                                                                      cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTgag                                                                                    >  1:30838/1‑188 (MQ=40)
                                                                                                                                                                                      cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAAt                                    >  1:954896/1‑236 (MQ=255)
                                                                                                                                                                                      cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTATACGCCACTGGTGCTGCATTCGAGCACGTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGaa                                                                                                                                    >  2:627025/1‑140 (MQ=255)
                                                                                                                                                                                      cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCAGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAcc                                                                              >  2:997735/1‑194 (MQ=255)
                                                                                                                                                                                          aGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTgag                                                                                    <  2:873640/184‑1 (MQ=40)
                                                                                                                                                                                           ggCGATGACGTCTATCAGGACGCTGGCCAGAGCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTga                                                                                     <  2:370133/182‑1 (MQ=255)
                                                                                                                                                                                             cGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCAc                              <  1:557631/235‑1 (MQ=255)
                                                                                                                                                                                                   cGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                               <  2:249826/180‑1 (MQ=255)
                                                                                                                                                                                                     tCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTg                     <  1:894247/236‑1 (MQ=255)
                                                                                                                                                                                                          cAGGATGCGGGCCAGATCAGCCTGGGCATCAGCTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTACGGTTGAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAAcc                                <  1:829052/220‑1 (MQ=255)
                                                                                                                                                                                                          cAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTTCTTCGGTTCGAGTCACCGACACCCTCGATGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACt                 >  1:616173/1‑235 (MQ=255)
                                                                                                                                                                                                          cAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATACCCACTCGGAGCTGTACtt                >  1:33950/1‑236 (MQ=255)
                                                                                                                                                                                                               cGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaac                                                         >  1:516259/1‑190 (MQ=40)
                                                                                                                                                                                                                 cTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGc                                                                                          >  2:1010896/1‑155 (MQ=255)
                                                                                                                                                                                                                 cTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGc                                                                                          <  1:1010897/155‑1 (MQ=255)
                                                                                                                                                                                                                    gCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCAcc                                                                  >  1:345845/1‑176 (MQ=255)
                                                                                                                                                                                                                    gCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGa                         <  1:757019/217‑1 (MQ=38)
                                                                                                                                                                                                                    gCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGa                         >  2:757018/1‑217 (MQ=38)
                                                                                                                                                                                                                    gCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCa                                                                                                                                                          <  2:458430/88‑1 (MQ=255)
                                                                                                                                                                                                                    gCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTCGAGCTGGGCGGTGCTGCTTCAGTTCa                                                                                                                                                          >  1:458431/1‑88 (MQ=14)
                                                                                                                                                                                                                    gCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGATGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTaca                                                                          >  1:860731/1‑168 (MQ=255)
                                                                                                                                                                                                                     ccAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                               <  2:64130/194‑1 (MQ=40)
                                                                                                                                                                                                                      cAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCCGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGa    <  1:445704/236‑1 (MQ=18)
                                                                                                                                                                                                                      cAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGGGCTGCATTCgag                                                                                                                                                                                          >  2:152254/1‑54 (MQ=255)
                                                                                                                                                                                                                        gATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCGTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACg  <  1:281011/236‑1 (MQ=37)
                                                                                                                                                                                                                               ttGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGt                                                                                                    <  2:1014299/131‑1 (MQ=18)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
CGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACG  >  AE015451/199327‑199778

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