Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,983,147 C→A 9.5% Q64K (CAA→AAA)  aslB → putative AslA‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,983,1470CA 9.5% 215.0 / 2.5 82Q64K (CAA→AAA) aslBputative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (56/17);  new base A (2/7);  total (58/24)
Fisher's exact test for biased strand distribution p-value = 2.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGC  >  NC_000913/3982911‑3983380
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          
aTTTCCCCCTTAACGTTTTTAAGGATGGCTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGGGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCaa                                                                                                                                                                                                                                          <  2:672340/239‑1 (MQ=255)
   tCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCCAGCAAATGGACGATGACACGCTGGAACCGCATGGCCGTCACAATATCGCAGCAAGCGAACGCCaaaaa                                                                                                                                                                                                                                       >  2:101926/1‑238 (MQ=255)
      cccTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTCTCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACCCTGGAAGCGTATGTCCGTCCCTATATCGCGGCCAGCGAACCGCAAAACGaa                                                                                                                                                                                                                                    >  1:561632/1‑239 (MQ=255)
      cccTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCACCGAACGGCAAAACGaa                                                                                                                                                                                                                                    >  1:834915/1‑239 (MQ=255)
      cccTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTACCCATATCGCGGCCAGCGAACCGCAAAACGaa                                                                                                                                                                                                                                    >  1:217510/1‑239 (MQ=255)
      cccTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACCCTGGAAGCGCATGTCCGTCACTATATCGCTGCCGGCGAACCGCAAAACGaa                                                                                                                                                                                                                                    >  1:342336/1‑239 (MQ=255)
           aaCGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTg                                                                                                                                                                                                                                 >  1:403001/1‑237 (MQ=255)
                     ggacaggcTTCGGGCAAGGGGCGACCATTCTTCAACCGGTTCCAAGGCGGGCTTTTAATTTGAGGGGGAACCCGAGCGGTTCCGATTGTAATCTGAACTGTGACTCCTTGTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTgg                                                                                                                                                                                                                     <  2:14820/232‑1 (MQ=255)
                         tGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGGGCGGTTCCGATTGTAATCTGAACTGTGACTACGGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTAACGCATATGGACGATGACACGCTGGAAGCGTAGTACCGTCACTATCTCGCTGCCGGCCAACCGAAAAAAGAAGTGGCTGTTAcacggcgcg                                                                                                                                                                                                                 >  2:426662/1‑231 (MQ=255)
                         tGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGCTGACACGCTTGAACCGAATGTCCGTCACTATATCGCTGCCAGCGACCCGCAAAACGAAGTGGCTTTTAACTTGCAgg                                                                                                                                                                                                                 >  1:517985/1‑239 (MQ=255)
                         tGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTCACGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACCCTGGAAACGTATGTCCGTCCCTATATCGCGGCCAGCGAACCGCCAAACGCAGAGGCTTTAACCTggacgg                                                                                                                                                                                                                 >  2:989114/1‑235 (MQ=255)
                           aCTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGCCCGCAACTATATCGATGCCAGCGAACCGCAAAACGCAGAGGCTTTTACCTGCCAgggg                                                                                                                                                                                                               >  2:942093/1‑238 (MQ=255)
                               cAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCACCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGgcggcg                                                                                                                                                                                                           >  1:965734/1‑239 (MQ=255)
                                      ggAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCaa                                                                                                                                                                                                     >  1:1081255/1‑238 (MQ=255)
                                       gAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGactgcc                                                                                                                                                                                                     >  1:14820/1‑232 (MQ=255)
                                       gAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAAc                                                                                                                                                                                                    >  1:281101/1‑238 (MQ=255)
                                        ggCGGCCATGCTTCAACAGGTTTCATAGCGTTCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGGGACTACTGTTTTTATCTGGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCaa                                                                                                                                                                                                     <  2:217510/235‑1 (MQ=255)
                                        aGCGACCATGCTTCAACAGGGTCCAACGCCGTCTTTTCATGTGGTGGCGAAACCGAGCGGTTCCGATTGTAATCCGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACg                                                                                                                                                                                                   <  2:281101/238‑1 (MQ=255)
                                            aCCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGGTGGCGAAACCGCGCGGTTCCGATTGTAATCTGAACTGTGGCTACTGTTTTTATCTCCAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGTACGCTGACACGACGGAACCGTATGCCCGTCACTATATTCCTGCCACCGAAACGCAAAACGAAGAGGCTTTTACCTGGCAGAGCGGCGCACAAAcacaac                                                                                                                                                                                               >  2:280136/1‑233 (MQ=255)
                                            aCCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCAAACCAACGCAAcc                                                                                                                                                                                              >  2:1014679/1‑238 (MQ=255)
                                                  ctgcAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGCAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGCAAGGGCGGCGAACCAACGcaccccgg                                                                                                                                                                                           >  1:135972/1‑229 (MQ=255)
                                                  ctgcAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCgg                                                                                                                                                                                           >  1:695350/1‑236 (MQ=255)
                                                       aCAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTCTATCCCTGCCAGCGAACCGCAAAACGAATGGGCTTTTACCTGGCAGGGGGGCGAAACAAACCAACACGGGCTggggg                                                                                                                                                                                   >  2:176840/1‑236 (MQ=255)
                                                           gTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGc                                                                                                                                                                                                                               <  1:571715/191‑1 (MQ=255)
                                                             tCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCCCTATATCGCTGCCACCGAACCGCAAAACGACGTGGCTTTTACCTGGCATGGCGGCGAACCAACCCCACCCGGGCTGGCGTTTTAcc                                                                                                                                                                             >  2:93504/1‑239 (MQ=255)
                                                             tCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCCCTGCCAGCGAACCGCCAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAAAGCTACCCGGGCTGGCGTTTTAc                                                                                                                                                                              >  1:26018/1‑238 (MQ=255)
                                                                 aCGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGCCCGTCACTATATCGCTGCCACCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGGGGACCAACGCAACACGGGGTGGCGTTTTAccgccg                                                                                                                                                                         >  2:553616/1‑239 (MQ=255)
                                                                      tGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGACGTGGCTTTTACCTGGCAGGGCGGCGATCCCACGCTACTCGGGCTGGCGTTTTCCCgccgtgccg                                                                                                                                                                    >  1:608601/1‑239 (MQ=255)
                                                                      tGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCgccgtgccg                                                                                                                                                                    >  1:39800/1‑239 (MQ=255)
                                                                              atggaagggCGAACCCGGGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAACCAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTa                                                                                                                                                           <  1:913948/232‑1 (MQ=255)
                                                                                   aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGATAAGCCAGTCACGCATATGGACGATGACACGCTGGAACCGTATGTCCGTCACCATATCGCTGCCAGCGAACCGCAAACCGAAGGGGCTTTTCCCTGGCAGGGCGGCGAACCAACCCTACTCGGGCTTGCGTTTTACCCCAGGGCCGGAgg‑gcccccggc                                                                                                                                                      >  2:473217/1‑232 (MQ=255)
                                                                                   aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTCTGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGACTTTTACCTTGCAGGGCGGCGAACCAACGCTACTCGTGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGc                                                                                                                                                      >  1:293654/1‑239 (MQ=255)
                                                                                   aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGAAACACGAATGGGCTTTTACCTGGCAGGGCGGCGCACCAACGCAACTCGGGCTGGCGTTTTACCGCCGGGCCGGAGC‑GCAACAGGc                                                                                                                                                      >  2:83034/1‑239 (MQ=255)
                                                                                    tgggcaaaCCGGGGGGTTCCAATTGTATCCTGAATTGTACCTTTGTTTTTTATCGCAAAAAACAACCCCTTACCCGCGAAAACCCAGTACCTAATAGGGCCGATGACCCGCTGGAGGCTGATGTCCGTCACTATACCGCTGCCAGCGAACCGAAAAACGAAGTGGCTTTTACCGGGCAGGGCGGCGAACCAACGTTCCCGGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCg                                                                                                                                                     <  2:164626/234‑1 (MQ=255)
                                                                                     ggCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTCGCAGGCCGGCGAACCAACGCTACCCGGGCTGCGGTTTTACCGCCGTGCCGAAGG‑CCCACAGGGGa                                                                                                                                                    >  2:902036/1‑239 (MQ=255)
                                                                                               aGCGGATCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTAACGCATATGGACGATGACACGCTGGAATCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGGACCAACGCCACTCGGGCTGGGGTTTTACCGCCGTGCCGGCGG‑TCAACAGGCGAAAAATGGGGc                                                                                                                                          >  2:907911/1‑239 (MQ=255)
                                                                                                gCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCATCGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTCCAGGCGAAATATggg                                                                                                                                            >  1:158885/1‑235 (MQ=255)
                                                                                                      ctgATTTTAATCTGCAGTGTCACTACTTTTTTTATCACAAAACACAACCCCTTACCCGCGAAAACCCAGTCACGCATGGAGATGATGACCGGCTGGAAGCTTGTGTCCTGCACTATATCGCGGCAAGCAACCCGCAAACGGAAGTGGCTTTTCCTGGGCGGGGCGGCGAACCAACGCTACGCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAgg                                                                                                                                   <  2:1052054/237‑1 (MQ=255)
                                                                                                         gaccaaaTCTGATCGGGTACTTCTTTTTTTATCTCAAAAACCAAGCCCTTCACCACGAAAAGCCAGTCACCAATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTCCCTGGCAGGGCGGCGAACCAACGCTCCTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGa                                                                                                                                  <  2:1051445/232‑1 (MQ=255)
                                                                                                          ttGTAATCTGAACTGGGACTACTGTTTTTATCCCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAg                                                                                                                                <  1:90703/238‑1 (MQ=255)
                                                                                                          ttGAAATCTGAACTGTGACTCCTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACCCGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAg                                                                                                                                <  1:176840/238‑1 (MQ=255)
                                                                                                              aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAACCAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGGATGTCCGTCACTATCTCGCTGCCAGCGAACCTCAAAACGAATTGGCTTTTCCCTGGACGGGCGGCGAACCAACACTACCCGGGATGGCGTGTTACCGCCGTGCCG‑AGCAGCTACAGACCAAACTTGGGGCTGGCAGTAAggaaaa                                                                                                                           >  2:361304/1‑234 (MQ=255)
                                                                                                              aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGG‑CCTACAGGCGACAGATGGGGCTGGCAGGAAGGTAAg                                                                                                                           >  1:756656/1‑239 (MQ=255)
                                                                                                              aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTCCAGGCGAAATATGGTGCTGGCAGGAAGATAAg                                                                                                                           >  1:134723/1‑239 (MQ=255)
                                                                                                              aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAg                                                                                                                           >  1:11000/1‑239 (MQ=255)
                                                                                                                         tgACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTAAAGGCGAAATATGGGGCTGGCAGGAAGATAAGTAACAgggc                                                                                                                  >  1:1051445/1‑234 (MQ=255)
                                                                                                                           actactGTTTTTATCTCGAAAAACAATCCCTATACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGCAATATGGTGCTGGCAGGAAGATAAGTAACAACTTCCAg                                                                                                              >  1:321045/1‑239 (MQ=255)
                                                                                                                             tcctGTTTTAATCTCGAAAAACAATCCTTTTCCCGCGAAAAGCCAGTCACGCATAGGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGACCCGAAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGTTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAAc                                                                                                                       <  1:460422/226‑1 (MQ=255)
                                                                                                                              actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGAAAGATAAGGAACAGCTTCCaacct                                                                                                           >  1:141029/1‑235 (MQ=255)
                                                                                                                              actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAAAATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACt                                                                                                           >  1:54570/1‑239 (MQ=255)
                                                                                                                              actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCCCCGTGCCGTAGC‑GCTAAAGGGCAAATAGGGGGCTGGGGGGGAGAAAAAGAACAGGTTCCAgtaa                                                                                                           >  2:1016045/1‑236 (MQ=255)
                                                                                                                                   ttttATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGGGCCGTAGC‑CCTCCAGGCGAAATTGGGGGCGGGCAGCAAGATAGATAACAGCTTcccaacaaaag                                                                                                       >  2:936633/1‑229 (MQ=255)
                                                                                                                                     ttacccagAAAACCAATCCCTTTACCGCGAAAAGCCAGTACCGCATTGGGACGATGACACGCTGGAAGCGTATTTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCg                                                                                                    <  1:1014679/232‑1 (MQ=255)
                                                                                                                                         ctcGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGAAAAAAGTTCCAAACAAAAGGCGggcg                                                                                                >  2:323812/1‑238 (MQ=255)
                                                                                                                                          tcGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGGGCTGGCAGGGAGAAAAGAAAAAGAGTCCAGACTAAAggggggggg                                                                                               >  2:1098846/1‑232 (MQ=255)
                                                                                                                                              aaaaCACTCCCTTTCCCGAAAAAACCCAGTCCCGCATATGGAGGATGACCCGCGGAAAGGTTTTGTCCGTCATTATATCGCTGCCAGCAAACCGCAAAACGAAGTGGCTTTTACCGGGCAGGGCGGCGAACCAACGCTACGCGGGGTGGCGTTTGACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTgg                                                                                           <  2:566285/239‑1 (MQ=255)
                                                                                                                                                          ttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATg                                                                                 >  1:566285/1‑237 (MQ=255)
                                                                                                                                                           tACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGGGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATTTGGGTCGGGCAGGAAGATAAGTAAAAGCTACAAGCAAAACGGCGAGATGCCGGAGGAGGAAggggg                                                                              >  2:493439/1‑237 (MQ=255)
                                                                                                                                                                cgAAAAGCCAGTCACGCATATGGACGATGACACGCTGGACGCGTATGTACGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCtt                                                                                                                                                                                                                             >  2:725145/1‑92 (MQ=255)
                                                                                                                                                                cgAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCtt                                                                                                                                                                                                                             <  1:725145/92‑1 (MQ=255)
                                                                                                                                                                cgAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGt                                                                         >  1:637424/1‑239 (MQ=255)
                                                                                                                                                                 gAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGTTGACGAATGGGGCGCGtt                                                                        >  1:361870/1‑239 (MQ=255)
                                                                                                                                                                          gTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                               >  1:288179/1‑239 (MQ=255)
                                                                                                                                                                           tCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTCCTGCTGGATGACGAATGGTGCGCGTTTTTCCCggcc                                                              >  1:1069415/1‑237 (MQ=255)
                                                                                                                                                                                    tGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTTACAGCTTCCCGACTAACGGCGGTCTGCTGGATGACGAAAGGGGCGCGTTTTTCCCgggagctcctttt                                                     >  2:469494/1‑228 (MQ=255)
                                                                                                                                                                                           tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATATGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGGGCGCGTTTCCCGCGGAGCATCATTTTCttgttg                                              >  1:424073/1‑239 (MQ=255)
                                                                                                                                                                                           tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgttg                                              >  1:305942/1‑239 (MQ=255)
                                                                                                                                                                                           tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATACTTTTCttgttg                                              >  1:938586/1‑239 (MQ=255)
                                                                                                                                                                                           tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGGAACAGCTTCCCGACTAAAGGCGTGGTGCTGGATGACGAATGGTGCGCGGTTCCCGCGGGGCATCATTTTCtttttg                                              >  2:1061349/1‑239 (MQ=255)
                                                                                                                                                                                                             tgtcctgctttattTCGCGCCCCGCGACCCTAAAACCGGAGTGGCTTTTACTGAGCAGGCGGGCGAACCACGGCTACCGGGGTGGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                            <  2:468053/229‑1 (MQ=255)
                                                                                                                                                                                                             tGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGGGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGGGCATCATTTTCTTGTTGGTTTTACGCTGGATgtc                             >  1:830146/1‑236 (MQ=255)
                                                                                                                                                                                                             tGTCCGTCACCATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGGGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGCAGATACGGAACAGCTTCCAGACTAACGGCGGGCTGCTGACTGGCGAAGGGGGCGCGTTTCTCGCGGGGCATCATTTTTTGGTTGGGTTATcagtggggggg                             >  2:972932/1‑228 (MQ=255)
                                                                                                                                                                                                                   tCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTCACGGCGTGCTGCTGGTTGACGAATGGGGGGCGTTTCTCGCGGGGCATACTTTTTTTGTTTGTTTTTCCATGGATGGCCCGcccg                      >  2:372959/1‑237 (MQ=255)
                                                                                                                                                                                                                         tatCGCCGCCAGCAAACCGAAAAACGAAGTGGCTTTTCCCTGGCAGGGCGGCGAACCAACGCTACTGGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                <  1:323812/239‑1 (MQ=255)
                                                                                                                                                                                                                         tagcggtcccACCGCACCGAAAACCGTGGTGGCTTTTCCCTGGCGGGGCGGCGACCCAACGTTACTCGGGCGGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                <  1:133531/231‑1 (MQ=255)
                                                                                                                                                                                                                         taccGCTGCCGGCAAACGGAAAAACGAATTGGTTTTTCCCGGGCGGGGCGGCGAACCAACGCTACTCGGGCTGGCTTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCGCGCCTGAGATcc                <  1:936633/236‑1 (MQ=255)
                                                                                                                                                                                                                         aaccGCTCCCAGCGAACCGCAAACCGAGGGTGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                <  1:101926/236‑1 (MQ=255)
                                                                                                                                                                                                                             gCTGCAAGCGAACCGCAAAGCGAAGGGGCTTGTCCCTGGCAGGGCGGCGAACCAATGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaa            <  1:907911/239‑1 (MQ=255)
                                                                                                                                                                                                                             gCTCCCGGCAAACCGCAAACCGAAGTGCTTTTTCCCTGGAAGGGCGGCAAACCAACGCTACGCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaa            <  1:942093/239‑1 (MQ=255)
                                                                                                                                                                                                                              cTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTaa                                                                                                                        >  1:188662/1‑130 (MQ=255)
                                                                                                                                                                                                                                        aCCGCAAAACGAAGTGGCTTTTACCTGGAAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc  <  1:181482/238‑1 (MQ=255)
                                                                                                                                                                                                                                        aCCGAAAACCGAGGGGGTTTTTCCTTGGCAGGGGGGCGAACCAACGCTACCCGGGCGGGCTTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCATACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc  <  1:93504/238‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          
ATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGC  >  NC_000913/3982911‑3983380

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑