Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,915 T→C 9.8% F1807F (TTT→TTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451199,9150TC 9.8% 222.5 / 8.4 70F1807F (TTT→TTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (2/5);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

TACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCG  >  AE015451/199703‑200147
                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
tACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAaacctgaacctg                                                                                                                                                                                                                           >  1:625585/1‑228 (MQ=255)
     cATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAg                                                                                                                                                                                                                     >  2:575358/1‑229 (MQ=255)
     cATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAg                                                                                                                                                                                                                     <  1:575361/229‑1 (MQ=255)
           ccTGACCACCAAGGACGGTCTGCCGATCAAGAACCACGCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAaacctg                                                                                                                                                                                                                                 <  2:163688/211‑1 (MQ=255)
                  aacaaAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGaaaa                                                                                                                                                                                                                                     >  2:319498/1‑200 (MQ=255)
                  aacaaAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGt                                                                                                                                                                                                 >  1:446460/1‑236 (MQ=255)
                    caaAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCa                                                                                                                                                                                               <  1:929230/236‑1 (MQ=255)
                        gACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACcg                                                                                                                                                                                            >  1:518560/1‑235 (MQ=255)
                              cTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAg                                                                                                                                                                                                          >  2:249535/1‑215 (MQ=40)
                               tGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                                                        <  1:735775/216‑1 (MQ=38)
                                  cGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAgg                                                                                                                                                                                  >  1:421748/1‑235 (MQ=255)
                                           aCCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGAcccc                                                                                                                                                                                                                <  1:344702/196‑1 (MQ=40)
                                             cgcTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAaacctg                                                                                                                                                                                                                                 <  1:684309/175‑1 (MQ=255)
                                                  ggAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGa                                                                                                                                                                   >  1:9692/1‑234 (MQ=255)
                                                  ggAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAg                                                                                                                                                                 >  1:921079/1‑236 (MQ=255)
                                                  ggAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGCAACGaa                                                                                                                                                                  >  1:537307/1‑235 (MQ=9)
                                                      cTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCg                                                                                                                                                              <  1:249538/235‑1 (MQ=255)
                                                      cTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCg                                                                                                                                                              <  1:593576/235‑1 (MQ=255)
                                                            ttCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGc                                                                                                                                                               >  1:575098/1‑228 (MQ=255)
                                                             tCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGt                                                                                                                                                       >  1:34829/1‑235 (MQ=255)
                                                                 gCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAg                                                                                                                                                                                   >  1:707275/1‑203 (MQ=255)
                                                                 gCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAgg                                                                                                                                                  >  1:580904/1‑236 (MQ=255)
                                                                   tGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAg                                                                                                                                                                                                   >  1:778833/1‑185 (MQ=255)
                                                                        gACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGc                                                                                                                                                                                                         <  1:238650/174‑1 (MQ=255)
                                                                        gACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGc                                                                                                                                                                                                         >  2:238648/1‑174 (MQ=255)
                                                                              accacTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCAca                                                                                                                                                                              >  2:655042/1‑194 (MQ=255)
                                                                                             gCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACACCGGTCAGCACCGAg                                                                                                                                                                                                   >  2:136671/1‑159 (MQ=255)
                                                                                                gCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTgg                                                                                                                                                                                                                          >  1:519795/1‑133 (MQ=255)
                                                                                                gCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTgg                                                                                                                                                                                                                          <  2:519792/133‑1 (MQ=255)
                                                                                                 ccAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGc                                                                                                                   <  1:459304/235‑1 (MQ=255)
                                                                                                  cAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCa                                                                                                                                                                                    >  1:629760/1‑169 (MQ=255)
                                                                                                  cAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGcc                                                                                                                  >  1:56110/1‑235 (MQ=255)
                                                                                                           cacGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTc                                                                                                         <  1:544790/235‑1 (MQ=255)
                                                                                                           cacGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTc                                                                                                         <  1:232768/235‑1 (MQ=255)
                                                                                                             cGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCaa                                                                                                       <  1:699996/235‑1 (MQ=255)
                                                                                                              cggttCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCaaa                                                                                                      >  1:400211/5‑235 (MQ=9)
                                                                                                               ggCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCaa                                                                                                       <  2:421745/233‑1 (MQ=255)
                                                                                                                gCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGCCACCGACGAGCCAGGCCCACCAGGCAACGAAg                                                                                                                                                                 >  2:204841/1‑174 (MQ=255)
                                                                                                                 cTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGa                                                                                                  <  1:169202/236‑1 (MQ=255)
                                                                                                                   cGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGc                                                                                                                                                                                 >  1:358833/1‑155 (MQ=21)
                                                                                                                          tGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCa                                                                                            <  2:346786/233‑1 (MQ=255)
                                                                                                                           gCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCa                                                                                            <  2:518557/232‑1 (MQ=255)
                                                                                                                            cAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCAc                                                                                           >  1:314317/1‑232 (MQ=255)
                                                                                                                             aaCCGCCCCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGc                                                                                                                                                                                      <  2:737112/140‑1 (MQ=38)
                                                                                                                                cGCCCCAGACAACGTCTATGTCGGCACCAACGAGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGc                                                                                                                                                                                 <  2:274735/142‑1 (MQ=255)
                                                                                                                                   aCCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGAc                                                                                                 <  1:156234/219‑1 (MQ=255)
                                                                                                                                       gACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAAc                                                                              >  2:741252/1‑234 (MQ=255)
                                                                                                                                               cTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGAc                                                                                                                                                            >  1:659374/1‑148 (MQ=255)
                                                                                                                                                         aCCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACg                                                                                                                                       >  2:628279/1‑159 (MQ=255)
                                                                                                                                                          ccAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCg                                                                                                                        >  2:156233/1‑173 (MQ=255)
                                                                                                                                                          ccAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCa                                                                                            >  1:836646/1‑201 (MQ=255)
                                                                                                                                                             aCGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGAcc                                                                                                <  2:34829/194‑1 (MQ=255)
                                                                                                                                                                              gCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACgcc                                                            >  2:219189/1‑211 (MQ=38)
                                                                                                                                                                               ccAGCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAAACTGGACAAGGCCCCGGTCAGCGCCGAGGTGACCGACGAGCCAGGCACTCCAGGCAACGAAGGCg                                                                                                                                                              <  2:179019/114‑1 (MQ=11)
                                                                                                                                                                                cATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGAACCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACGCCAGGCa                                                                                                                                                                       >  2:647852/1‑104 (MQ=255)
                                                                                                                                                                                  tCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGGCAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                           <  1:746968/226‑1 (MQ=37)
                                                                                                                                                                                   cGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCAc                                                                                           >  1:820431/1‑177 (MQ=39)
                                                                                                                                                                                   cGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCAc                                                                                           <  2:820428/177‑1 (MQ=40)
                                                                                                                                                                                    gccgcgGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTgg                                                                                                                                                                                                                          <  2:893894/47‑1 (MQ=255)
                                                                                                                                                                                    gACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAcc                                                <  2:575095/219‑1 (MQ=39)
                                                                                                                                                                                     aCGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAc                                                               <  2:675211/203‑1 (MQ=14)
                                                                                                                                                                                                ggcgCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCa                    <  1:874931/235‑1 (MQ=38)
                                                                                                                                                                                                       aTGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCTAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAc                                                               >  2:341666/1‑185 (MQ=39)
                                                                                                                                                                                                       aTGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCTAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAc                                                               <  1:341669/185‑1 (MQ=255)
                                                                                                                                                                                                       aTGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTACACGTCAACCAGcc                                                                         <  2:437017/175‑1 (MQ=9)
                                                                                                                                                                                                          cGTGGAAGTTTGCAACCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCAc                                                                                           <  1:548541/154‑1 (MQ=255)
                                                                                                                                                                                                          cGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGAACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACGGTGCACGTCAAGCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAAGTCTCAGGTCACCATCAAGGCTGGCg           >  2:883679/1‑234 (MQ=255)
                                                                                                                                                                                                             ggAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCa                               <  2:56110/211‑1 (MQ=39)
                                                                                                                                                                                                                aGTTCGAACACCTGAATCTGGACAAGACCCCGGTCAGCACGGAGGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAAGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCa                    <  2:823380/219‑1 (MQ=11)
                                                                                                                                                                                                                 gTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg  <  1:337891/236‑1 (MQ=255)
                                                                                                                                                                                                                   ttGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg  <  2:580901/234‑1 (MQ=255)
                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
TACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCG  >  AE015451/199703‑200147

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