Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,306,035 C→T G188D (GGC→GAC)  yjcS ← linear primary‑alkylsulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,0350CT98.8% 246.7 / NA 81G188D (GGC→GAC) yjcSlinear primary‑alkylsulfatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/47);  minor base G (0/1);  total (33/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCC  >  NC_000913/4305788‑4306278
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
tCTGACGATAGTTTTCGTCGGTGCAATAATGCTCAGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACACACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGCACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGtcac                                                                                                                                                                                                                                                   >  2:273422/1‑251 (MQ=255)
   gACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCAt                                                                                                                                                                                                                                                >  1:180042/1‑251 (MQ=255)
   gACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCAt                                                                                                                                                                                                                                                >  1:180377/1‑251 (MQ=255)
          gtttcgggcgGTGCAATAAGGCCGGGGTCGCCCGTGGCCAGCGTCCCGCCAAGGCCATTCCCCACATTGCCTTGCGCGTTGTGCGGCAGTACCACACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTc                                                                                                                                                                                                                                         <  2:466287/243‑1 (MQ=255)
          gTTTCCGTCGGTGCAATAATGCCCGGGTCGCCCTTGGCCAGCGTCACGCCAAGGCAATTACCCCCATTGCCTTGCGCGTTGTGCGGCAGAAACAGCCCGTAAGAGTAGAGCGCACGGGGGCTCATGATTTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTc                                                                                                                                                                                                                                         <  1:452643/251‑1 (MQ=255)
            gtccgtcgGTGCAAAAATGCTAGGGTCGCCAGTTGCCAGCGTCACGCCAAGGCCATTCACCACATTGCCTGGCGGGTTGTGCGGCAGTAAAAGACCGTAAGAGTAGAGCGCACGGCGGCGAAGGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTAGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCg                                                                                                                                                                                                                                        <  1:497668/247‑1 (MQ=255)
             ttCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGCGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGgt                                                                                                                                                                                                                                      >  1:105170/1‑251 (MQ=255)
             ttCGTCGGTGAAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCAAAGGCCATTACCCACATTGCCTTGGGGGTTGTGGGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCgg                                                                                                                                                                                                                                       <  2:219349/250‑1 (MQ=255)
                  gggTGCAATAATGCTCGGCTCGCCCGTTGCCAGCGTCACCCTAAGCCCTTTACCCACATTGCCTGGCGCGTTGTGCGGCAGTAACAGACCGGAAGAGTAGAGCGCACGGCGGCTCAGGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTTCATAAAGCCTGCAGGGGCAATCACCTGAGCTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGc                                                                                                                                                                                                                                 <  1:558777/250‑1 (MQ=255)
                    gTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCtg                                                                                                                                                                                                                               >  1:477617/1‑251 (MQ=255)
                    gTGCAATAATGCTCGGGTCGCCCGTTCCCAGCGCCCCGCCAAGGCCTTTACCCACATTGCCTTGCGCGTTGTGCGGCAGAAACAGACCGTAAGAGTAGAGCGCAGGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCtg                                                                                                                                                                                                                               <  1:335840/251‑1 (MQ=255)
                    gTGCAATAATGCCCGGGTGGCCGGTGGCCAGCTCCACGCCAAGGCCTTTACCCACATTGCCTTGGGCGTTGTGCGGGAGTAACAGACCGATAGAGTAGAGCGCACGGCTGGTCAGGATTTTGCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCTATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCtg                                                                                                                                                                                                                               <  1:602882/251‑1 (MQ=255)
                            attcTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCGTGCGCGTTGGGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTa                                                                                                                                                                                                                       <  2:215451/248‑1 (MQ=255)
                              gCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAgtgt                                                                                                                                                                                                                        >  1:554858/1‑248 (MQ=255)
                              gCTCGGGACGCCCGGTGACCGCGTCACGCCAAGGCCATTACCCACATTGACTTGAGCGTTGTGCGGCAGTAACAGACCGTAAGCGAAGCGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGATGATGGCATCGTCCATAAAGCCTGCAGGGGCAATCACCTGCACTTTGCCGAATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTAGGTGTGGCTGTGAgtgt                                                                                                                                                                                                                        >  2:557436/1‑248 (MQ=255)
                                 cGGGTCGCCCGTTGCCAGAGTCAAGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCCGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTCCAGGGGCAATCCCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTAGCTGTGGGTGTAGATa                                                                                                                                                                                                                   >  2:574333/1‑250 (MQ=255)
                                  gggTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTCCCCACATTGCTTTGCGCTTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGAGGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGAGTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAAc                                                                                                                                                                                                                 <  1:269162/251‑1 (MQ=255)
                                    gTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAAAGCCTTTCACGTCCCCATCGTGGTCGGTGTGGCTGTGAGTGTAGATACCgg                                                                                                                                                                                                               >  2:509475/1‑251 (MQ=255)
                                     tCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGc                                                                                                                                                                                                              >  1:212112/1‑251 (MQ=255)
                                     tCGCCCGTTGACAGCGTCAAGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGTCAGTAAAAGACCGTCAGAGTAGAGCGCACGGCGGCACATGATGTTACCCGACAACACGTTTTCGCGCATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCCGAACATTGCCGGATTTAACAGCGGCGTCAGGGATAATGCCTTTCAAGTCACCATATTAGTCGGAGTGACTGTGAGTGTAGACAActgc                                                                                                                                                                                                              >  1:297331/1‑248 (MQ=255)
                                                        cGCCAAGGCCATTACCCCCATGGCCTTGCGGGTTGTGCGGCAGTAACAGACGGTAAGGGATGGGCCCTGGGCGGCCTATGAGTTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGa                                                                                                                                                                                           <  2:395299/251‑1 (MQ=255)
                                                        cGCCAAGCCCAGTCCCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGAAGAGCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGa                                                                                                                                                                                           <  1:315040/251‑1 (MQ=255)
                                                            aaGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGAt                                                                                                                                                                                       >  2:316368/1‑251 (MQ=255)
                                                                   ttACCCACGTGGCCTGGCGGTTTGTGCGGCAGTAACGAACCGTAAGAGTGGAGCGCCCGGCGGCTAATGATTTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCt                                                                                                                                                                                    <  2:491549/247‑1 (MQ=255)
                                                                   ttACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCt                                                                                                                                                                                    >  1:491549/1‑247 (MQ=255)
                                                                    tccccACATTGCCTTGCGCGTTGTGCGGCAGTAACAGCCCGTAAGATTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGaa                                                                                                                                                                                <  1:211848/248‑1 (MQ=255)
                                                                       cccctgtgcccggcgcGGTGTGCGGCAGTAACAGACCGTAAGGGTAGAGCGCACGGCGGCTCATGATTTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTaa                                                                                                                                                                            <  2:344142/239‑1 (MQ=255)
                                                                                 ttgagcgTTGTGCGGCAGTAACAGACCGTCAGAGTAGAGAGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGCCGATGGTGGACGTAAAGGTCAAGTg                                                                                                                                                                  >  2:305307/5‑251 (MQ=255)
                                                                                  ggcgcGTGGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTATGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGc                                                                                                                                                                 <  1:519002/250‑1 (MQ=255)
                                                                                    cgggTTGTGGGGCAGTAAAAGCCCTAAAGAGTAGAGCGCCCGGCGGCGCATGATTTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCgg                                                                                                                                                               <  1:577714/248‑1 (MQ=255)
                                                                                     gcgTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATTATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTGTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGc                                                                                                                                                              <  1:509475/251‑1 (MQ=255)
                                                                                     gcgTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAAAGGCAACACTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCACGTGCGGc                                                                                                                                                              >  1:243565/1‑251 (MQ=255)
                                                                                      cgTTGTGCGGCAGTAACAGCCCGTAAGAGTAGCGCGCACGTCGGCTCATTATTTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCt                                                                                                                                                             <  1:441665/251‑1 (MQ=255)
                                                                                       gTTGTGCGGCAGTAACAGACCGAAAGAGTAGAGCGCACGGGGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCtt                                                                                                                                                            <  1:539185/251‑1 (MQ=255)
                                                                                          gggCGGCAGTAACAGACCGTAAGAGGAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGTACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTg                                                                                                                                                         <  1:177863/249‑1 (MQ=255)
                                                                                          gggCGGCAGAAACAGACCGTAAGAGGAGAGCGCACGGCGCCTCATGAGGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTg                                                                                                                                                         <  2:11829/249‑1 (MQ=255)
                                                                                                  gTAACAGACCTGAGGAGTAGAGCGCACGGCGGTTCGTGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGGGCGTCTTTTGCGGCAggc                                                                                                                                                 <  2:391677/251‑1 (MQ=255)
                                                                                                    accacaccacaagggaagagCGCACGGCGCCGCAAGATGTGACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTCTGCGGCAggcgg                                                                                                                                               <  2:245890/235‑1 (MQ=255)
                                                                                                        gACCGAAAGATTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTc                                                                                                                                           <  1:210403/251‑1 (MQ=255)
                                                                                                                  gTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGt                                                                                                                                 <  1:316368/251‑1 (MQ=255)
                                                                                                                    agagCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTGAcccagaccgg                                                                                                                                >  2:429124/1‑242 (MQ=255)
                                                                                                                    agagCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCa                                                                                                                                       <  1:429124/243‑1 (MQ=255)
                                                                                                                            aggcggtTCATGTTGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTTAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaat                                                                                                                       <  1:223984/250‑1 (MQ=255)
                                                                                                                             ggcggcTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCGACAATCGGTTTTTTCGGCCAATGCTGGAAGTAAGGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAAAaat                                                                                                                       >  2:593614/1‑250 (MQ=255)
                                                                                                                                 ggatatgatgtacgccggcAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAAGGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGGGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTTCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATgt                                                                                                                  <  2:25195/233‑2 (MQ=255)
                                                                                                                                    cattattttccccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGTCGTCACCAGCGGGTCGATAACAATAATGCCt                                                                                                                <  1:520158/240‑1 (MQ=255)
                                                                                                                                        attatccccGCCAGCACGTTTTCGCTGATGGCTACGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTtctc                                                                                                           <  1:344367/245‑1 (MQ=255)
                                                                                                                                        atgTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTAACCAGCGGGTCGATAACAATAATGCCTTtctc                                                                                                           >  1:30347/1‑251 (MQ=255)
                                                                                                                                            tccccGCCGGCACGTTTTCGCTGATGGCTTCGTCCATAAGGCCTGCAGGGGCAATCACCTTAACCTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGGCACCATAGTGGTCGGTGTGGCTGGGAGTGTAGATAACGGCAACAAGCGGTTTTTGCGGACGGTGCTGGAAGTAAAGGTCAAGTGCGGTTTTTGCGGCAGGCGGCGTCACCAGAGGGTCGATAACAATAATGCCTTTCTCGcc                                                                                                        <  2:489505/248‑1 (MQ=255)
                                                                                                                                            tccccGCCAGCACGTTTTCGCTGAGGGCTGTGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGcc                                                                                                        <  1:319205/248‑1 (MQ=255)
                                                                                                                                            tACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGcc                                                                                                        <  1:239635/250‑1 (MQ=255)
                                                                                                                                                         gTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTa                                                                                           >  1:64135/1‑250 (MQ=255)
                                                                                                                                                           tttCGCTGATGGTTTCGTCCATAAAGCTTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt                                                                                         <  1:372755/250‑1 (MQ=255)
                                                                                                                                                              cgcagATGGCTACGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTa                                                                                     <  2:282639/247‑1 (MQ=255)
                                                                                                                                                               gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCt                                                                                                                                                                                                                                >  1:341991/1‑111 (MQ=255)
                                                                                                                                                               gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCt                                                                                                                                                                                                                                <  2:341991/111‑1 (MQ=255)
                                                                                                                                                               gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTGTGCCGGATTTAACATCGGCTTCAGAGATAATTCCTTTCACGTCACCATAGTGGTCGGTGGGGCTGTGAGTGTAAATAACGGCAACAATCGGTTTTTGCGGACGAGGCTGTAAGTAAAGGTCAAGTTCTTCTTTTGCGTCAGGCGGCGTCACCATCGGGTCGATAACAATAATGCCTTTCTCGCCCCCAACGATCGTAATGTTag                                                                                    <  2:387651/251‑1 (MQ=255)
                                                                                                                                                                     ggCTTCGTCCATAAAGCCTGCAGTGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGGACGTAATGTTAGAGatat                                                                              >  1:24096/1‑251 (MQ=255)
                                                                                                                                                                      gCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTGTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGTCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGGAATGTTAGAGATATc                                                                             <  2:100445/251‑1 (MQ=255)
                                                                                                                                                                         ccggccATAAAGCCTGCAGGAGCAATCACCTTAACTTTGCCGGATTTAACATCGGCTTCATAGATAATGCCTTTCTCGTCCCCGTAGTGGTCGGTGGGGCTTTGAGTGTAGATAACGGCAGCAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTTGGGCAGGCGGCGTCAGCAGCGGGTCGATACCAATAATTCCTTTCTCGCCCTGAAGGTACGTATTGTCAAAGATATCCt                                                                           <  2:420194/246‑1 (MQ=255)
                                                                                                                                                                               aTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCcgcg                                                                    <  2:212112/251‑1 (MQ=255)
                                                                                                                                                                                 aaaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCCGCTTCCGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATTTTCGAGATATCCTGGCCgcgcc                                                                  >  1:314843/1‑250 (MQ=255)
                                                                                                                                                                                        gCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGa                                                             <  2:477617/249‑1 (MQ=255)
                                                                                                                                                                                          aGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATAc                                                          >  1:68430/1‑250 (MQ=255)
                                                                                                                                                                                              gCAATCACCTGAACTGTGCCGGATTTAACATCGGCTTCAGAGATAATGCGTTGCACGTCACCATAGTTGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGTCAGGCGTGGTCACCAGCTGGTCGATAACAATAATGCCTTTCTCGCTCTCAACGAGCGGAATGTTAGAGATATCGGTGCCGCTTACCGGATACAtttt                                                    <  2:261741/251‑1 (MQ=255)
                                                                                                                                                                                                aaTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAAAAATAATGCCTTTCTCGCCCTCCACGAACGGAATTTTAGAGATCTCCCGGCCCCGCACCTGATACCTTTTAt                                                  >  2:287458/1‑251 (MQ=255)
                                                                                                                                                                                                   aaCCTGATCTTTGTCGGAGTTAACATCGTCTTCAGAGATAATGCCTTTCACGTCCGTATATTGGTCGGTGTTGCTGTGAGTGTAGATAACGGCAACAATGGGTTTTTGCGGACGATGCTGGAGGTAAAGGTCATGTGCGGCTTTTGCGTCAGGCGGGGTCACCAGCGGGTCGATAACAATAAAGCCTTCCTCCCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTCATACATTTTATCg                                                <  2:332811/249‑1 (MQ=255)
                                                                                                                                                                                                    cccTGAACTTTGCCGGATTTAACACCGGCTTCAGAGAGAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCGGGAAGTAAAGTGCAAGTGCGGCTTTGGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCGTTCGCGCCCTCAACGAACGTAAGGTGATAGATATCCTGTCCGCGCACCTGATAGATTTTGTCGGt                                              <  2:562549/250‑1 (MQ=255)
                                                                                                                                                                                                    aCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGt                                              <  1:305591/251‑1 (MQ=255)
                                                                                                                                                                                                      cTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTg                                             >  1:175199/1‑250 (MQ=255)
                                                                                                                                                                                                      cTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCCGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTg                                             >  1:174851/1‑250 (MQ=255)
                                                                                                                                                                                                       tGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCCGCGGGTGGGAAACAATAATGCCCTTTTCGCCCTCAACGGACGGAAGGCTAGAGATATCCCGTCCCCGCACCTTCAACATTTTTTCGGTGAc                                           >  2:281964/1‑251 (MQ=255)
                                                                                                                                                                                                          aCTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTGCGCTTCACCATAGTGGTCGGTGTGGCTGTGAGTGGAGATAACGTCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCGTTTGCGGCAGTCGGCGTCACCAGCGGGGCGAGAACGATGATGGCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTTGCCGCGCACCTGATACATTTTATCGGTGACttt                                        <  1:597873/251‑1 (MQ=255)
                                                                                                                                                                                                                     tttAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGGTACATTTTCTTGGTGCCTTTGAACAGGcc                               >  1:32932/1‑249 (MQ=255)
                                                                                                                                                                                                                        aaCATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATACCGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGACGATAACAATACTGCCTTTCTCGCCCTCAACGAACGGACTGTTAGAGATATCCTGGCCGCGCCCCCGATACATTTTATCGGTGACTTTGAACATGCCAGaa                           >  1:373315/1‑250 (MQ=255)
                                                                                                                                                                                                                           aTCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGGGGCGTTACCAGAGGGTCGATAAAAATAATGCCTTTTTCGCCCCCCACGACCGTAATGTTTGAAATATACTGGCCGCGCACCTGATACATTTTTTCTGTGGCTTTggaccggaccggaacac                       >  2:566522/1‑235 (MQ=255)
                                                                                                                                                                                                                           aTCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATAc                       >  1:136809/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCCGAAATACCGt                    >  1:576708/1‑251 (MQ=255)
                                                                                                                                                                                                                                      agaTAATGCCTTTCACGTCTCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCAGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGc            <  2:524317/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCACCAATCGGTTTTTGCGGACGATGCTGGAATTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGCCGGTCGATAACAATAATACCTTTCTCGCCCCCCACGAACGTAATGTTGAAGATATCCTTGCCGCGCACTTGTTACCATTTTTCGGTGAACTTGAACAAGGCAGAAATAcccgttattttggcc         >  2:336417/1‑238 (MQ=255)
                                                                                                                                                                                                                                          aaTGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCGCGCCCTCAACGAACGTAATGTTATAGATATCCTGGCCGCGCACCTGATACATGTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGAc         <  1:273422/250‑1 (MQ=255)
                                                                                                                                                                                                                                                tttCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGcc  >  1:518077/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
TCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCC  >  NC_000913/4305788‑4306278

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