Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,997,786 T→C K173E (AAG→GAG)  PA14_RS22840 ← DUF853 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,997,7860TC100.0% 1890.6 / NA 475K173E (AAG→GAG) PA14_RS22840DUF853 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (230/245);  total (230/245)

Only 100 of 475 total aligned reads displayed.

GCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGG  >  NC_008463/4997639‑4997932
                                                                                                                                                   |                                                                                                                                                  
 cccccGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGACGCGCCGGTGAACAGCGCCCGGTCCTCGcc                                                                                                                                                 >  1:6644996/1‑150 (MQ=255)
    ccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAg                                                                                                                                              <  2:4133731/150‑1 (MQ=255)
     cgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAgg                                                                                                                                             >  2:3049671/1‑150 (MQ=255)
       cccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAt                                                                                                                                           >  1:1982696/1‑150 (MQ=255)
          gTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAtctc                                                                                                                                        >  2:5476094/1‑150 (MQ=255)
          gTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAtctc                                                                                                                                        <  1:393702/150‑1 (MQ=255)
          gTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAtctc                                                                                                                                        <  2:6872760/150‑1 (MQ=255)
           tCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTcc                                                                                                                                       <  2:3662922/150‑1 (MQ=255)
           tCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTcc                                                                                                                                       <  1:1379887/150‑1 (MQ=255)
             cGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCgg                                                                                                                                     <  1:7149748/150‑1 (MQ=255)
              ggCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGa                                                                                                                                    >  1:7769018/1‑150 (MQ=255)
              ggCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGa                                                                                                                                    >  1:9187447/1‑150 (MQ=255)
                 tCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTg                                                                                                                                 >  2:5238556/1‑150 (MQ=255)
                  cGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGCAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGt                                                                                                                                <  2:1622916/150‑1 (MQ=255)
                    gggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCTTGCGGCGCATCAGCGACTGCCCGGAATCGCCGGTGTACAGCGCCCGGTCCTCGCCGTGGATCTCCGGATTGTcc                                                                                                                              <  2:9707544/150‑1 (MQ=255)
                        ggaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCa                                                                                                                          <  1:8136307/150‑1 (MQ=255)
                         gaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGCACAGCGCCCGGTCCTCGCCGAGGATCTCCGGCTTGTCCTTCAg                                                                                                                         >  2:9330250/1‑150 (MQ=255)
                          aggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAgg                                                                                                                        >  2:39445/1‑150 (MQ=255)
                           ggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggt                                                                                                                       <  2:10453049/150‑1 (MQ=255)
                              tGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGTCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggtggt                                                                                                                    <  1:9865155/150‑1 (MQ=255)
                                      aGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTGGagcagc                                                                                                            >  2:886664/1‑150 (MQ=255)
                                       gTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCTgcg                                                                                                           >  1:3068138/1‑150 (MQ=255)
                                        ttGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAgcgc                                                                                                          <  2:2064822/150‑1 (MQ=255)
                                          gcagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCTGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcct                                                                                                        <  2:4284073/150‑1 (MQ=255)
                                          gcagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcct                                                                                                        <  2:2669513/150‑1 (MQ=255)
                                           cagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcctt                                                                                                       <  2:6129652/150‑1 (MQ=255)
                                            agggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCCGGTGGTGGAGCAGCGccttc                                                                                                      >  1:5759683/1‑150 (MQ=255)
                                             gggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttca                                                                                                     <  1:269301/150‑1 (MQ=255)
                                               gcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACCGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcagg                                                                                                   <  1:6714032/150‑1 (MQ=255)
                                                 aggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCGTCAGGTGGTTGAGCAGCGccttcaggtc                                                                                                 >  1:3592064/1‑150 (MQ=255)
                                                      cGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttca                                                                                            >  1:3851020/1‑150 (MQ=255)
                                                      cGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttca                                                                                            <  1:1208104/150‑1 (MQ=255)
                                                       gCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcag                                                                                           >  1:8622279/1‑150 (MQ=255)
                                                           aagaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCg                                                                                       >  1:10028127/1‑150 (MQ=255)
                                                              aagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagc                                                                                    <  1:7971349/150‑1 (MQ=255)
                                                                ggCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcag                                                                                  >  1:4570650/1‑150 (MQ=255)
                                                                   ctaggCGCCCTGCTGCTCCAGAGTGGCCAGGCGTCGCAGCAGCGCCTGCCCGGAAGCGCCGGAGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCATGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGTTCGagcagcag                                                                               >  1:10442685/4‑150 (MQ=255)
                                                                   cTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCTGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagcag                                                                               <  1:10689374/150‑1 (MQ=255)
                                                                      ggCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGcagcagcag                                                                            >  1:5837888/1‑150 (MQ=255)
                                                                      ggCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGcagcagcag                                                                            <  2:5019182/150‑1 (MQ=255)
                                                                      ggCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGcagcagcag                                                                            <  1:4122852/150‑1 (MQ=255)
                                                                       gcgcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCATCAGCAgg                                                                           <  2:1256449/150‑1 (MQ=255)
                                                                       gcgcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAgg                                                                           <  1:7366290/150‑1 (MQ=255)
                                                                          cccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                        <  2:1841810/150‑1 (MQ=255)
                                                                          cccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAGTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                        <  1:922921/150‑1 (MQ=255)
                                                                             tgctgcTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTcc                                                                     <  2:2832912/150‑1 (MQ=255)
                                                                              gctgctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTccc                                                                    >  2:9515812/1‑150 (MQ=255)
                                                                               ctgctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCg                                                                   <  2:23067/150‑1 (MQ=255)
                                                                               ctgctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATGGTCCTTCAGGTGGTTGAGCAGCGCCTGCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCg                                                                   >  2:4601667/1‑150 (MQ=255)
                                                                                 gctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGt                                                                 <  2:5775474/150‑1 (MQ=255)
                                                                                 gctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGt                                                                 <  2:9596419/150‑1 (MQ=255)
                                                                                 gctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGt                                                                 >  1:9416729/1‑150 (MQ=255)
                                                                                    ccAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCgg                                                              >  1:4390915/1‑150 (MQ=255)
                                                                                    ccAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCgg                                                              <  1:5416068/150‑1 (MQ=255)
                                                                                       gagTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCa                                                           <  1:5912041/150‑1 (MQ=255)
                                                                                         gTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCAcc                                                         >  1:10281786/1‑150 (MQ=255)
                                                                                          tGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCt                                                        >  1:1578867/1‑150 (MQ=255)
                                                                                          tGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCt                                                        >  1:7721289/1‑150 (MQ=255)
                                                                                                 gcggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggc                                                 <  2:4507415/150‑1 (MQ=255)
                                                                                                  cggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcg                                                <  2:10702583/150‑1 (MQ=255)
                                                                                                    gcgcAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcggc                                              <  1:2354307/150‑1 (MQ=255)
                                                                                                     cgcAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGgcggcg                                             <  2:5950878/150‑1 (MQ=255)
                                                                                                      gcagcaGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGa                                            >  1:4901452/1‑150 (MQ=255)
                                                                                                      gcagcaGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGa                                            >  2:8714919/1‑150 (MQ=255)
                                                                                                       cagcagCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGaa                                           >  2:4684840/1‑150 (MQ=255)
                                                                                                       cagcagCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGaa                                           >  2:7226817/1‑150 (MQ=255)
                                                                                                        agcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAAc                                          <  1:3345219/150‑1 (MQ=255)
                                                                                                         gcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACa                                         <  2:1437121/150‑1 (MQ=255)
                                                                                                         gcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTGCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCAACCTGGAAGGCGGCGAACa                                         >  2:6678562/1‑150 (MQ=255)
                                                                                                            gcgcCTGCCCGGAAGCGCCGGAGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTGCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCTGCAGGCCTTCCCGGTCGGCCACCTGGAGGGCGGCGAACAgcg                                      >  1:7259229/1‑150 (MQ=255)
                                                                                                             cgcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgcgc                                     >  2:7153242/1‑150 (MQ=255)
                                                                                                              gcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGTATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGcgcg                                    <  1:192448/150‑1 (MQ=255)
                                                                                                              gcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGcgcg                                    >  1:766967/1‑150 (MQ=255)
                                                                                                              gcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCGTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGcgcg                                    >  2:6324584/1‑150 (MQ=255)
                                                                                                               ccTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCgg                                   <  2:10105571/150‑1 (MQ=255)
                                                                                                                 tGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGcc                                 >  2:7763839/1‑150 (MQ=255)
                                                                                                                  gCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCct                                <  1:8972923/150‑1 (MQ=255)
                                                                                                                   cccGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctg                               >  1:2547354/1‑150 (MQ=255)
                                                                                                                    ccGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgc                              <  2:10146187/150‑1 (MQ=255)
                                                                                                                    ccGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgc                              >  1:7098217/1‑150 (MQ=255)
                                                                                                                    ccGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgc                              >  1:10580780/1‑150 (MQ=255)
                                                                                                                              cGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCgg                    >  1:649810/1‑150 (MQ=255)
                                                                                                                               ggTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGt                   >  1:5667249/1‑150 (MQ=255)
                                                                                                                               ggTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGt                   >  1:9746203/1‑150 (MQ=255)
                                                                                                                                gTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAGTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTc                  <  1:2818738/150‑1 (MQ=255)
                                                                                                                                 tGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCa                 <  1:7504257/150‑1 (MQ=255)
                                                                                                                                  gAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCaa                >  1:5280701/1‑150 (MQ=255)
                                                                                                                                  gAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCaa                <  1:9503916/150‑1 (MQ=255)
                                                                                                                                   aaCAGCGCCCGGTCCTCGCCGAGGATCTACGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGTCCTTCCCGGTCGGCCACCTGGAAGGCGTCGAACAGCGCGGCCTGCTGGCTATCGGTCAAc               <  1:6781300/150‑1 (MQ=255)
                                                                                                                                       gcgcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCa           >  2:8299134/1‑150 (MQ=255)
                                                                                                                                       gcgcCCGGTCCTCGCCGAGGATCTCCGGATTGTACTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTCGAAGGCGGCGAACAGAGCGGCCTGCTGGCTATCGGTCAACTCCa           >  2:3369154/1‑150 (MQ=255)
                                                                                                                                        cgcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCag          <  2:10304459/150‑1 (MQ=255)
                                                                                                                                          cccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagga        <  1:6094430/150‑1 (MQ=255)
                                                                                                                                          cccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagga        <  1:8423793/150‑1 (MQ=255)
                                                                                                                                            cgctCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGAGGTTGAGCACCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg      <  2:7588056/147‑1 (MQ=255)
                                                                                                                                            cGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg      <  1:9990421/150‑1 (MQ=255)
                                                                                                                                              gTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGtt    >  2:8361453/1‑150 (MQ=255)
                                                                                                                                              gTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGtt    >  1:9712665/1‑150 (MQ=255)
                                                                                                                                               tCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGACCAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTg   >  1:8555426/1‑150 (MQ=255)
                                                                                                                                                ccTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTgg  >  2:6477309/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGG  >  NC_008463/4997639‑4997932

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