Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 202,111 G→A 9.5% E2539E (GAG→GAA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451202,1110GA 9.5% 154.6 / 4.9 52E2539E (GAG→GAAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base A (2/3);  total (24/28)
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.03e-01

GAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCT  >  AE015451/201887‑202344
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
gAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTgcg                                                                                                                                                                                                                                  <  2:354523/234‑1 (MQ=255)
                    gACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGAcc                                                                                                                                                                                                                                        <  1:483807/208‑1 (MQ=255)
                        ccGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGAcgc                                                                                                                                                                                                                 <  1:684338/227‑1 (MQ=255)
                           gTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGAt                                                                                                                                                                                                       <  1:271107/234‑1 (MQ=255)
                             cAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGt                                                                                                                                                                                                                                     <  1:315130/202‑1 (MQ=255)
                                 aCTACCGTCACTGACGATCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTa                                                                                                                                                                                               <  1:395596/236‑1 (MQ=255)
                                          aCTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGc                                                                                                                                                                                      <  1:472887/236‑1 (MQ=255)
                                             gACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTgg                                                                                                                                                                                   >  1:300692/1‑236 (MQ=40)
                                               cGAGCCGGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGAcc                                                                                                                                                                                                                                        <  2:346278/181‑1 (MQ=255)
                                                  gCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTgc                                                                                                                                                                                                                                   >  1:397684/1‑183 (MQ=40)
                                                  gCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGAc                                                                                                                                                                                                    <  1:614400/214‑1 (MQ=255)
                                                  gCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGAc                                                                                                                                                                                                    >  2:614399/1‑214 (MQ=40)
                                                  gCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGcc                                                                                                                                                                                 >  1:309365/1‑233 (MQ=38)
                                                  gCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCa                                                                                                                                                                                <  2:300691/234‑1 (MQ=255)
                                                   ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCa                                                                                                                                                                                <  2:636526/233‑1 (MQ=255)
                                                   ccAGGCACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAg                                                                                                                                                                                                            <  2:6398/205‑1 (MQ=255)
                                                           tCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAgg                                                                                                                                                                                          <  1:363159/215‑1 (MQ=255)
                                                           tCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTg                                                                                                                                                                      <  1:477392/235‑1 (MQ=255)
                                                               ggCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCaa                                                                                                                                                                                                             >  2:554412/1‑192 (MQ=40)
                                                                gCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGCCTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGAcgc                                                                                                                                                                                                                 <  1:149452/187‑1 (MQ=255)
                                                                   aCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACACACGACGCGCaa                                                                                                                                                                                                             >  1:157488/1‑188 (MQ=255)
                                                                        ggCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                         >  1:734329/1‑235 (MQ=255)
                                                                        ggCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACACGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCa                                                                                                                                                                                                                       >  1:451981/1‑173 (MQ=255)
                                                                            aCATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGAGGACGTCTGCCAGGAcc                                                                                                                                                                                       >  2:280505/1‑200 (MQ=255)
                                                                                  tCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGAt                                                                                                                                                 <  2:734328/233‑1 (MQ=255)
                                                                                                    ccGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAgg                                                                                                                                                                                          >  1:265172/1‑174 (MQ=255)
                                                                                                    ccGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAgg                                                                                                                                                                                          <  2:265171/174‑1 (MQ=255)
                                                                                                     gaaCCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCAc                                                                                                                                                              <  2:378252/199‑1 (MQ=255)
                                                                                                         cAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                         >  2:362351/1‑202 (MQ=255)
                                                                                                         cAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCt                                                                                                                       >  1:639397/1‑236 (MQ=255)
                                                                                                                     gCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGac                                                                                                                                                 >  1:445824/1‑197 (MQ=255)
                                                                                                                          gaACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACg                                                                                                      <  1:345679/236‑1 (MQ=255)
                                                                                                                           aaCGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCCTCAAGGCCGGTGAGACCAg                                                                                                                                                                                                                                      >  2:1176/1‑107 (MQ=255)
                                                                                                                                               cGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACg                                                                                                      >  2:266282/1‑215 (MQ=255)
                                                                                                                                                                  cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTGCACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                         <  1:107289/145‑1 (MQ=255)
                                                                                                                                                                       tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAg                                                          >  1:489721/1‑235 (MQ=255)
                                                                                                                                                                                   cGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCGGTGGAt                                                                                                                                                 >  1:327557/1‑136 (MQ=255)
                                                                                                                                                                                    gTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGtt                                                                                              >  2:177958/1‑186 (MQ=255)
                                                                                                                                                                                      gACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCt                                           <  1:700180/235‑1 (MQ=255)
                                                                                                                                                                                       aCCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTcg                                                                           <  1:602558/202‑1 (MQ=255)
                                                                                                                                                                                              gCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                         >  2:151970/1‑117 (MQ=255)
                                                                                                                                                                                               caacaaCGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGCCGCGCAAGGCGATGACTTCTCCCAAGatg                                                                                                                                                                                       >  2:64778/1‑84 (MQ=14)
                                                                                                                                                                                               caacaaCGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGTTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGTGGCAACGCTTCGGTTCAAGTGACCGACACCCTcg                                                                           <  2:156714/194‑1 (MQ=255)
                                                                                                                                                                                               caacaaCGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGTTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGGGGCAACGCTTCGGTTCAAGTGACCGACACCCTcg                                                                           >  1:156714/1‑194 (MQ=255)
                                                                                                                                                                                                aacaacGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                          <  2:155661/114‑1 (MQ=19)
                                                                                                                                                                                                aacaacGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCAAGTcc                                                                                                                                                          >  1:155661/1‑114 (MQ=2)
                                                                                                                                                                                                     cGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAg                              <  2:323016/233‑1 (MQ=255)
                                                                                                                                                                                                         cAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGCCGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGtt                                                                                              <  1:151970/165‑1 (MQ=255)
                                                                                                                                                                                                          aGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACGCCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAg                                                                      >  1:254682/1‑188 (MQ=255)
                                                                                                                                                                                                               aCCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATc                   <  2:325109/234‑1 (MQ=255)
                                                                                                                                                                                                                ccATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCAc                 <  1:554413/235‑1 (MQ=255)
                                                                                                                                                                                                                         ccGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCAt        >  1:202044/1‑235 (MQ=255)
                                                                                                                                                                                                                              gagaCCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCt  >  1:320657/1‑236 (MQ=255)
                                                                                                                                                                                                                              gagaCCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCt  >  1:521997/1‑236 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
GAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCT  >  AE015451/201887‑202344

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