Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 200,386 C→T 6.5% D1964D (GAC→GAT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451200,3860CT 6.5% 231.4 / 3.9 75D1964D (GAC→GATPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base T (3/2);  total (40/35)
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.70e-01

CACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGT  >  AE015451/200155‑200613
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   
cacCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                                                                                                 <  1:779197/236‑1 (MQ=37)
           cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaa                                                                                                                                                                                                                       >  1:594235/1‑235 (MQ=40)
                 gCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACg                                                                                                                                                                                                                                    <  1:209249/216‑1 (MQ=39)
                 gCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCAc                                                                                                                                                                                                                <  1:605921/236‑1 (MQ=255)
                      gACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCg                                                                                                                                                                                                                            <  1:350084/219‑1 (MQ=39)
                            tATCCGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCt                                                                                                                                                                                                                                <  1:61396/209‑1 (MQ=255)
                                aGGACGCTGGCCAGATCAGCCTGGGCATCAGCTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTc                                                                                                                                                                                                 <  1:1083032/236‑1 (MQ=255)
                                       tGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaacaa                                                                                                                                                                                                                    >  2:725069/1‑210 (MQ=255)
                                         gCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaa                                                                                                                                                                                                                       >  2:896277/1‑205 (MQ=40)
                                         gCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaa                                                                                                                                                                                                                       <  1:896280/205‑1 (MQ=255)
                                          ccAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGAcc                                                                                                                                                                                        >  1:78004/1‑235 (MQ=40)
                                            aGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACg                                                                                                                                                                                                                                    <  2:866430/189‑1 (MQ=255)
                                            aGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGaa                                                                                                                                                                                                                          >  2:1083029/1‑198 (MQ=255)
                                               tCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACgg                                                                                                                                                                                   <  1:982254/235‑1 (MQ=40)
                                                      gggCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                                                                                  >  2:801388/1‑197 (MQ=14)
                                                      gggCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                                                                                  <  1:801391/197‑1 (MQ=14)
                                                          aTCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgtcgt                                                                                                                                                                       >  1:367797/1‑236 (MQ=255)
                                                                tCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCt                                                                                                                                                                                                        >  1:229785/1‑197 (MQ=40)
                                                                     aGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGt                                                                                                                                                                                                      <  1:935576/194‑1 (MQ=39)
                                                                     aGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGcc                                                                                                                                                             >  1:635535/1‑235 (MQ=40)
                                                                                   gTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                                                                  <  1:688126/168‑1 (MQ=39)
                                                                                           ttCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTCACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACTCGGAGc                                                                                                                                                                                                         >  2:3669/1‑168 (MQ=255)
                                                                                              gagaACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAg                                                                                                                                                                                              >  2:775181/1‑177 (MQ=35)
                                                                                                gaACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGtt                                                                                                                                             >  2:372765/1‑224 (MQ=39)
                                                                                                       gAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCAcc                                                                                                                          <  1:906162/236‑1 (MQ=255)
                                                                                                          cTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCa                                                                                                                         <  2:635532/234‑1 (MQ=14)
                                                                                                          cTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCAcc                                                                                                                          <  2:426070/233‑1 (MQ=25)
                                                                                                           tGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                                                                         <  2:594232/233‑1 (MQ=255)
                                                                                                             gcactgtggttgttgatcggcagaccatctttgttggtcagggtgatggtgtaggtgatctcaccaCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGCTGGTCTGCCGATCAACAACCAc                                                                                                                                                                                                                <  1:294872‑M1/78‑1 (MQ=255)
                                                                                                                    ggttattgatcggcaagccgtctttgttggtcagggtaatggtgtaagtgatctcaccgCCTTCGGTGACCGAAGGCGGTGAGATCACTTATACCATTACCCTGACCAACAAAGACg                                                                                                                                                                                                                                    <  2:912280‑M1/58‑1 (MQ=255)
                                                                                                                     cGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTCCTTCAAGCTGACCGACGGCACCACtgtcgt                                                                                                                                                                       >  2:174132/1‑177 (MQ=2)
                                                                                                                     cGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATg                                                                                                            >  1:927612/1‑236 (MQ=37)
                                                                                                                       cTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGt                                                                                                                                                                                                      <  2:888745/144‑1 (MQ=32)
                                                                                                                         tCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGt                                                                                                           <  2:407881/233‑1 (MQ=38)
                                                                                                                               cAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGCCTAt                                                                                                             >  1:107323/1‑225 (MQ=37)
                                                                                                                                     aCCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACgcgc                                                                                             <  1:658737/235‑1 (MQ=38)
                                                                                                                                          cacCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGt                                                                                                                                              >  2:509187/1‑181 (MQ=35)
                                                                                                                                          cacCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGt                                                                                                                                              <  1:509190/181‑1 (MQ=37)
                                                                                                                                            ggttggtcagggtaatggtgtaagtgatctcaccgCCTTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACTACTGTCGTCGTGGCAGCCAACAGca                                                                                                                                                      >  1:451587‑M1/36‑171 (MQ=255)
                                                                                                                                             ccTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCggt                                                                                         <  2:927609/231‑1 (MQ=38)
                                                                                                                                                cgacgaAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGCTCGGCCACTGCAACCGCACCAg                                                                                                                        >  2:1067893/1‑197 (MQ=14)
                                                                                                                                                cgacgaAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACGGGtt                                                                                                                                             >  2:798824/1‑176 (MQ=17)
                                                                                                                                                     aaGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCAcc                                                                                                                          <  2:78002/190‑1 (MQ=39)
                                                                                                                                                          gtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGc                                                                                >  2:820581/1‑227 (MQ=39)
                                                                                                                                                          gtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCCATGTCGGCACCAACGCGCCg                                                                                           >  2:532125/1‑216 (MQ=39)
                                                                                                                                                          gtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCAGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGCTCTGCCCATCAACAACCACTc                                                                                                                                                                                                              >  2:64252/1‑101 (MQ=255)
                                                                                                                                                            ggCCAAGCTCACTCCCACCCCTTCGGTGACCGAAGGCGGTGAGCTCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCg                                                                           <  1:124266/230‑1 (MQ=255)
                                                                                                                                                                 aggtgatctcaccgCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGCCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATc                                                                            <  1:454162‑M1/210‑1 (MQ=255)
                                                                                                                                                                            aCCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGt                                                                    <  2:405995/221‑1 (MQ=25)
                                                                                                                                                                               CCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGAt                                                                                                                                                                                                                          >  2:1042113‑M1/1‑68 (MQ=255)
                                                                                                                                                                               CCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGAt                                                                                                                                                                                                                          <  1:1042116‑M1/68‑1 (MQ=255)
                                                                                                                                                                                 TTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATg                                                                                 >  1:728956‑M1/1‑203 (MQ=255)
                                                                                                                                                                                   CGGTGACCGAAGGCGGTGAGATCACCTACAGCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGa                                                                                                                                                                                                           >  2:31616‑M1/1‑79 (MQ=255)
                                                                                                                                                                                          cGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGaaa                                       <  1:820584/236‑1 (MQ=34)
                                                                                                                                                                                              ggCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCa                                                                                                   >  2:728953/1‑172 (MQ=33)
                                                                                                                                                                                                                  cATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTgg                           >  1:115582/1‑224 (MQ=35)
                                                                                                                                                                                                                  cATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGGTCGAAAACCTGAACCTGGACAAGAccc                  >  2:1010668/1‑233 (MQ=37)
                                                                                                                                                                                                                   aTCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctgaa                                <  2:1015123/218‑1 (MQ=35)
                                                                                                                                                                                                                   aTCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctgaa                                >  1:1015126/1‑218 (MQ=35)
                                                                                                                                                                                                                   aTCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCAGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGt              >  1:871104/1‑236 (MQ=18)
                                                                                                                                                                                                                    tCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCg                                                                                           >  2:15696/1‑158 (MQ=33)
                                                                                                                                                                                                                    tCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                 <  2:372494/232‑1 (MQ=35)
                                                                                                                                                                                                                    tCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                 >  1:372497/1‑232 (MQ=35)
                                                                                                                                                                                                                      tcccTGACCAACAAAGACGGTCTGCCGATCAACAACCTCTCGGAGCTGTACTTCCAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAGCAGCACCACCGGTTCGGACACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAAcgcg                                                                                              <  1:107471/152‑1 (MQ=255)
                                                                                                                                                                                                                       cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAAcgc                                                                                               >  2:749599/1‑151 (MQ=33)
                                                                                                                                                                                                                       cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTc             >  2:662909/1‑233 (MQ=35)
                                                                                                                                                                                                                        ccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa         >  1:1052065/1‑236 (MQ=37)
                                                                                                                                                                                                                             ccAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCActaat                                                                                                 <  2:578206/143‑5 (MQ=12)
                                                                                                                                                                                                                             ccAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCgg                                                                     >  2:470522/1‑171 (MQ=35)
                                                                                                                                                                                                                               aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCAGCAGACAACGTCTATGTCGGCACCAACGCGCCggtggt                                                                                      >  2:599517/1‑152 (MQ=18)
                                                                                                                                                                                                                               aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCggtggt                                                                                      <  1:599520/152‑1 (MQ=33)
                                                                                                                                                                                                                               aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAg    >  1:608757/1‑234 (MQ=35)
                                                                                                                                                                                                                               aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGt  >  1:101097/1‑236 (MQ=35)
                                                                                                                                                                                                                               aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGt  >  1:678563/1‑236 (MQ=37)
                                                                                                                                                                                                                                   aaGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGt                                                                                                           <  1:551763/127‑1 (MQ=35)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   
CACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGT  >  AE015451/200155‑200613

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