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% 1960.5 / NA 496K173E (AAG→GAG) PA14_RS22840DUF853 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (235/261);  total (235/261)

Only 100 of 496 total aligned reads displayed.

GCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGG  >  NC_008463/4997639‑4997932
                                                                                                                                                   |                                                                                                                                                  
 cccccGCCCGTCCGGCTGGGGGCGGAGGATGTCTTCCAGGTGCAGGGCAGGCTCGCCGACGAAGGCCGCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGACGCGCCGGTGAACAGCGCCCGGTCCTCGcg                                                                                                                                                 >  2:8373975/1‑149 (MQ=255)
   cccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCCGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGa                                                                                                                                               >  2:6273067/1‑150 (MQ=255)
     cgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCTGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAgg                                                                                                                                             >  2:367453/1‑150 (MQ=255)
     cgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAgg                                                                                                                                             <  1:4784327/150‑1 (MQ=255)
      gcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGa                                                                                                                                            >  1:7833444/1‑150 (MQ=255)
       cccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAt                                                                                                                                           <  1:6358640/150‑1 (MQ=255)
       cccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCACGGTCCTCGCCGAGGAt                                                                                                                                           >  1:8759327/1‑150 (MQ=255)
         cgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAtct                                                                                                                                         >  1:293656/1‑150 (MQ=255)
         cgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAtct                                                                                                                                         <  2:6959882/150‑1 (MQ=255)
           tCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTcc                                                                                                                                       <  2:2866904/150‑1 (MQ=255)
            ccGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCg                                                                                                                                      <  2:5492179/150‑1 (MQ=255)
            ccGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCg                                                                                                                                      <  1:4458830/150‑1 (MQ=255)
              ggCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGa                                                                                                                                    <  2:1327536/150‑1 (MQ=255)
              ggCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGa                                                                                                                                    <  1:9394525/150‑1 (MQ=255)
               gCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAt                                                                                                                                   >  1:1388591/1‑150 (MQ=255)
                 tCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGACCAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTg                                                                                                                                 >  1:3098654/1‑150 (MQ=255)
                   ggggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTc                                                                                                                               >  2:6674113/1‑150 (MQ=255)
                     ggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCt                                                                                                                             <  1:6565524/150‑1 (MQ=255)
                        ggaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCa                                                                                                                          <  1:10660123/150‑1 (MQ=255)
                         gaggagGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAg                                                                                                                         >  2:1650798/1‑150 (MQ=255)
                         gaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAg                                                                                                                         >  2:1214511/1‑150 (MQ=255)
                          aggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAgg                                                                                                                        <  1:10175743/150‑1 (MQ=255)
                          aggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAgg                                                                                                                        >  2:2295220/1‑150 (MQ=255)
                               gTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGtt                                                                                                                   >  1:4933775/1‑150 (MQ=255)
                                tCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCGCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTg                                                                                                                  <  2:8633657/150‑1 (MQ=255)
                                 cTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGa                                                                                                                 <  1:4705995/150‑1 (MQ=255)
                                  ttCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGag                                                                                                                >  2:4203743/1‑150 (MQ=255)
                                   tCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCNAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagc                                                                                                               >  1:5669152/1‑150 (MQ=255)
                                    ccAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagca                                                                                                              >  2:9533671/1‑150 (MQ=255)
                                     cAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagcag                                                                                                             >  1:5342392/1‑150 (MQ=255)
                                      aGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagcagc                                                                                                            <  1:736740/150‑1 (MQ=255)
                                       gTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAgcg                                                                                                           <  2:10376638/150‑1 (MQ=255)
                                           cagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcctt                                                                                                       >  1:9072118/1‑150 (MQ=255)
                                            agggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttc                                                                                                      <  1:9692138/150‑1 (MQ=255)
                                             gggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttca                                                                                                     <  2:3090675/150‑1 (MQ=255)
                                             gcgcaggCTCGCCGAAGAAGGCCTCGGCGCCCTCCTGCTCCAGAGTGCCCAGGCGGCGCAGCAGCGCCTGCCCGGACGCGCCGGTGAACAGCGCCCGGTCCACGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttca                                                                                                     <  2:9165839/148‑1 (MQ=255)
                                                     tCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttc                                                                                             <  2:5074934/150‑1 (MQ=255)
                                                     tCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGCTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttc                                                                                             >  1:2657347/1‑150 (MQ=255)
                                                      cGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttca                                                                                            >  1:252627/1‑150 (MQ=255)
                                                      cGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttca                                                                                            >  2:2365457/1‑150 (MQ=255)
                                                        ccGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcagg                                                                                          >  2:7696855/1‑150 (MQ=255)
                                                         ctaagatgGCCTCGGCGCCCTGCTGCTCCAGTGTGGCCAGGGGGCGCAGCAGCGCCTGCCGGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGTATCTCCGGATTGTCCTTCTGGTGGTTGAGAAGCGccttctggtccttcaggg                                                                                         <  2:1821101/148‑2 (MQ=255)
                                                          gaagaaGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCcttcaggtccttcaggtc                                                                                        >  1:885380/1‑150 (MQ=255)
                                                          gaagaaGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCcttcaggtccttcaggtc                                                                                        <  1:10037543/150‑1 (MQ=255)
                                                          gaagaaGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCcttcaggtccttcaggtc                                                                                        >  2:7415760/1‑150 (MQ=255)
                                                          gaagaaGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATGGTCCTTCAGGTGGTTGAGCAGCGCcttcaggtccttcaggtc                                                                                        <  2:7511844/150‑1 (MQ=255)
                                                           aagaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCg                                                                                       >  1:3621435/1‑150 (MQ=255)
                                                           aagaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCg                                                                                       >  1:3399740/1‑150 (MQ=255)
                                                           aagaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTGCAGGTCCTTCAGGTCg                                                                                       >  2:5728582/1‑150 (MQ=255)
                                                            agaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGa                                                                                      <  2:3470372/150‑1 (MQ=255)
                                                            agaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGa                                                                                      <  1:10477371/150‑1 (MQ=255)
                                                                 gCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagc                                                                                 <  1:7232088/150‑1 (MQ=255)
                                                                   cTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagcag                                                                               <  2:7882511/150‑1 (MQ=255)
                                                                    tCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagcagc                                                                              <  1:7570054/150‑1 (MQ=255)
                                                                     cGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAgcagcagca                                                                             <  2:5771308/150‑1 (MQ=255)
                                                                       gcgcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAgg                                                                           >  2:8239468/1‑150 (MQ=255)
                                                                          cccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                        >  1:2049474/1‑150 (MQ=255)
                                                                          cccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                        <  2:9335382/150‑1 (MQ=255)
                                                                          cccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                        <  1:9143553/150‑1 (MQ=255)
                                                                          cccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                        <  2:4679392/150‑1 (MQ=255)
                                                                           ccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGACGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCGTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCtt                                                                       >  2:3867461/1‑150 (MQ=255)
                                                                            ctgctgCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTc                                                                      <  1:1540558/150‑1 (MQ=255)
                                                                               ctgctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCg                                                                   <  2:8280786/150‑1 (MQ=255)
                                                                                  ctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTc                                                                >  1:9836955/1‑150 (MQ=255)
                                                                                   tCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCg                                                               >  1:5212961/1‑150 (MQ=255)
                                                                                      agagTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGcc                                                            <  1:2359393/150‑1 (MQ=255)
                                                                                       gagTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCa                                                           >  1:6728006/1‑150 (MQ=255)
                                                                                       gagTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCGTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCa                                                           >  1:5803482/1‑150 (MQ=255)
                                                                                           ggCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTg                                                       <  1:5545195/150‑1 (MQ=255)
                                                                                              cAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGaa                                                    >  2:2582447/1‑150 (MQ=255)
                                                                                                    gcgcAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcggc                                              >  2:5832352/1‑150 (MQ=255)
                                                                                                    gcgcAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcggc                                              <  1:8244327/150‑1 (MQ=255)
                                                                                                        agcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAGTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAAc                                          <  2:2028491/150‑1 (MQ=255)
                                                                                                           agcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgc                                       >  2:2672480/1‑150 (MQ=255)
                                                                                                           agcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgc                                       <  2:4047659/150‑1 (MQ=255)
                                                                                                             cgcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgcgc                                     >  1:2765005/1‑150 (MQ=255)
                                                                                                               ccagCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTAAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCAGGAAGGCGGCGAACAGCGCgg                                   >  1:6903742/4‑150 (MQ=255)
                                                                                                                cTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGc                                  <  2:3302549/150‑1 (MQ=255)
                                                                                                                 tGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGcc                                 >  1:10575712/1‑150 (MQ=255)
                                                                                                                     cGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgct                             <  2:585183/150‑1 (MQ=255)
                                                                                                                     cGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgct                             <  2:5977055/150‑1 (MQ=255)
                                                                                                                     cGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgct                             <  2:7219223/150‑1 (MQ=255)
                                                                                                                       gAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTgg                           <  1:3226769/150‑1 (MQ=255)
                                                                                                                        aaGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGc                          <  1:9577493/150‑1 (MQ=255)
                                                                                                                         aGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGCCCTGCTGGCt                         >  2:8970112/1‑150 (MQ=255)
                                                                                                                           cgcCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTAt                       <  1:1835967/150‑1 (MQ=255)
                                                                                                                           cgcCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTAt                       <  2:1478695/150‑1 (MQ=255)
                                                                                                                              cGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTCTCgg                    >  2:5242191/1‑150 (MQ=255)
                                                                                                                               ggTGACCAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGt                   <  1:6099578/150‑1 (MQ=255)
                                                                                                                               ggTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGt                   <  2:6803944/150‑1 (MQ=255)
                                                                                                                                gTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTc                  <  1:8375831/150‑1 (MQ=255)
                                                                                                                                 tGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCa                 >  1:1929628/1‑150 (MQ=255)
                                                                                                                                   caCAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAAc               <  2:6110278/149‑1 (MQ=255)
                                                                                                                                         gcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagg         >  1:4045453/1‑150 (MQ=255)
                                                                                                                                          cccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagga        >  2:10473113/1‑150 (MQ=255)
                                                                                                                                            cGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg      >  1:4203412/1‑150 (MQ=255)
                                                                                                                                              gTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGtt    <  2:4417355/150‑1 (MQ=255)
                                                                                                                                               tCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTg   <  2:5344166/150‑1 (MQ=255)
                                                                                                                                               tCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTg   <  2:9755515/150‑1 (MQ=255)
                                                                                                                                                ccTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTgg  >  2:6548550/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGG  >  NC_008463/4997639‑4997932

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