Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,598,833 G→A 8.6% intergenic (‑72/+59) dctB ← / ← bglX C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,598,8330GA 8.6% 231.0 / 7.8 81intergenic (‑72/+59)dctB/bglXC4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base A (2/5);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

CCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAA  >  NC_002947/1598685‑1598981
                                                                                                                                                    |                                                                                                                                                     
ccGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTg                                                                                                                                                       <  4:234751/149‑1 (MQ=255)
ccGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTg                                                                                                                                                       <  6:94562/149‑1 (MQ=255)
 cGGCGCTGAGTGCGGTACCCAGCAGGATCAGCAGGGTGATCAGCCCCAGACGGAGGGCGCGAAAGGAAAAGGGCATCGGGTGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTTCTCTgc                                                                                                                                                      <  6:291175/149‑1 (MQ=255)
 cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgc                                                                                                                                                      <  6:129202/149‑1 (MQ=255)
  ggCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgca                                                                                                                                                     <  5:54779/149‑1 (MQ=255)
      cTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAc                                                                                                                                                 <  7:14715/149‑1 (MQ=255)
         aGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAcc                                                                                                                                                <  1:297961/147‑1 (MQ=255)
         aGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTg                                                                                                                                              >  2:77387/1‑149 (MQ=255)
         aGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTg                                                                                                                                              >  2:383616/1‑149 (MQ=255)
          gTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGt                                                                                                                                             <  7:191183/149‑1 (MQ=255)
             cGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                           <  2:30783/148‑1 (MQ=38)
              ggTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGgc                                                                                                                                         >  5:324338/1‑149 (MQ=37)
              ggTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGgc                                                                                                                                         <  3:276563/149‑1 (MQ=37)
                  cccAGCCCGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGGGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGCGCTGCAGAGTACCTGTAGGcgcg                                                                                                                                      <  8:361722/148‑1 (MQ=37)
                   ccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGgcgacgt                                                                                                                                    >  2:341938/1‑145 (MQ=255)
                    cAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGtt                                                                                                                                   >  7:48111/1‑149 (MQ=255)
                     aGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAgg                                                                                                                                          >  6:31208/1‑141 (MQ=255)
                     aGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAgg                                                                                                                                          <  5:31208/141‑1 (MQ=255)
                         cGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTAcc                                                                                                                              <  6:324338/149‑1 (MQ=25)
                            tCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCcgc                                                                                                                           <  4:327653/149‑1 (MQ=255)
                              agcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCGCCTGTAGGCGCGGGTTTACCcgcg                                                                                                                          >  3:310160/1‑148 (MQ=255)
                              agcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCcgcg                                                                                                                          >  3:367257/1‑148 (MQ=255)
                               gcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGa                                                                                                                         >  2:349769/1‑148 (MQ=255)
                                  gggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTAcc                                                                                                                              <  8:65099/140‑1 (MQ=255)
                                  gggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTAcc                                                                                                                              >  7:65099/1‑140 (MQ=255)
                                   ggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaaa                                                                                                                      >  5:199645/1‑147 (MQ=255)
                                   ggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaaa                                                                                                                      <  6:199645/147‑1 (MQ=255)
                                   ggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAgg                                                                                                                    >  4:215723/1‑149 (MQ=255)
                                   ggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAgg                                                                                                                    >  8:213477/1‑149 (MQ=255)
                                     tGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGt                                                                                                                                    <  6:173482/131‑1 (MQ=255)
                                     tGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGt                                                                                                                                    >  5:173482/1‑131 (MQ=255)
                                         cAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGaa                                                                                                              >  5:105322/1‑149 (MQ=255)
                                           gCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACa                                                                                                            <  6:27951/149‑1 (MQ=255)
                                             cccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAgg                                                                                                          <  3:215723/149‑1 (MQ=255)
                                               cAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcg                                                                                                       >  4:79650/1‑149 (MQ=255)
                                                 gACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTAcac                                                                                                                             >  6:16508/1‑124 (MQ=255)
                                                 gACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTAcac                                                                                                                             <  5:16508/126‑3 (MQ=255)
                                                         gcgcGAAAGGTAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             <  3:294817/149‑1 (MQ=255)
                                                         gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             <  8:370058/149‑1 (MQ=255)
                                                             gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTc                                                                                          >  1:235844/1‑148 (MQ=255)
                                                               aaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                       <  3:79650/149‑1 (MQ=255)
                                                                    aaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAt                                                                                  <  4:279496/149‑1 (MQ=255)
                                                                       gggCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAGCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTc                                                                                >  5:296699/1‑148 (MQ=255)
                                                                       gggCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCa                                                                               <  1:77387/149‑1 (MQ=255)
                                                                         gCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                             <  2:235844/149‑1 (MQ=255)
                                                                             cGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGctc                                                                         >  3:217230/1‑149 (MQ=255)
                                                                              gggggCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGctc                                                                         <  4:217230/148‑1 (MQ=255)
                                                                               ggggCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctg                                                                       <  8:48111/149‑1 (MQ=255)
                                                                                         cGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCAcc                                                             <  1:383616/149‑1 (MQ=255)
                                                                                           gtggtgGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCgc                                                                                                                <  5:123575/97‑2 (MQ=255)
                                                                                           gtggtgGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCgc                                                                                                                >  6:123575/1‑96 (MQ=255)
                                                                                           gtggtgGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGc                                                           >  8:376504/1‑149 (MQ=255)
                                                                                            tggtggGATGTCGGGGCCTGCCTCTTCGCGTGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGcc                                                          <  4:310160/149‑1 (MQ=255)
                                                                                               tggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                                       <  4:367257/149‑1 (MQ=255)
                                                                                               tggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                                       <  7:213477/149‑1 (MQ=255)
                                                                                               tggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                                       >  3:65707/1‑149 (MQ=255)
                                                                                                  gATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGaa                                                    <  6:105322/149‑1 (MQ=255)
                                                                                                       cGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCCGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCTAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAg                                               >  3:13382/1‑149 (MQ=255)
                                                                                                              tgcctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCGTTTTCAGTCAAAGCAATTGAAAGCTCTGCTGGTGCCCCGCCTGGGAATCCAGCCCGAcc                                        >  5:40546/1‑149 (MQ=38)
                                                                                                              tgcctCTTCGCGGGTAAACCCGCGCCTAAAGCTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAACGCCAGAACAGGCGGCCATTTTAAGTCAAAACAATTCAAAGCTCTGCTGCTGCACCGACAGGGAATCCAGCCCTAcc                                        >  2:372053/1‑149 (MQ=38)
                                                                                                               gcctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCt                                       <  1:349769/149‑1 (MQ=255)
                                                                                                                 ctctTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCGTTTTCAGTCAAAGCAATTAAAAGCTCTGCAGCTGCACCGCCAGGGAATCCAGCCCTAccag                                      >  8:346054/1‑146 (MQ=255)
                                                                                                                  tctTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGaa                                                    >  6:298262/1‑133 (MQ=255)
                                                                                                                  tctTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGaa                                                    <  5:298262/133‑1 (MQ=255)
                                                                                                                      cgcgGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGAc                                   <  7:346054/146‑1 (MQ=255)
                                                                                                                      cgcgGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACa                                  <  4:65707/147‑1 (MQ=255)
                                                                                                                         gtgTAAACCCGTGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGaa                             <  1:341938/147‑1 (MQ=255)
                                                                                                                          ggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                             <  2:141453/100‑1 (MQ=255)
                                                                                                                          ggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                             >  1:141453/1‑100 (MQ=255)
                                                                                                                          ggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAAc                            >  2:271439/1‑149 (MQ=255)
                                                                                                                          ggTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGaa                             >  2:3125/1‑148 (MQ=25)
                                                                                                                                 ccGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTg                     >  4:335349/1‑149 (MQ=255)
                                                                                                                                     gcCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTc                 <  1:3125/149‑1 (MQ=37)
                                                                                                                                       cTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAtt                                                                                 <  1:382896/83‑1 (MQ=255)
                                                                                                                                       cTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAtt                                                                                 >  2:382896/1‑83 (MQ=255)
                                                                                                                                          cAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCaa            >  6:345477/1‑149 (MQ=255)
                                                                                                                                            ggTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCTGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAAcc          <  1:372053/149‑1 (MQ=255)
                                                                                                                                                  ctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                                       >  5:347372/1‑98 (MQ=255)
                                                                                                                                                  ctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                                       <  6:347372/98‑1 (MQ=255)
                                                                                                                                                   tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                             <  7:50226/75‑1 (MQ=255)
                                                                                                                                                   tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                             >  2:210610/1‑75 (MQ=255)
                                                                                                                                                   tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                             >  8:50226/1‑75 (MQ=255)
                                                                                                                                                   tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                             <  1:210610/75‑1 (MQ=255)
                                                                                                                                                   tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc    >  5:67877/1‑148 (MQ=255)
                                                                                                                                                   tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGCCCGGAACAGGCGGGCCTTTTTAAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTCGACATTGAACTCCCCTGGTTCGGCAACCCGCTg     >  1:95205/1‑148 (MQ=255)
                                                                                                                                                   tGCAGCACCTGTAGGCGCGGTATTACCAGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCGGCCCTACCTGGACAATGAACTCCCCTGGTTCGGCACCACGCTGc    >  6:257954/1‑148 (MQ=255)
                                                                                                                                                    gcagcaCCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAATGCAATTCACAGCTCTGCGGCTGCACCGCCTGGGCTTCCAGCCCTACCGGGCCAGTGAGCGGCCCGGGGTCGGCAACCCGGTGcgg  >  4:304810/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAA  >  NC_002947/1598685‑1598981

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