Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 126,115 A→G 8.9% intergenic (+25/+25) PP_0120 → / ← thrB zinc ABC transporter/homoserine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947126,1150AG 8.9% 249.0 / 11.6 90intergenic (+25/+25)PP_0120/thrBzinc ABC transporter/homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base G (6/2);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

GTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGC  >  NC_002947/125974‑126258
                                                                                                                                             |                                                                                                                                               
gTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcgc                                                                                                                                          >  5:180956/1‑149 (MQ=255)
   cggcTTGCCGAGCTTGATGCCTTTGGTTGCACCGATCCGGTGGATGGCAAATGGTATGAGCGCTTGCTGGAGAATCTGGGCGGTGACCTGGAAGTCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCga                                                                                                                                        <  2:56687/145‑1 (MQ=255)
   cGCCTTGCCGATCTGTTTGACTTGTGGGGGACCGATCCGGTGGATGTCAATGTGTATTAGAGCTTGCTGGAGAAGCTGGTCGGTGACCTGGCAGGCTGCCTGTAGCATTTTTAAGGCTGTGCCGGCCCCTTCGCGGGCGTGCCCGCgaa                                                                                                                                       <  8:330506/149‑1 (MQ=255)
    gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcgcg                                                                                                                                         >  1:196570/1‑146 (MQ=255)
     ccTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaaga                                                                                                                                     <  3:126681/149‑1 (MQ=255)
      cTTGCTGAGCTGGATGCCTTGGGTGGCAGCGATGCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaagaa                                                                                                                                    <  8:151533/149‑1 (MQ=255)
      cTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaagaa                                                                                                                                    <  7:224745/149‑1 (MQ=255)
      cTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaagaa                                                                                                                                    <  3:304081/149‑1 (MQ=255)
       ttGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGaagaag                                                                                                                                   >  3:183627/1‑149 (MQ=255)
         gCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGc                                                                                                                                  <  4:364022/148‑1 (MQ=255)
          ccGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAg                                                                                                                                <  7:193833/149‑1 (MQ=255)
          ccGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAg                                                                                                                                <  8:114915/149‑1 (MQ=255)
          ccGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAg                                                                                                                                <  4:45557/149‑1 (MQ=255)
           cGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGa                                                                                                                               >  3:288397/1‑149 (MQ=255)
            gAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGAGGTGACCTGGCAGGCTGCATGGAGAAGTTGTAAGCCTGTGCCGGCCCCTTCTCGGGCACGCCCGCga                                                                                                                                        >  3:93091/1‑139 (MQ=255)
                tGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCgg                                                                                                                          >  7:9420/1‑149 (MQ=255)
                 ggATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGt                                                                                                                         >  1:4878/1‑149 (MQ=255)
                       cTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATTAGAGCTTGCTGGAGAAGATGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCACTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGCTCTACAg                                                                                                                   >  1:13539/1‑149 (MQ=255)
                       cTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAg                                                                                                                   >  1:3787/1‑149 (MQ=255)
                       cTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAg                                                                                                                   >  8:244162/1‑149 (MQ=255)
                       cTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAg                                                                                                                   >  7:316408/1‑149 (MQ=255)
                       cTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAg                                                                                                                   >  7:275311/1‑149 (MQ=255)
                        ttGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGAGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAgc                                                                                                                  >  7:204154/1‑149 (MQ=255)
                          ggggggCACCGAGCCGCGGGCTGGCAAAGGGGACCTGCGCGTGCCGGAGCAGCTGGGCGGTGCCCTGGCAGGCTGCCGGGAGCAGTTGGAAGCCGGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAgcgc                                                                                                                <  8:242713/149‑1 (MQ=255)
                          gggTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAgc                                                                                                                  <  4:243323/147‑1 (MQ=255)
                           ggTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGcgcg                                                                                                               <  4:291591/149‑1 (MQ=255)
                               gCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAAc                                                                                                            <  2:4878/148‑1 (MQ=255)
                                cACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACgg                                                                                                          <  2:3787/149‑1 (MQ=255)
                                cACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACgg                                                                                                          >  5:206658/1‑149 (MQ=255)
                                 aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCGTGCCCGCgaagg                                                                                                                                     >  7:195653/1‑120 (MQ=255)
                                 aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCGTGCCCGCgaagg                                                                                                                                     <  8:195653/121‑2 (MQ=255)
                                 aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACgg                                                                                                          >  6:34531/1‑148 (MQ=255)
                                 aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGAAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGc                                                                                                         >  5:291159/1‑149 (MQ=255)
                                  ccGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCa                                                                                                        <  4:288397/149‑1 (MQ=255)
                                          gTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAAt                                                                                                <  8:153526/149‑1 (MQ=255)
                                          gTGGATGGCAAAGGGGATCAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGGAGTTGTAAGCCTGTGCCGGCCCCTTCGCTGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAAt                                                                                                <  6:295520/149‑1 (MQ=255)
                                                  cAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGaa                                                                                                             >  4:150093/1‑128 (MQ=255)
                                                  cAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGaa                                                                                                             <  3:150093/128‑1 (MQ=255)
                                                   aaaGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACctgc                                                                                       >  3:321994/1‑149 (MQ=255)
                                                    aatggTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCTGGCTGCCTGTAGCTGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCTTAGAAGCAGACGCGGTCTACAGCGCGAACGTCAGGTGCAATGCCACctgct                                                                                      <  2:13539/146‑1 (MQ=255)
                                                     aGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACctgctg                                                                                     <  4:321994/149‑1 (MQ=255)
                                                       ggTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGgc                                                                                   <  2:342710/149‑1 (MQ=255)
                                                           tGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTg                                                                               >  8:62785/1‑149 (MQ=255)
                                                           tGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTg                                                                               >  1:137243/1‑149 (MQ=255)
                                                             cgcggTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTCGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCGGGCGCTgc                                                                              <  5:295590/147‑1 (MQ=255)
                                                               cgcTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                            >  1:230568/1‑148 (MQ=255)
                                                                  ttGCTGGAGAAGCTGGGCGGTGACCTGGCAGTCTGGCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGa                                                                        <  3:310246/149‑1 (MQ=255)
                                                                    gCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGAc                                                                       <  4:183627/148‑1 (MQ=255)
                                                                    gCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGAc                                                                       <  8:275311/148‑1 (MQ=255)
                                                                      tGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAg                                                                                                       <  7:276085/114‑1 (MQ=255)
                                                                      tGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAg                                                                                                       >  8:276085/1‑114 (MQ=255)
                                                                       ggAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaa                                                                                                                                       >  8:219469/1‑81 (MQ=255)
                                                                       ggAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaa                                                                                                                                       <  7:219469/81‑1 (MQ=255)
                                                                       ggAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCAc                                                                                           >  6:302942/1‑125 (MQ=255)
                                                                       ggAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCAc                                                                                           <  5:302942/125‑1 (MQ=255)
                                                                          gaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTc                                                                >  8:175725/1‑149 (MQ=255)
                                                                             gCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGaa                                                             <  2:196570/149‑1 (MQ=255)
                                                                                gggCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActc                                                          >  2:124377/1‑149 (MQ=255)
                                                                                gggCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCCGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCACCCGCCCTTCGAActc                                                          >  8:284703/1‑149 (MQ=255)
                                                                                 ggCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActc                                                          <  7:38853/148‑1 (MQ=255)
                                                                                   cGGTGACCTGGCAGGCTGCCTGTAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCTAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgct                                                       <  4:193247/149‑1 (MQ=255)
                                                                                    ggTGACCTGGGATGCTGGCGTGTCGAGTTGTAATCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCCGTCTACAGCGCGAACTGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                      <  8:9420/149‑1 (MQ=255)
                                                                                    ggTGACCTGGCAGGCTGCCTTGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                      <  3:56764/149‑1 (MQ=255)
                                                                                    ggTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                      <  7:175725/149‑1 (MQ=255)
                                                                                        aCCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGAt                                                  <  7:62785/149‑1 (MQ=255)
                                                                                          cTGGCAGGCTGCCTGGAGCTGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGGGAACGGCAGGGGCAATGCCAGCTGGTGGCGCTGCGCCTGTCGCGCGTCGAGCTGGCTCGGAtgg                                                >  8:239085/1‑147 (MQ=255)
                                                                                               aGGCTGCCTGGAGCAGTTCTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAACAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCCTGGa                                            >  5:220007/1‑148 (MQ=255)
                                                                                                ggCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGt                                                                                                     >  6:444/1‑90 (MQ=255)
                                                                                                     ccTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatcatc                                       <  6:180956/147‑1 (MQ=255)
                                                                                                     ccTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATcac                                     <  5:34531/149‑1 (MQ=255)
                                                                                                     ccTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATcac                                     <  2:230568/149‑1 (MQ=255)
                                                                                                      cTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATcaca                                    <  8:316408/149‑1 (MQ=255)
                                                                                                          aGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcatcca                                >  5:80043/1‑149 (MQ=255)
                                                                                                             aGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACatccatcc                             >  6:271649/1‑149 (MQ=255)
                                                                                                               ttGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCTGACGCGGTCTAc                                                                                                                     >  4:147834/1‑59 (MQ=255)
                                                                                                               ttGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatcat                                        <  4:189533/136‑1 (MQ=255)
                                                                                                               ttGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatcat                                        >  3:189533/1‑136 (MQ=255)
                                                                                                               ttGTAAGCCTGTGCCGGCCCCTTCGCGGGCACCCCACCGAAACACCATACTCCGTATACACCACGAACGGCACGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACATCGAACTCGCTAGGATCGTTGATCATCACATCCATCCCg                           >  8:202400/1‑149 (MQ=255)
                                                                                                                 gTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGc                         >  7:23837/1‑149 (MQ=255)
                                                                                                                  tAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGc                         <  7:244162/148‑1 (MQ=255)
                                                                                                                     gCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAAc                     >  5:114650/1‑149 (MQ=255)
                                                                                                                               gCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCg           <  5:271649/149‑1 (MQ=255)
                                                                                                                                 ttcTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGc          >  7:326523/3‑148 (MQ=255)
                                                                                                                                 cccTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGc          >  2:301021/1‑148 (MQ=255)
                                                                                                                                   cTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAAt       >  5:340561/1‑149 (MQ=255)
                                                                                                                                    ttCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAAt       >  7:344567/1‑148 (MQ=255)
                                                                                                                                     tCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAAt       >  4:347930/1‑147 (MQ=255)
                                                                                                                                     tCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCAGGCGAACGCCTCGGCCGCAATcc     >  4:353444/1‑148 (MQ=255)
                                                                                                                                      cgcgGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATc      <  1:301021/147‑1 (MQ=255)
                                                                                                                                        cgGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGgc  <  1:99881/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGC  >  NC_002947/125974‑126258

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