Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,341,704 Δ1 bp pseudogene (664/957 nt) PA14_RS10890 → DUF58 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,341,7020C.97.7% 2267.8 / 38.3 485pseudogene (662/957 nt)PA14_RS10890DUF58 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base . (245/227);  total (252/233)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 497 total aligned reads displayed.

GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848
                                                                                                                                                        |                                                                                                                                                  
  gCTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAcca                                                                                                                                                  >  2:6914903/1‑150 (MQ=255)
   cTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccag                                                                                                                                                 >  1:7892844/1‑149 (MQ=255)
   cTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccag                                                                                                                                                 >  1:1435225/1‑149 (MQ=255)
     gggTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagcc                                                                                                                                               <  2:4052523/151‑5 (MQ=255)
       gTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCgg                                                                                                                                            >  1:756325/1‑151 (MQ=255)
       gTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCgg                                                                                                                                            <  1:3256862/151‑1 (MQ=255)
        tCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCggg                                                                                                                                           >  2:1510940/1‑151 (MQ=255)
        tCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCggg                                                                                                                                           <  2:6087413/151‑1 (MQ=255)
         cGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGc                                                                                                                                          <  2:3910157/151‑1 (MQ=255)
         cGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGc                                                                                                                                          >  2:9493635/1‑151 (MQ=255)
          gACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCg                                                                                                                                         >  2:4485357/1‑151 (MQ=255)
          gACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCg                                                                                                                                         >  1:7461454/1‑151 (MQ=255)
             ccagcgccAGCAGGCGCTGGACTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGAc                                                                                                                                       <  1:4700347/145‑1 (MQ=255)
               tgcgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTc                                                                                                                                     >  1:8477676/5‑151 (MQ=255)
                  gccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCCGGGACTGCGCGCCTACCAGCCGGGCGACTcgcg                                                                                                                                  >  1:6816722/2‑151 (MQ=255)
                  gccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTcgcg                                                                                                                                  >  2:6553140/2‑151 (MQ=255)
                   ccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcgg                                                                                                                                 >  2:488990/1‑151 (MQ=255)
                       cagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggc                                                                                                                             <  2:5192534/151‑1 (MQ=255)
                        agGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggcg                                                                                                                            <  1:7336027/151‑1 (MQ=255)
                            gcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTg                                                                                                                        >  2:1822116/1‑151 (MQ=255)
                             cTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGc                                                                                                                       <  1:1004998/151‑1 (MQ=255)
                              tGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCa                                                                                                                      <  1:6298708/151‑1 (MQ=255)
                              tGGTCGACCCTCGGCCGCTGGCCGGCGCGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCa                                                                                                                      <  2:52246/151‑1 (MQ=255)
                              aggTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCa                                                                                                                      <  1:3529828/150‑1 (MQ=255)
                                 tCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTg                                                                                                                   >  1:1769417/1‑151 (MQ=255)
                                    aCCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGaa                                                                                                                >  2:33022/1‑151 (MQ=255)
                                     cccTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAg                                                                                                               <  1:6216884/151‑1 (MQ=255)
                                     cccTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAg                                                                                                               <  2:8474257/151‑1 (MQ=255)
                                        tCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGCGGATGCCGTGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCg                                                                                                            >  1:1227460/1‑151 (MQ=255)
                                         cGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGt                                                                                                           >  1:498775/1‑151 (MQ=255)
                                         cGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGt                                                                                                           >  1:6716679/1‑151 (MQ=255)
                                         cGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGt                                                                                                           <  2:8111305/151‑1 (MQ=255)
                                          ggCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtt                                                                                                          <  2:298077/151‑1 (MQ=255)
                                             cGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTc                                                                                                       <  2:5531867/151‑1 (MQ=255)
                                             cGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTAACAGGGACGGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTc                                                                                                       >  1:9547495/1‑151 (MQ=255)
                                               cTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTccc                                                                                                     >  2:2544263/1‑151 (MQ=255)
                                                 ggCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGt                                                                                                   <  2:9766224/151‑1 (MQ=255)
                                                 ggCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGt                                                                                                   <  2:8660552/151‑1 (MQ=255)
                                                    cGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGc                                                                                                >  1:4842647/1‑151 (MQ=255)
                                                     ggCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGCCTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGcc                                                                                               >  1:2234925/1‑151 (MQ=255)
                                                     ggCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGcc                                                                                               >  1:5292703/1‑151 (MQ=255)
                                                     ggCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGcc                                                                                               >  1:2533137/1‑151 (MQ=255)
                                                        cgagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGACGGCGTTGTCCCGTGGCCAg                                                                                             >  2:5107341/5‑151 (MQ=255)
                                                         gagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAgg                                                                                            <  2:7909550/148‑1 (MQ=255)
                                                         gagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGAATACCAGGGACTGCGCGCCGACCAGCCGGGCGACTCGCGGCGGCGCCTGAACTGGAAGGCGTTTTCCCGGGGCCAgg                                                                                            >  2:5898814/4‑151 (MQ=255)
                                                         gagccGGTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAgg                                                                                            <  2:965319/148‑1 (MQ=255)
                                                            ccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgc                                                                                         <  1:6040438/151‑1 (MQ=255)
                                                             cGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgct                                                                                        >  2:6400748/1‑151 (MQ=255)
                                                              gTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACAAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctg                                                                                       >  1:5123200/1‑151 (MQ=255)
                                                                  cTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGCCTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTgg                                                                                   <  2:4280550/151‑1 (MQ=255)
                                                                  cTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTgg                                                                                   <  2:7362810/151‑1 (MQ=255)
                                                                   tCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGt                                                                                  <  1:10840740/151‑1 (MQ=255)
                                                                     gcgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGa                                                                                <  2:770860/151‑1 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGaa                                                                               >  2:2031171/1‑151 (MQ=255)
                                                                        ccGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAgg                                                                             >  1:1623104/1‑151 (MQ=255)
                                                                          gggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGAc                                                                           >  2:8258503/1‑151 (MQ=255)
                                                                                ggaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgcc                                                                     <  1:802215/151‑1 (MQ=255)
                                                                                    gaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCgcgc                                                                 >  1:3926664/1‑151 (MQ=255)
                                                                                    gaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCgcgc                                                                 >  2:2208588/1‑151 (MQ=255)
                                                                                     aTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCt                                                                >  2:9027088/1‑151 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                               >  1:2790742/1‑151 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                               <  1:5443031/151‑1 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            <  2:9274790/151‑1 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            <  1:3196328/151‑1 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            <  1:417618/151‑1 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            <  2:1326231/151‑1 (MQ=255)
                                                                                           aGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggcgg                                                          <  1:4343680/151‑1 (MQ=255)
                                                                                            ggACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCggg                                                         <  2:10550444/151‑1 (MQ=255)
                                                                                             gACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCgggg                                                        <  1:8107676/151‑1 (MQ=255)
                                                                                                    ggCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCt                                                 <  2:8699110/151‑1 (MQ=255)
                                                                                                            cTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGAGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGAATTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTgg                                         >  1:291675/1‑151 (MQ=255)
                                                                                                             tGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                        <  2:2825785/151‑1 (MQ=255)
                                                                                                              ggCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGaa                                       >  1:5741775/1‑151 (MQ=255)
                                                                                                               gCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAt                                      <  1:269562/151‑1 (MQ=255)
                                                                                                               gCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAt                                      >  1:1283966/1‑151 (MQ=255)
                                                                                                                cGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGACCTGTGGCTGGACtt                                     >  2:2388360/1‑151 (MQ=255)
                                                                                                                cGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAtt                                     <  1:10576678/151‑1 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGACCTGTGGCTGGAATTc                                    >  1:8331583/1‑151 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                    >  1:220970/1‑151 (MQ=255)
                                                                                                                  aGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGACGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCg                                   >  2:5723689/1‑151 (MQ=255)
                                                                                                                  aGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCg                                   >  1:6422967/1‑151 (MQ=255)
                                                                                                                  aGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTGTCCCGTGGCCAGGGGCTGCTGGTGCAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCg                                   >  2:350712/1‑151 (MQ=255)
                                                                                                                   gCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGa                                  <  2:10611993/151‑1 (MQ=255)
                                                                                                                   gCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGCGGAAGACCTGTGGCTGGAATTCGa                                  <  1:7032733/151‑1 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGACCTGTGGCTGGAATTCGAgc                                >  2:720214/1‑151 (MQ=255)
                                                                                                                      cGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgcg                               >  2:10029717/1‑151 (MQ=255)
                                                                                                                        ggCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGcc                             >  1:9677948/1‑151 (MQ=255)
                                                                                                                        ggCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGcc                             >  1:1541801/1‑151 (MQ=255)
                                                                                                                         gCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAATGCGTTTTCCCGTGGCCTGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGTATTCGTGCGCCt                            <  2:3117907/151‑1 (MQ=255)
                                                                                                                         gCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCt                            >  1:6527564/1‑151 (MQ=255)
                                                                                                                            tCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcgg                         >  2:7530371/1‑151 (MQ=255)
                                                                                                                              gacgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcg                       >  2:9085740/1‑151 (MQ=255)
                                                                                                                              gacgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcg                       >  1:10248835/1‑151 (MQ=255)
                                                                                                                                cgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGACTTCGAGCGCCTCGGCggg                     >  2:4211231/1‑151 (MQ=255)
                                                                                                                                cgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCggg                     >  1:8989621/1‑151 (MQ=255)
                                                                                                                                cgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCggg                     <  1:7672419/151‑1 (MQ=255)
                                                                                                                                 gacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCgggg                    >  1:7142857/1‑151 (MQ=255)
                                                                                                                                      ccAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGGTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGCATTCGAGCGCCTCGGCGGGGAgctc               >  2:8899687/1‑151 (MQ=255)
                                                                                                                                            aCTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGt         >  1:9448459/1‑151 (MQ=255)
                                                                                                                                                   gcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCGGGGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTgg  >  2:10211944/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848

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: