Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 90,432 A→G E113E (GAA→GAG)  rsmH → 16S rRNA m(4)C1402 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091390,4320AG98.8% 258.9 / ‑4.4 82E113E (GAA→GAG) rsmH16S rRNA m(4)C1402 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (39/41);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

GTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAG  >  NC_000913/90191‑90676
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
gTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCttt                                                                                                                                                                                                                                              <  1:473746/251‑1 (MQ=255)
    cgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCt                                                                                                                                                                                                                                          >  1:124069/1‑251 (MQ=255)
    cgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCt                                                                                                                                                                                                                                          <  1:299005/251‑1 (MQ=255)
      gtggtCACTCACGTCTGATCCTCTCGCAGCTTGGCGGAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCCATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGAGCTGGCCGAATACGTTGCCGAGCGCGATATTATCGGCAAGGTCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCttt                                                                                                                                                                                                                                        >  1:624133/1‑251 (MQ=255)
      gtggtCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCttt                                                                                                                                                                                                                                        >  1:45042/1‑251 (MQ=255)
           cactcacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGCTCGCTCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGTGCGCGATCTTATCGGCCAGCTCGACGGCATTC‑CACTCGATCTTGGCGTCTCTTCCACGCAAATTGATGATGCTGAGCGCGGCTTTTCCTTTATgcg                                                                                                                                                                                                                                   >  2:63337/1‑251 (MQ=255)
               cacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGAt                                                                                                                                                                                                                               >  1:648552/1‑251 (MQ=255)
                ccGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATg                                                                                                                                                                                                                              <  2:17358/250‑1 (MQ=255)
                acGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGGGCGCGGCTTTTCCTTTATGCGCGATg                                                                                                                                                                                                                              >  2:292778/1‑251 (MQ=255)
                  gTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGt                                                                                                                                                                                                                            >  1:386198/1‑251 (MQ=255)
                  gTCTGAGCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGt                                                                                                                                                                                                                            <  2:187351/251‑1 (MQ=255)
                             tcGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACAt                                                                                                                                                                                                                 >  1:483831/1‑251 (MQ=255)
                               gCAGCTTGGCGAAGGGGGGCGTTTGCGGGCGAGCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCAGCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGc                                                                                                                                                                                                               <  2:388518/251‑1 (MQ=255)
                                 aGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGGAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCtttt                                                                                                                                                                                                                                             >  2:25433/1‑219 (MQ=255)
                                 aGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCtttt                                                                                                                                                                                                                                             <  1:25433/219‑1 (MQ=255)
                                  gctgggCGAAGAGGGGCGTTTGCTGGGGATCGATCGCGACCCGCAGGCTATCGGCATTGCGAGGACTATT‑GAGGATCCGCGCTTCTCCATCAGCCACGGACCTTTTGCCGCGCTGGTCGAATACGTTGCCGAGCGCGATCTGATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTa                                                                                                                                                                                                            <  2:578879/247‑1 (MQ=255)
                                        cGAAGAGGTGCGTTTGCTAGCGATCGATCGCGACGCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCGTCTACAGCATCCACGGACCCTTCGCCGCGCTGGCAGAATACGTTGCCGAGCGCGATCTTATCGGCAAGAGCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCGTTTCCTTGAGGCGCGATGGTCCGCTGGACATGCGTATGGAcc                                                                                                                                                                                                      <  2:276977/251‑1 (MQ=255)
                                          aaGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGAccc                                                                                                                                                                                                     >  2:193198/1‑250 (MQ=255)
                                            gagGGGCGTTTGCGGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAAc                                                                                                                                                                                                  <  2:488431/251‑1 (MQ=255)
                                             atgggCGTGTGCTGGCGAGCGATCGCGACCCGCAGGCTATCGACGTTGCGAGGACTATT‑GATGACCCTCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTCGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTCTCCGTGATGCGCGATGGTCCGCTGGACATTCGTATGGTCCCAAcc                                                                                                                                                                                                 <  2:104558/249‑1 (MQ=255)
                                              gggggttttGATGGCTATCGATCGCGCCCCGCAGGCTATCGCCGTTGCGAAGACTAGT‑GAGGATCCGCGCTTCTCCATCATCCACGGACCTTGCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAAccc                                                                                                                                                                                                <  2:686845/245‑1 (MQ=255)
                                               gggCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCg                                                                                                                                                                                               <  1:372519/251‑1 (MQ=255)
                                                   gTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTggg                                                                                                                                                                                           >  1:629906/1‑251 (MQ=255)
                                                    gtttcTGGCGATCGAGAGCGACCCGTAGGCTTTCGTAGTTGTGAAGACTATT‑GATGATCCGCGCTTCTCCATAATCCACGGCCCGTTCTCCGCGCGGGGCGAAGACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCCCCGCATCTTGATGATGCTGAGCGCGGCTTTTCCTGTATGCACTATGGCCCGTTGGACATGCGTATGGACCCAACCCGTgg                                                                                                                                                                                            <  2:533071/245‑1 (MQ=255)
                                                     ttGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGc                                                                                                                                                                                          <  1:519821/250‑1 (MQ=255)
                                                      tGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGc                                                                                                                                                                                          >  1:288159/1‑249 (MQ=255)
                                                        gtGGCGATCGGTCGCGACCCGTAGGCTACCGCCGTTGCGAAGACTTAT‑GATGACCCGCGCTTCTCCATCAGTCACGGACCTTTCGCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGAGCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCTCCGCAACTTGATTATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGGCCGCTTGACATGCGTATGGACCCACCCCGTGGTCAGTc                                                                                                                                                                                      <  2:124069/250‑1 (MQ=255)
                                                            cgatcgatCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGcc                                                                                                                                                                                  <  1:49648/251‑1 (MQ=255)
                                                             gatcgatcGCGACCCGCAGGCTATCGCCGTTGCTAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCATCCCGTGGGCAGTCAGCCg                                                                                                                                                                                 <  2:648552/251‑1 (MQ=255)
                                                             gatcgatcGCGAACAGCAGGCTATCGACGTTGAGAAGACTATT‑GAAGATACGCGCTTATCCATCATCCACGGACCTATCTCCGTGCTGGGAGAATAAGGTGCCGAGCGCGCTCGTATCGGAAAGATCGCCGGCAATC‑TCCTCGATCTTGGCGTCTCTTCACCGCACCTTGATTCTGCTGAGCGCGGCATTTCCTTTATGCGCGATGGTACGCTGGAAAGACGTATTGACCCAACCCGTGGGCAGtgaaccg                                                                                                                                                                                 >  2:120666/1‑245 (MQ=255)
                                                                cgatcgCGGCCCGCGGGCTATCGCCGTGGCGAATACTATT‑GATGGTCCGCGCTTCTCCATCATCCACGGACCTTGCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑GCCTCCATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTg                                                                                                                                                                              <  2:434977/251‑1 (MQ=255)
                                                                      cgACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGCTCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCACCCCGTGGGCAGTCAGCCGCTGAATGGc                                                                                                                                                                        >  2:47219/1‑251 (MQ=255)
                                                                      cgACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCAGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCACGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGGTGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACCTGCGTATGGACCCAACCCGTGGGCAGCCCTCCGCTGAGTGGc                                                                                                                                                                        >  2:83344/1‑251 (MQ=255)
                                                                       gACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCt                                                                                                                                                                       >  1:57625/1‑251 (MQ=255)
                                                                       gACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GAGGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCt                                                                                                                                                                       >  1:183089/1‑251 (MQ=255)
                                                                        aCCCCCAGGCTTTCGGCGGTGCGAAGACCATTTGA‑GATCCGCGCCTCCTCATCAGACCCGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGAGCTTATCGGCAAGATCGGAGGCCTTCCTCC‑CGATCTTGGCGTCTCTTCACCGTAACCTGATGATGCCGAGCGGCGTCTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTa                                                                                                                                                                      <  2:80628/251‑1 (MQ=255)
                                                                           cGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACa                                                                                                                                                                    >  2:133734/1‑250 (MQ=255)
                                                                            gCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCTCTTCTCCATCATCCACGGACCTTGCTCCGCGCTGGGCGAATACGGTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTGGGCGTCTCTTCACCGCAGCTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGTTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACaaa                                                                                                                                                                  <  2:57051/251‑1 (MQ=255)
                                                                            gCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACaaa                                                                                                                                                                  <  1:398117/251‑1 (MQ=255)
                                                                            ccAGGCTATCGCCGGTGCGAGGTATATT‑GATGGTACGCGCTTATCCATCATCCACGGACCGTTCTCCGCGCTGGGCGAATACGGTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAATTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCGGGACATGTGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACaaa                                                                                                                                                                  <  2:139852/250‑1 (MQ=255)
                                                                                gCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGc                                                                                                                                                              <  1:193198/251‑1 (MQ=255)
                                                                                             cGAAGACTATT‑GATGATCCGCGCTTCTCCATCAGCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGat                                                                                                                                                 <  1:47219/251‑1 (MQ=255)
                                                                                              gAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGata                                                                                                                                                <  1:393685/251‑1 (MQ=255)
                                                                                              gAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGat                                                                                                                                                 >  1:138467/1‑250 (MQ=255)
                                                                                               aaGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGCCCTTTCTCCGCGCTGGGCGCATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCAGCTGGACATGCGTATGTACCCAACCCGTGGGAAGTCAGCCGCGGAATGGCTACAAACCGCAGAAGAAGCCGctat                                                                                                                                               >  2:638997/1‑251 (MQ=255)
                                                                                               aaGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGatat                                                                                                                                               <  1:291725/251‑1 (MQ=255)
                                                                                                       t‑gatgaTCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGGAGCCGAGATCGCCTggg                                                                                                                                       >  1:133514/1‑251 (MQ=255)
                                                                                                          atgatcggtgcctctcCATCAGCCTCGTACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGTGCGATCTGATCGGCAAGATCGACGGCAGTC‑TCCTCGATCTTGGCGTCTCTTGACCGAAACTGGATGATGCTGATCGCGGCTTTTCTTTTAGGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGGGGTCAGGCATCCGCTGAATGGCTACAAATCGCCTAAGAAGCCGATACCGTCTGGGt                                                                                                                                      <  2:457632/238‑1 (MQ=255)
                                                                                                          aggttCCGTGCTCCTCCCTCATCCACTGTCCTTTCCCCGCGCTGGGTGAATACTTTTCCGAGCGCGATCTTACCGTCAAGATCGACGGCATGC‑TCCCCGAGCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCTCGGGACATTCGTATGGGCCCAACCCGTGTGCAGTCAGCCTCTGAATGGCTACAAACCGCAGAAGTAGCCGATATCGCCGGGGTa                                                                                                                                     <  2:638987/247‑1 (MQ=255)
                                                                                                             atCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTg                                                                                                                                  >  1:123751/1‑251 (MQ=255)
                                                                                                                cgcgCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCACCCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGaa                                                                                                                                >  2:183155/1‑250 (MQ=255)
                                                                                                                  cgcTTCTCCACCATCCTCGGACCGTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTAGCGGCAAGATCGACGGCATTC‑CCCTCGATCTTGGCGTCTCTTCACCGCAACGTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAAc                                                                                                                             <  2:553456/251‑1 (MQ=255)
                                                                                                                   ccTTCTCCATCATACCCGTACCTGTCTCCGTGCGGGGCGAATTCGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGAGGGTCCGCTGGACATGCGTATGGCCCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAAcc                                                                                                                            <  2:386229/250‑1 (MQ=255)
                                                                                                                     ttccccATCATCCCCGGACCTTTCTCCGCGCGGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGaagaag                                                                                                                                                      <  2:53044/219‑1 (MQ=255)
                                                                                                                     ttCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGaagaag                                                                                                                                                      >  1:53044/1‑223 (MQ=255)
                                                                                                                         ccATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGAGCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGTTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATgg                                                                                                                       <  1:83344/250‑1 (MQ=255)
                                                                                                                              atcCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGa                                                                                                                    >  1:453738/1‑248 (MQ=255)
                                                                                                                                 cACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACAGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCg                                                                                                              >  1:278464/1‑251 (MQ=255)
                                                                                                                                   cGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGCTATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCggt                                                                                                            >  2:19779/1‑249 (MQ=255)
                                                                                                                                    ggACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGGGAAGAGCGttt                                                                                                           >  2:207276/1‑251 (MQ=255)
                                                                                                                                      aCCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTg                                                                                                         >  1:632228/1‑251 (MQ=255)
                                                                                                                                                           aTACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGcc                                                                                    >  1:30079/1‑250 (MQ=255)
                                                                                                                                                                    cGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGc                                                                          <  1:68661/251‑1 (MQ=255)
                                                                                                                                                                                    ggCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGGGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACacac                                                          >  2:84375/1‑248 (MQ=255)
                                                                                                                                                                                     gCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCg                                                         >  2:173367/1‑251 (MQ=255)
                                                                                                                                                                                         gATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTCTGGCCCCAACCCGTGGGCAGTCAGCCGCTGCATGGCTACAAACCGACGTAGAAGCCGATATCGCCGGGGTATTGAAATCCTAGGGTGAAGCGCGTTTTGCCCAACCGC‑TTGTCCGCGCCATTTTCGCGCTCAACCGCGAACATCCggtg                                                      >  2:565082/1‑247 (MQ=255)
                                                                                                                                                                                         gATCGACGGAATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGCGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAACCCGAAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGAC‑AACCGCATTGCCAGCGCAATTGTCGAGCGTAACCGCGAACAGCCGATg                                                      >  1:68159/1‑250 (MQ=255)
                                                                                                                                                                                          aTCGACGGCTTTC‑TCCCCGAGTTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGAc                                                    <  2:288159/251‑1 (MQ=255)
                                                                                                                                                                                          aTCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTg                                                                                                                     >  1:484223/1‑187 (MQ=255)
                                                                                                                                                                                          aTCGACGGCATTC‑TCCACGATCTTGGCGTCTCTTCCCCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATTGCCTGGGTATTGAAAACCTATGGTg                                                                                                                     <  2:484223/187‑1 (MQ=255)
                                                                                                                                                                                               cGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCa                                               <  1:500017/251‑1 (MQ=255)
                                                                                                                                                                                                   aTTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATAATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCa                                            >  1:200850/1‑250 (MQ=255)
                                                                                                                                                                                                        ccctcGATCTTGGTTTCTCTTCACCGCAATTTGAGGATTCTGAGCGCGGCTTTCCCTTTATGCGCGATGGCCCGCTGGACATGCGTATGGACCCAACCCGTGGGCCGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGaa                                       <  2:629906/250‑1 (MQ=255)
                                                                                                                                                                                                               tCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCgg                                >  1:498775/1‑251 (MQ=255)
                                                                                                                                                                                                                  tggggTCTCTGCACCGCAACTGGGTGATGCTGAGCGGGGTTTTGCCTTTATGCGCGTTGGTCCGCTGGACATGCGAATGGACCCAACCCGGGGGCAGTCAGACGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGGGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGa                               <  2:512010/245‑1 (MQ=255)
                                                                                                                                                                                                                  tGGCGTCTCTGCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTGATGCGCGATGGTCCGCTGGACAGGCGGATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGc                                  <  2:66074/246‑1 (MQ=255)
                                                                                                                                                                                                                        ctctTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTTGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTgg                       >  1:673671/1‑251 (MQ=255)
                                                                                                                                                                                                                         tctTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGgc                      >  1:192337/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCAGAAGTCGTGGCTGCTGCa               <  1:34314/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tgatgCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGTGCGTAACCGCGAACAGCCGATGACCCGCACCACAGAACTGGCGGAAGGCGAGGCTGCTGCAACGACGGTgg     >  2:575939/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gatgCTGAGCGCGGCTTTTCCTTTATGCGCGATGGGCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGa     <  1:571671/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAg  <  1:19779/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATT‑GATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTC‑TCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCC‑AAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAG  >  NC_000913/90191‑90676

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