Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,550 T→C 15.2% 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,5500TC15.2% 216.0 / 12.4 82L1686L (TTG→CTG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (27/42);  new base C (4/9);  total (31/51)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

CATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGAC  >  AE015451/199336‑199777
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                   
cATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTg                                                                                                                                                                                                                 <  1:693682/235‑1 (MQ=255)
          gTCACCATCACGGCTGCCCAGACTTCGGTAGCCGAGACCGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTg                                                                                                                                                                                                       <  1:287391/235‑1 (MQ=255)
                           agagacagcGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAAc                                                                                                                                                                                      <  1:291595/227‑1 (MQ=255)
                            cAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGtgcgg                                                                                                                                                                                                 >  2:658734/1‑221 (MQ=255)
                                                 gTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcggcggag                                                                                                                                                                                                                 >  2:291592/1‑184 (MQ=255)
                                                 gTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCt                                                                                                                                                                              <  2:644371/221‑1 (MQ=255)
                                                  tCAAGCCGCCCTTCACCGGGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGc                                                                                                                                                                                                            <  2:350310/190‑1 (MQ=255)
                                                      gccgccCTTCCCCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTggg                                                                                                                                                                                                                                 <  2:332728/165‑1 (MQ=255)
                                                      gCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGtgc                                                                                                                                                                     >  2:772219/1‑225 (MQ=255)
                                                      gCCGACCTCCCCCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCgaga                                                                                                                                                                                        <  1:168971/206‑1 (MQ=255)
                                                              tCACCGTGCACGGCAACCAGCCGCTGGCCCACGACTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcggcg                                                                                                                                                                                                                    <  1:93448/169‑1 (MQ=255)
                                                                       aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTc                                                                                                                                         >  1:663395/1‑236 (MQ=255)
                                                                            aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCGcgac                                                                                                                                      >  1:888748/1‑234 (MQ=255)
                                                                                 gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtg                                                                                                                                >  1:1038683/1‑235 (MQ=255)
                                                                                   cGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTcgacg                                                                                                                                     >  2:1058908/1‑228 (MQ=255)
                                                                                   cGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtggtg                                                                                                                             <  1:772222/236‑1 (MQ=255)
                                                                                       ggCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCa                                                                                                                         >  1:509162/1‑236 (MQ=255)
                                                                                                 cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGAc                                                                                                                  <  2:1038680/233‑1 (MQ=255)
                                                                                                      cGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCAccc                                                                                                          >  1:93471/1‑236 (MQ=255)
                                                                                                            ccTGAGCAACAACGCCCAGGTCGCCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGtgc                                                                                                                                                                                                   <  2:113155/141‑1 (MQ=37)
                                                                                                             cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcggca                                                                                                                                                                                                                    >  1:547673/1‑122 (MQ=255)
                                                                                                             cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTcggca                                                                                                                                                                                                                    <  2:547670/123‑2 (MQ=255)
                                                                                                             cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTcggca                                                                                                                                                                                                                    >  1:600939/1‑122 (MQ=255)
                                                                                                                 gCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGAcc                                                                                               >  1:917690/1‑236 (MQ=255)
                                                                                                                        cGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTAGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCCGACTGCCACCCCTTCGGTGACCGAAg                                                                                           <  2:663392/233‑1 (MQ=255)
                                                                                                                        cGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTgccgct                                                                                                                                                                 >  2:688123/1‑163 (MQ=255)
                                                                                                                         gCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCgg                                                                                                    >  1:897261/1‑223 (MQ=255)
                                                                                                                           ccAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                      <  1:170354/235‑1 (MQ=255)
                                                                                                                           ccAGGTCACCATCAAGGCTGGCGAGACCAGCGCACCGTACACTCACGACGCGCAAGGCGATGACGTCTATCAGGACGCGGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                      <  1:391404/235‑1 (MQ=255)
                                                                                                                           ccAGGTCACCATCAAGGACGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                      <  1:305770/235‑1 (MQ=255)
                                                                                                                            ccagcgtcctgatagacgtcatcgccttgcgcggcgtgggtgtacggcgcgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCa                                                                                                                                                                                                                            <  1:600156‑M2/48‑1 (MQ=255)
                                                                                                                             aGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTgag                                                                                   >  1:44291/1‑236 (MQ=255)
                                                                                                                               gTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCaa                                                                                                                        <  2:857172/197‑1 (MQ=11)
                                                                                                                                tCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGGGGACGCCACTGGtgctg                                                                                                                                                                                                 >  2:124264/1‑123 (MQ=25)
                                                                                                                                  aCCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGCCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCa                                                                               <  1:527893/235‑1 (MQ=255)
                                                                                                                                   ccATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                              <  1:970986/235‑1 (MQ=255)
                                                                                                                                   ccATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                              <  1:400430/235‑1 (MQ=255)
                                                                                                                                    cATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                              <  2:740063/234‑1 (MQ=11)
                                                                                                                                     aTCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTg                                                                                     <  1:517182/226‑1 (MQ=255)
                                                                                                                                     aTCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGACGCGGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGtct                                                                                                                                                          >  2:391401/1‑155 (MQ=11)
                                                                                                                                          ggCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                              <  1:330244/228‑1 (MQ=255)
                                                                                                                                             gtcatcgccttGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCAt                                                                     <  2:824316‑M1/223‑1 (MQ=255)
                                                                                                                                             gtcatcgccttGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATc                                                                    >  1:824319‑M1/12‑235 (MQ=255)
                                                                                                                                                 gagaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCCTTCGAGAACCTTGAGCTggg                                                                                                                                                                           >  1:343444/1‑128 (MQ=255)
                                                                                                                                                   gaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTacac                                                                        <  2:44290/225‑1 (MQ=40)
                                                                                                                                                      cAGCGCGCCGTACACCCACGACGCGCAGGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                      <  1:899767/208‑1 (MQ=255)
                                                                                                                                                       aGCTCGCCGTACACCCACGACGCGCGAGGCGATGACGTCTATCAGGGCGCTGGCCGGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTTGTGCTGCATTCGAGAACCtt                                                                                                                                                                                   <  2:107088/114‑1 (MQ=11)
                                                                                                                                                            GCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCGGGCCAGATCAGCCTGGGCATCAACTCCGCGGTGGAcgt                                                                                                                                                                                                            <  2:637878‑M1/84‑2 (MQ=255)
                                                                                                                                                             tcgtgggtgtacggcgcccAAGGCGATGACGTCTATCAGGACGGTGGCCAGATCAGCTTGGGCATCAACTCGGCGGt                                                                                                                                                                                                                  <  2:915750‑M2/58‑1 (MQ=255)
                                                                                                                                                             tcgtgggtgtaccgcgcgcATGGCGATGTCGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcg                                                                                                                                                                                                                       <  1:485609‑M2/53‑1 (MQ=255)
                                                                                                                                                             ccGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAg                                                    <  1:813943/235‑1 (MQ=255)
                                                                                                                                                             ccGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAg                                                    <  1:526030/235‑1 (MQ=40)
                                                                                                                                                                 acacCCACGACGCGCAAGGCGAGGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAg                                                                                                                                                                                >  1:810340/1‑107 (MQ=38)
                                                                                                                                                                  cacCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTg                                                                                     <  2:509159/197‑1 (MQ=255)
                                                                                                                                                                  cacCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                              <  2:93469/204‑1 (MQ=255)
                                                                                                                                                                   acCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACACAGAc                                                  >  1:354609/1‑231 (MQ=255)
                                                                                                                                                                        cggtgcgcGAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAg                                                                                           <  1:953210‑M2/177‑1 (MQ=255)
                                                                                                                                                                           cgcgcAGGGCGATGACGTCTATCAGGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCt                                                                                                                                                                              <  2:366990‑M2/94‑1 (MQ=255)
                                                                                                                                                                           cgcgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTg                                                                                                                                                                             >  2:1088083‑M2/6‑100 (MQ=255)
                                                                                                                                                                           cgcgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGTTCa                                     >  1:136113‑M2/6‑236 (MQ=255)
                                                                                                                                                                             cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGAc                                                                                                                                                                                                              >  2:191592‑M2/4‑65 (MQ=255)
                                                                                                                                                                              gcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCg                                                                                              >  1:669936‑M2/3‑176 (MQ=255)
                                                                                                                                                                              gcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCg                                                                                              <  2:669933‑M2/174‑1 (MQ=255)
                                                                                                                                                                              gcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                              <  2:1017445‑M2/190‑1 (MQ=255)
                                                                                                                                                                                aaGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaa                                                         >  1:115937/1‑211 (MQ=39)
                                                                                                                                                                                 aGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAAcc                               <  1:1011827/236‑1 (MQ=40)
                                                                                                                                                                                 aGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                               <  1:634903/236‑1 (MQ=21)
                                                                                                                                                                                  ggCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAAcc                               >  1:631748/1‑235 (MQ=255)
                                                                                                                                                                                  ggCGATGACGTCTATCAGGACGCGGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAAcc                               <  1:269781/235‑1 (MQ=255)
                                                                                                                                                                                        gACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAAcc                               <  1:138727/229‑1 (MQ=37)
                                                                                                                                                                                         acagcTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATc                                      <  1:1047092/217‑1 (MQ=255)
                                                                                                                                                                                         aCGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCg                          <  2:917687/233‑1 (MQ=255)
                                                                                                                                                                                              gaTCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGa                                                             <  2:1059245/192‑1 (MQ=255)
                                                                                                                                                                                                 cAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCAccc                                                                >  2:701119/1‑187 (MQ=255)
                                                                                                                                                                                                  aGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAg                                                                                                                                  <  2:425608/120‑1 (MQ=255)
                                                                                                                                                                                                        gtGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGacac                                                                                                                                             <  1:578521/102‑1 (MQ=255)
                                                                                                                                                                                                         tagCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCt                                             <  2:516974/196‑1 (MQ=255)
                                                                                                                                                                                                         tGGCGAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCg                                                                                                     >  2:588002/1‑142 (MQ=255)
                                                                                                                                                                                                         tGGCCAGATCAGCTTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCTGCATTCGAGAACCTCGAGCTGGGCGGTGCTGCTTCGGTTCATGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCAcg                                                                 <  2:356750/178‑2 (MQ=11)
                                                                                                                                                                                                            ccAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGGGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGCAGGCGGTGAGATCACCTACACCATCACCCTGACcaaca                                                       >  1:237992/1‑185 (MQ=255)
                                                                                                                                                                                                             cAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCgaga                                                                                                                                                                                        <  2:987458/55‑1 (MQ=35)
                                                                                                                                                                                                                   aGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGAc  >  1:938916/1‑231 (MQ=40)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                   
CATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGAC  >  AE015451/199336‑199777

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