Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 200,473 T→C 7.7% G1993G (GGT→GGC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451200,4730TC 7.7% 242.1 / 2.4 79G1993G (GGT→GGCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (41/31);  new base C (3/4);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACC  >  AE015451/200247‑200706
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         
tCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCa                                                                                                                                                                                                                                    >  2:651576/1‑234 (MQ=40)
       ccTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCg                                                                                                                                                                                                                           <  1:491040/236‑1 (MQ=255)
        cTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAAc                                                                                                                                                                                                                             <  2:456578/233‑1 (MQ=255)
             gCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAg                                                                                                                                                                                                                     >  1:951801/1‑236 (MQ=255)
               tGGGCGGTTCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                                                                                                                                                                      <  2:965007/233‑1 (MQ=38)
               tGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                                                                                                                                                                      <  2:732665/233‑1 (MQ=255)
                  gCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAAc                                                                                                                                                                                                                <  1:815420/236‑1 (MQ=255)
                      ttCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGCTCGGCCACTGCa                                                                                                                                                                                                                               >  2:894553/1‑217 (MQ=37)
                                         aCCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAg                                                                                                                                                                                                                     >  2:1134310/1‑208 (MQ=38)
                                                           ctggtggCCAAGCTGACTGCCACCCCTGCGGTGACCGAAGGCGGTGAGATCACCTACCCCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCAc                                                                                                                                                                                                                        <  1:253453/186‑1 (MQ=255)
                                                             ggtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCAcc                                                                                                                                                                                                                       <  2:1125600/186‑1 (MQ=14)
                                                             ggtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCAcc                                                                                                                                                                                                                       >  1:1125601/1‑186 (MQ=18)
                                                                  ccAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCaa                                                                                                                                                                                               >  2:447647/1‑205 (MQ=38)
                                                                  ccAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTc                                                                                                                                                                <  1:866197/236‑1 (MQ=39)
                                                                    aaGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCa                                                                                                                                                               >  1:387691/1‑235 (MQ=39)
                                                                                  gccTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCAt                                                                                                                                                                          <  1:20383/209‑1 (MQ=21)
                                                                                    cTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTg                                                                                                                                               <  1:651577/235‑1 (MQ=37)
                                                                                         gTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGtt                                                                                                                                         <  1:533535/236‑1 (MQ=37)
                                                                                               gAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGCTGCGTGGAAGTTCGaaaa                                                                                                                                   >  1:311251/1‑236 (MQ=18)
                                                                                               gAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACCGCAACCGCACCAGACAACGTCTATGTCGGCACCAAcg                                                                                                                                                                                             >  2:1097864/1‑178 (MQ=18)
                                                                                                          gaTCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctg                                                                                                                               >  2:966791/1‑229 (MQ=34)
                                                                                                                 tACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAg                                                                                                                   <  2:856398/234‑1 (MQ=37)
                                                                                                                      cATCACCCTGACCACCAAAGCCGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGt                                                                                                                                                                 <  2:310078/183‑1 (MQ=25)
                                                                                                                        tCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCg                                                                                                                                                                        >  2:768837/1‑174 (MQ=34)
                                                                                                                         cccccTGACCAACAAAGACGGCCTGCCGATCAACACCCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCa                                                                                                                                                                           <  2:281822/168‑1 (MQ=18)
                                                                                                                           cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCgg                                                                                                            >  2:594885/1‑231 (MQ=37)
                                                                                                                           cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGCGGCCGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCATCGCGCCGGTGGTCAATGCCATCGACGCGGTCAgc                                                                                                                                                             >  1:366391/1‑182 (MQ=25)
                                                                                                                                   agcaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCa                                                                                                         <  1:841384/224‑1 (MQ=14)
                                                                                                                                   aacaacGACGGTCTGCCGATCAACAACCACGCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                                                                                                                                                                      <  2:94309/117‑1 (MQ=16)
                                                                                                                                   aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGAc                                                                                                                                                                                                                   >  2:545942/1‑120 (MQ=35)
                                                                                                                                   aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGt                                                                                                                                                >  2:79067/1‑187 (MQ=35)
                                                                                                                                   aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAg                                                                                                 >  2:756561/1‑234 (MQ=35)
                                                                                                                                   aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCACGGCCAGCACCGAg                                                                                                 >  2:1099882/1‑234 (MQ=37)
                                                                                                                                   aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCAcaa                                                                                                                                                                                                >  1:492880/1‑137 (MQ=21)
                                                                                                                                        aGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAACCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAg                                                                                                                                                              >  1:999001/1‑168 (MQ=25)
                                                                                                                                          aCGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACcgac                                                                                        <  1:452687/236‑1 (MQ=40)
                                                                                                                                             gTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGcc                                                                                                                                                                                         <  1:350224/136‑1 (MQ=34)
                                                                                                                                             gTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACcgacga                                                                                      >  1:1025418/1‑235 (MQ=40)
                                                                                                                                                 gCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGAccctgtctctta                                                                                 >  1:242811/1‑226 (MQ=255)
                                                                                                                                                  ccGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCa                                                                                                                                                       >  2:488123/1‑165 (MQ=35)
                                                                                                                                                   cGATCAACAACCACTCGTAGCTGTACTTCAAGCTGACCGCCGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCa                                                                                                         <  2:242810/210‑1 (MQ=25)
                                                                                                                                                       caacaaCCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCACCGCGCCGGTGGTCAATGCCATCGACGCGGTACGCGGCGCAGAt                                                                                                                                                    >  2:65594/1‑163 (MQ=255)
                                                                                                                                                         acaacCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGGTCGAAaacctgaacct                                                                                                                          >  2:470631/1‑187 (MQ=38)
                                                                                                                                                             ccACTCGGAGCTGTACTTCAAGCTGACCGCCGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGcc                                                                                   <  2:1105959/222‑1 (MQ=255)
                                                                                                                                                              cACTCGGCGCTGTACTTCAAGCTCACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                                                                                                              <  1:253045/194‑1 (MQ=37)
                                                                                                                                                               aCTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTa                                                                     <  2:480810/234‑1 (MQ=255)
                                                                                                                                                               aCTCGGAGCCGTACTTCAATCCGACTGACGCGACCCCGATCGTCGTGGCAGCGAACATCACCACCGGTTCGGCCACTGCAAAAGAACCAGACAACGTCTATGTCGGCACCAAcg                                                                                                                                                                                             >  2:1093192/1‑114 (MQ=11)
                                                                                                                                                                 ccGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGa                                                                 <  2:951800/235‑1 (MQ=255)
                                                                                                                                                                  cGGAGCTGTCCTTCAAGCTGACCGCCGGCACCACCGTCGTCGTGGCAGCCACCAGCACCTCCGGTTCCGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcggcg                                                                                                                                                         <  2:282870/147‑1 (MQ=11)
                                                                                                                                                                   ggAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                                                                                                                                                                      >  2:920589/1‑85 (MQ=35)
                                                                                                                                                                   ggAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                                                                                                                                                                                      <  1:920590/85‑1 (MQ=35)
                                                                                                                                                                   ggAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACg                                                                                                                                                                                                               >  2:482846/1‑92 (MQ=37)
                                                                                                                                                                   ggAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACg                                                                                                                                                                                                               <  1:482847/92‑1 (MQ=37)
                                                                                                                                                                   ggAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGCCCGCACCGAGGTGACCGACGCGCCAGGCACACCAGGt                                                                      >  2:380957/1‑229 (MQ=255)
                                                                                                                                                                     aGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAgg                                                              >  1:1111189/1‑235 (MQ=255)
                                                                                                                                                                      gCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGa                                                                                                  <  2:987365/198‑1 (MQ=39)
                                                                                                                                                                           aCTTCAAGCTGACCGACGGCCCCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGt                                                                                                                                                                 <  2:197692/130‑1 (MQ=38)
                                                                                                                                                                                  gCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAAtgcg                                                                                                                                                                                           >  1:803363/1‑97 (MQ=16)
                                                                                                                                                                                  gCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAAtgcg                                                                                                                                                                                           <  2:803362/97‑1 (MQ=16)
                                                                                                                                                                                  gCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCACTGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCGGGACAAGCCCCCGGTCAGCACCGAgg                                                                                                >  1:851842/1‑188 (MQ=255)
                                                                                                                                                                                  gCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                                                                                                              >  1:86397/1‑174 (MQ=38)
                                                                                                                                                                                  gCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAg                                                 >  1:751726/1‑235 (MQ=255)
                                                                                                                                                                                     gACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGGGGTCAACGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGCAACGAAg                                                               >  1:617533/1‑218 (MQ=255)
                                                                                                                                                                                       ccGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATATCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                      >  1:183206/1‑177 (MQ=38)
                                                                                                                                                                                         gACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAg                                                                                                                                           >  2:465237/1‑138 (MQ=37)
                                                                                                                                                                                         gACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                      >  1:904117/1‑175 (MQ=38)
                                                                                                                                                                                         gACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAg                                                                                 >  2:713600/1‑196 (MQ=40)
                                                                                                                                                                                               accaccGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGc                                                             <  1:447648/210‑1 (MQ=255)
                                                                                                                                                                                                           gtGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                                                                                                              <  2:171195/149‑1 (MQ=39)
                                                                                                                                                                                                           gtGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGc                                                                                    >  2:1025301/1‑175 (MQ=255)
                                                                                                                                                                                                            tGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCtt                        <  2:1111188/234‑1 (MQ=255)
                                                                                                                                                                                                            tGGCAGCCAACAGCACCACCGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                                                                                                              <  1:767648/148‑1 (MQ=37)
                                                                                                                                                                                                              gCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAAc                                                                   >  2:834721/1‑189 (MQ=255)
                                                                                                                                                                                                              gCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGt                    >  1:811743/1‑236 (MQ=255)
                                                                                                                                                                                                               cAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCCAGGTGACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACCGCCGACCa                             >  1:648808/1‑226 (MQ=255)
                                                                                                                                                                                                                       aGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGaaaa                                                                                                                                   <  2:22126/116‑1 (MQ=12)
                                                                                                                                                                                                                       aGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCgaga            <  1:1134311/235‑1 (MQ=255)
                                                                                                                                                                                                                               cGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCgc              >  2:575123/1‑224 (MQ=255)
                                                                                                                                                                                                                                 gTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAAcc  >  1:634989/1‑235 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         
TCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACC  >  AE015451/200247‑200706

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