Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,550 T→C 5.8% 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 5.8% 283.9 / 4.1 85L1686L (TTG→CTG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (32/47);  new base C (2/3);  total (35/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 = 9.63e-01

AGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACC  >  AE015451/199330‑199774
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                
aGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcg                                                                                                                                                                                                                    <  1:579184/235‑1 (MQ=255)
      cATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTg                                                                                                                                                                                                              <  1:975182/235‑1 (MQ=255)
        tCGCCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGt                                                                                                                                                                                                               <  2:570373/232‑1 (MQ=255)
          gTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGc                                                                                                                                                                                                         <  1:748636/236‑1 (MQ=255)
                              tGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCgaga                                                                                                                                                                                     <  1:728801/236‑1 (MQ=255)
                                  cAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGtt                                                                                                                                                                                          >  1:507044/1‑227 (MQ=19)
                                  cAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGtgctgcgtt                                                                                                                                                                                          <  2:507041/227‑4 (MQ=0)
                                        tCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAg                                                                                                                                                                             >  1:1156324/1‑234 (MQ=255)
                                                                   ttCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGAGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCAc                                                                                                                                                <  1:1209977/236‑1 (MQ=255)
                                                                     cACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGtt                                                                                                                                                        <  2:1226763/226‑1 (MQ=255)
                                                                       ccGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCAtt                                                                                                                                                                                          <  1:807544/190‑1 (MQ=255)
                                                                         gTGCACGTCAACCAGCTGCTGGCCCACGACCTGGTTGTGACCCTGAGCAACAACGCCCAGTTCACCATCAACGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGacac                                                                                                                                          <  1:674992/236‑1 (MQ=255)
                                                                               gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCg                                                                                                                                                                      >  2:963618/1‑202 (MQ=255)
                                                                                  aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcggcg                                                                                                                                                                                                                 <  1:1126615/156‑1 (MQ=40)
                                                                                  aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGCGCTGGGCGGTGCCGCTTCGGTTCa                                                                                                                                                      >  2:823120/1‑215 (MQ=255)
                                                                                     cAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGaa                                                                                                                                <  2:1156321/234‑1 (MQ=255)
                                                                                      aGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAg                                                                                                                               <  1:292348/234‑1 (MQ=255)
                                                                                          gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACt                                                                                                                                                                                                                      <  1:1280864/143‑1 (MQ=39)
                                                                                              gctcACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcggcgg                                                                                                                                                                                                                >  1:386214/4‑145 (MQ=255)
                                                                                              gCCCACGACCTGGTCGTGACCCTGCGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGCCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCTTGGGCATCAACt                                                                                                                                                                                                                      <  2:1174839/139‑1 (MQ=255)
                                                                                              gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGggctgcct                                                                                                                                                                                           >  1:182688/1‑164 (MQ=255)
                                                                                                 cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGacac                                                                                                                                          <  1:489628/212‑1 (MQ=255)
                                                                                                   cGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGtgctg                                                                                                                                                                <  1:762864/188‑1 (MQ=255)
                                                                                                   cGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGtgctg                                                                                                                                                                >  2:762861/1‑188 (MQ=255)
                                                                                                     aCCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAGACCAGCGCACCGTACACCCACGACGCGCAGGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCa                                                                                                                                                      <  1:309714/196‑1 (MQ=255)
                                                                                                       cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCt                                                                                                                  >  2:883545/1‑230 (MQ=255)
                                                                                                        gggTCGTGACCCTGAGCAACAACGCCCAGGTCACCAGCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCAtt                                                                                                                                                                                          <  2:366058/156‑1 (MQ=255)
                                                                                                          gTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGcc                                                                                                                       <  2:940177/222‑1 (MQ=255)
                                                                                                            cGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGtgctg                                                                                                                                                                <  1:374860/179‑1 (MQ=255)
                                                                                                             gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGCCGCCACTGGTGCTGCAtt                                                                                                                                                                                          <  1:340572/152‑1 (MQ=38)
                                                                                                                   cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGt                                                                                                >  1:1148741/1‑236 (MQ=255)
                                                                                                                        caacaaCGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCGTGGGCATCAACTCGGCGGCGGACGCCCCTGGGGCTGCATTCGAGAACCTTGa                                                                                                                                                                              >  1:271068/1‑153 (MQ=255)
                                                                                                                         aacaacGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTgg                                                                                                                                                                                                   >  2:374857/1‑131 (MQ=37)
                                                                                                                          acaacGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGAcc                                                                                                                                                                                                          >  2:340571/1‑122 (MQ=38)
                                                                                                                            aacGACCAGGGCACCATCAAGGATGGGTAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcggcg                                                                                                                                                                                                                 <  1:498970/114‑1 (MQ=14)
                                                                                                                                   aGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGa                                                                                                                                                                                                            >  2:351479/1‑112 (MQ=35)
                                                                                                                                   aGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtgg                                                                                                                            >  2:479880/1‑192 (MQ=255)
                                                                                                                                        aCCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                           <  1:739234/236‑1 (MQ=255)
                                                                                                                                          cctaaaGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGGCTATCAGTGCGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGa                                                                                                                                             <  2:372357/164‑1 (MQ=255)
                                                                                                                                          cATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtggtg                                                                                                                          <  1:177768/187‑1 (MQ=255)
                                                                                                                                                  cTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAGCTCGGCGGTGGACGc                                                                                                                                                                                                         <  2:348055/100‑1 (MQ=25)
                                                                                                                                                  cTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATc                                                                 <  1:883548/236‑1 (MQ=255)
                                                                                                                                                  cTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCAACATTCga                                                                                                                                                                                       >  1:1176910/1‑118 (MQ=255)
                                                                                                                                                    ggTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCgg                                                                                    >  2:1120539/1‑215 (MQ=255)
                                                                                                                                                    ggTGAGACCACCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                           <  2:1121949/224‑1 (MQ=255)
                                                                                                                                                      gtagaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                           <  1:250131/220‑1 (MQ=255)
                                                                                                                                                       gagaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGGGCTGCTTCGGGTCAAGTCACCgg                                                                                                                                             >  2:177768/1‑154 (MQ=255)
                                                                                                                                                       gagaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCCGGCCAGATCAGCTTGGGCATCAACTCGGCGGCGGAc                                                                                                                                                                                                           >  2:366173/1‑93 (MQ=16)
                                                                                                                                                                  gcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACaaa                                                  >  1:419288/1‑235 (MQ=255)
                                                                                                                                                                       acacCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTa                                                                        >  1:1168643/1‑208 (MQ=255)
                                                                                                                                                                           ccACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                   <  2:155559/193‑1 (MQ=255)
                                                                                                                                                                            cACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGAc                                                                                                               <  2:518205/164‑1 (MQ=255)
                                                                                                                                                                              cgacgCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGAcc                                                        >  1:636461/1‑217 (MQ=40)
                                                                                                                                                                               gacgCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                           <  2:1159089/197‑1 (MQ=255)
                                                                                                                                                                                 cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCgaga                                                                                                                                                                                     <  1:887335/89‑1 (MQ=33)
                                                                                                                                                                                 cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCgaga                                                                                                                                                                                     >  2:887332/1‑89 (MQ=33)
                                                                                                                                                                                 cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtggtg                                                                                                                          <  2:737776/148‑1 (MQ=255)
                                                                                                                                                                                 cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAg                                                 >  2:1092808/1‑221 (MQ=255)
                                                                                                                                                                                 cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAg                                                 <  1:1092811/221‑1 (MQ=255)
                                                                                                                                                                                   cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                   <  2:870718/185‑1 (MQ=40)
                                                                                                                                                                                   cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCaac                                >  1:72387/1‑235 (MQ=38)
                                                                                                                                                                                    gcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGtgc                                                                                                                                                                  >  2:262061/1‑105 (MQ=255)
                                                                                                                                                                                      aaGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtg                                                                                                                             >  2:949994/1‑140 (MQ=255)
                                                                                                                                                                                      aaGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAAc                             >  1:1098179/1‑236 (MQ=40)
                                                                                                                                                                                       aGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGCGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                           <  1:870648/221‑1 (MQ=39)
                                                                                                                                                                                             tGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCgg                      <  1:1120542/236‑1 (MQ=255)
                                                                                                                                                                                               aCGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGa                     <  1:531552/235‑1 (MQ=255)
                                                                                                                                                                                                cGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTga                                                                                 <  2:1168640/174‑1 (MQ=255)
                                                                                                                                                                                                cGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACt                         >  1:1168167/1‑230 (MQ=255)
                                                                                                                                                                                                     aTCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGa                                                                                                                >  2:93215/1‑138 (MQ=255)
                                                                                                                                                                                                       cAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCCCTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGGGGCCAAGCTGACTGCCACCCCCTCGGt                                                                                                >  2:290474/1‑152 (MQ=255)
                                                                                                                                                                                                          gACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCGACTGGTGCTGCATTCGAGAACCTTGGGCTGGGCGGtgctg                                                                                                                                                                <  2:46135/85‑1 (MQ=255)
                                                                                                                                                                                                          gACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                           <  2:715140/170‑1 (MQ=255)
                                                                                                                                                                                                             gCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCTCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTg                                                           >  2:632224/1‑183 (MQ=255)
                                                                                                                                                                                                               tGGGCAGATCAGCTGGGGCAGCAAGTCGGCGGGGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCACGTCACCGACACCCTCGACGAAGTGGTGGCCACGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTga                                                                                 <  2:46757/159‑1 (MQ=14)
                                                                                                                                                                                                                 gCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                           <  2:588651/195‑1 (MQ=255)
                                                                                                                                                                                                                 gCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                           >  1:588654/1‑195 (MQ=255)
                                                                                                                                                                                                                  ccAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTaca                                                                      >  2:827808/1‑167 (MQ=255)
                                                                                                                                                                                                                   cAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGAcc  <  2:1098176/234‑1 (MQ=40)
                                                                                                                                                                                                                   cAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGac                                                                                                                                            <  1:481755/96‑1 (MQ=255)
                                                                                                                                                                                                                   cAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTCGAGCTGGGCGGTGCTGCTTCGGTGCAAGTCACCGac                                                                                                                                            >  2:481752/1‑96 (MQ=255)
                                                                                                                                                                                                                     gATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCAcc                                                                                            <  2:882742/142‑1 (MQ=14)
                                                                                                                                                                                                                      aTCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGa                     <  1:90035/212‑1 (MQ=40)
                                                                                                                                                                                                                         aGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTga                                                                                 >  2:1163476/1‑149 (MQ=255)
                                                                                                                                                                                                                         aGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTg                                                                                  <  2:86130/147‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                
AGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACC  >  AE015451/199330‑199774

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