Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,744 T→C 13.7% N1750N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451199,7440TC13.7% 200.0 / 20.2 73N1750N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (4/6);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGG  >  AE015451/199509‑199971
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAAt                                                                                                                                                                                                                                     >  1:954896/1‑236 (MQ=255)
       cGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCAc                                                                                                                                                                                                                               <  1:557631/235‑1 (MQ=255)
               tCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTg                                                                                                                                                                                                                      <  1:894247/236‑1 (MQ=255)
                    cAGGATGCGGGCCAGATCAGCCTGGGCATCAGCTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTACGGTTGAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAAcc                                                                                                                                                                                                                                 <  1:829052/220‑1 (MQ=255)
                    cAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTTCTTCGGTTCGAGTCACCGACACCCTCGATGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACt                                                                                                                                                                                                                  >  1:616173/1‑235 (MQ=255)
                    cAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATACCCACTCGGAGCTGTACtt                                                                                                                                                                                                                 >  1:33950/1‑236 (MQ=255)
                              gCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGa                                                                                                                                                                                                                          >  2:757018/1‑217 (MQ=38)
                              gCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGa                                                                                                                                                                                                                          <  1:757019/217‑1 (MQ=38)
                                cAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCCGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGa                                                                                                                                                                                                     <  1:445704/236‑1 (MQ=18)
                                  gATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCGTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACg                                                                                                                                                                                                   <  1:281011/236‑1 (MQ=37)
                                           gggCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGcacca                                                                                                                                                                                             >  2:517344/1‑233 (MQ=37)
                                                  aaCTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTgtcg                                                                                                                                                                                       <  2:732290/232‑1 (MQ=255)
                                                   aCTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTg                                                                                                                                                                                  <  1:986045/236‑1 (MQ=255)
                                                   aCTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTg                                                                                                                                                                                  <  1:517345/236‑1 (MQ=18)
                                                     tCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTg                                                                                                                                                                                  <  1:424366/234‑1 (MQ=40)
                                                          ggTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTcgtcgt                                                                                                                                                                                   >  1:853269/1‑228 (MQ=40)
                                                               aCGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGAc                                                                                                                                                                                                        >  1:105181/1‑202 (MQ=35)
                                                                                ttCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGc                                                                                                                                                                                                            <  2:509552/181‑1 (MQ=14)
                                                                                ttCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGc                                                                                                                                                                                                            >  1:509553/1‑181 (MQ=14)
                                                                                ttCGAGAACCTTGAGCTGGGCGGTGCTGCGTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCCGACCCACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTTCTTCAAGc                                                                                                                                                                                                            >  2:698115/1‑181 (MQ=255)
                                                                                 tCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACt                                                                                                                                                    >  1:199728/1‑236 (MQ=40)
                                                                                           tGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCgtcg                                                                                                                                                                                       >  1:1130812/1‑191 (MQ=18)
                                                                                              gCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCg                                                                                                                                                                                                                            >  1:973657/1‑151 (MQ=14)
                                                                                              gCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCg                                                                                                                                                                                                                            <  2:973656/151‑1 (MQ=14)
                                                                                               cTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAAcc                                                                                                                                              <  1:76176/228‑1 (MQ=40)
                                                                                               cTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAg                                                                                                                                       <  1:550109/235‑1 (MQ=255)
                                                                                                     ggTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGt                                                                                                                                >  1:967786/1‑236 (MQ=40)
                                                                                                                  tgcAAGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGcacca                                                                                                                                                                                             <  1:509341/160‑1 (MQ=9)
                                                                                                                   tCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGAc                                                                                                                                                                                                        <  1:652379/150‑1 (MQ=32)
                                                                                                                    cAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGca                                                                                                                                                                                                >  2:809086/1‑157 (MQ=31)
                                                                                                                    cAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGca                                                                                                                                                                                                <  1:809087/157‑1 (MQ=31)
                                                                                                                    cAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGAc                                                                                                                                                     >  2:460578/1‑200 (MQ=39)
                                                                                                                     aaGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCActa                                                                                                                  <  1:679600/234‑3 (MQ=14)
                                                                                                                      cgTCACCGCCACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCg                                                                                                                                             <  1:448941/205‑1 (MQ=255)
                                                                                                                         cACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAAcc                                                                                                                                              >  1:439022/1‑202 (MQ=38)
                                                                                                                               cacCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCggcg                                                                                                       >  1:329257/1‑233 (MQ=37)
                                                                                                                                           ggtggtgGCCAAGCTGACTGCCACCCCGTCGGTAACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCa                                                                                                                  <  1:853117/211‑1 (MQ=255)
                                                                                                                                              tgtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTCGCAGCCAACAGCACCACCGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAAcgcg                                                                                                             <  1:822117/213‑1 (MQ=21)
                                                                                                                                                  gCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGt                                                                                   <  1:812640/236‑1 (MQ=39)
                                                                                                                                                     aaGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGccc                                                                                                          >  2:548504/1‑209 (MQ=38)
                                                                                                                                                             tGCCACCCCGTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGa                                                                                                                                      <  2:32387/174‑1 (MQ=255)
                                                                                                                                                                    ccTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCaa                                                                                                                                                <  2:792108/157‑1 (MQ=255)
                                                                                                                                                                    ccTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCaa                                                                                                                                                >  1:792109/1‑157 (MQ=255)
                                                                                                                                                                     cTTCGGTGACCGAAGGCGGTGCGATCACCTACACCATCACCCTGACCAACAACGACGGTCTGCCGATCAATAACCACTCGGCGCTGTACTTCAAGCTGACCGACGGCACCACTGTCCTCGTGGCAGCCAACAGCACCACg                                                                                                                                                                >  2:200352/1‑140 (MQ=255)
                                                                                                                                                                        cGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAAcc                                                                                                                                              <  2:785775/155‑1 (MQ=40)
                                                                                                                                                                               cGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGTAGCTGTACTTCAAGCTGTCCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGCCTATGTCGGCACCAAcgcc                                                                                                             >  1:1070736/1‑180 (MQ=255)
                                                                                                                                                                               cGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGaaa                                                      <  1:121303/236‑1 (MQ=40)
                                                                                                                                                                                        cgagaTCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATc                                                                                           <  2:1029506/189‑1 (MQ=255)
                                                                                                                                                                                          agaTCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAaacctgaacctg                                           >  1:706253/1‑236 (MQ=255)
                                                                                                                                                                                               aCCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTg                                                                 <  1:385265/209‑1 (MQ=40)
                                                                                                                                                                                                    cacCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGAc                                   <  2:695333/234‑1 (MQ=37)
                                                                                                                                                                                                            cgcTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGCGCGGCGACTGCAGCCGCCCCAGACAACGTCt                                                                                                                              >  2:188886/3‑135 (MQ=11)
                                                                                                                                                                                                            cccTGACTAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACAGGCTCGGCGACTGCAACCGCACCAGACAACGTCTacgg                                                                                                                          >  2:253044/1‑136 (MQ=255)
                                                                                                                                                                                                            cccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTg                                                                                                                                                                                  <  2:255109/83‑1 (MQ=38)
                                                                                                                                                                                                            cccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTg                                                                                                                                                                                  >  1:255110/1‑83 (MQ=38)
                                                                                                                                                                                                            cccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGAt                                                                                                     >  2:502321/1‑160 (MQ=255)
                                                                                                                                                                                                            cccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGAt                                                                      >  2:894392/1‑191 (MQ=255)
                                                                                                                                                                                                            cccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATg                                                                     >  1:735233/1‑192 (MQ=255)
                                                                                                                                                                                                            cccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCAACACt                                                                                                                                                                                           >  1:241141/1‑74 (MQ=25)
                                                                                                                                                                                                             ccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCgg                                                                  >  1:310079/1‑193 (MQ=40)
                                                                                                                                                                                                                  ccAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCACTGCCATCGACGCGGTCAGCGGCGCa                                                                         >  1:1055727/1‑182 (MQ=37)
                                                                                                                                                                                                                  ccAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGt                             >  2:527526/1‑226 (MQ=255)
                                                                                                                                                                                                                   caacaaAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATg                                                                                                <  2:380786/158‑1 (MQ=255)
                                                                                                                                                                                                                    aacaaAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAg                   >  2:492044/1‑234 (MQ=37)
                                                                                                                                                                                                                    aacaaAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACCGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCgg                                                                                                                                                         >  2:496004/1‑100 (MQ=37)
                                                                                                                                                                                                                        aaGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAACCCTGAACCTGGACAAGACCCCGGt                             <  2:616172/220‑1 (MQ=255)
                                                                                                                                                                                                                              gTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCg                                                                                                          >  2:767647/1‑137 (MQ=255)
                                                                                                                                                                                                                               tCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGc      <  1:1106086/236‑1 (MQ=255)
                                                                                                                                                                                                                               tCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGc      <  1:527527/236‑1 (MQ=255)
                                                                                                                                                                                                                                 tGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCa    >  2:686243/1‑236 (MQ=255)
                                                                                                                                                                                                                                  gCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGcc     >  1:1090048/1‑234 (MQ=255)
                                                                                                                                                                                                                                   ccGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGt                             <  2:33950/209‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGcc     <  2:967785/233‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAgg  >  1:375672/1‑236 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
CGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGG  >  AE015451/199509‑199971

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: