Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,677,312 C→A 10.0% intergenic (+43/+61) PP_4137 → / ← chrR outer membrane siderophore receptor/chromate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,677,3120CA10.0% 262.8 / 16.9 100intergenic (+43/+61)PP_4137/chrRouter membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand):  ref base C (21/63);  new base A (5/5);  total (32/68)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 101 total aligned reads displayed.

CAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGC  >  NC_002947/4677163‑4677458
                                                                                                                                                         |                                                                                                                                                  
cAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGc                                                                                                                                                    <  1:546905/150‑1 (MQ=255)
    gACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                                <  1:4779669/150‑1 (MQ=255)
     aCCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGCGGGG‑TTTACCC‑CGAATTCGTCAGT‑GAATTTTCCGAcc                                                                                                                                               >  1:437016/1‑149 (MQ=255)
     aCCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCCCGCATTCCGCATT‑TAAATTTCAGACg                                                                                                                                               >  1:436750/1‑150 (MQ=255)
      ccGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGc                                                                                                                                              <  2:1309216/150‑1 (MQ=255)
        gTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCAt                                                                                                                                            <  2:3665596/150‑1 (MQ=255)
            cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTcgc                                                                                                                                        <  1:4707577/150‑1 (MQ=255)
            cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTcgc                                                                                                                                        <  1:3763219/150‑1 (MQ=255)
            cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTcgc                                                                                                                                        <  1:3763054/150‑1 (MQ=255)
            cGCCGGCCGACCCCGACGCCGGCCAGGACCCGGCCCGTTCGACCCCGGGCGAAGGATGTGGCGTGTAGGCCGGGGTTT‑CCTACCGCTTCCGCTGACGGGTGGGCGC‑GGG‑TGTACCCGCGCATGCGTCATT‑GAATGTGCCGACGCATTcgc                                                                                                                                        <  2:1556516/150‑1 (MQ=255)
             gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTcgcg                                                                                                                                       <  2:1367917/150‑1 (MQ=37)
              ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                                      <  2:3402003/150‑1 (MQ=37)
              ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                                      <  2:443606/150‑1 (MQ=37)
              ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                                      <  2:2634065/150‑1 (MQ=37)
              ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                                      <  2:2052829/150‑1 (MQ=255)
              ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                                      <  2:3446300/150‑1 (MQ=255)
              ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                                      <  2:4565245/150‑1 (MQ=255)
              ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                                      <  2:5300967/150‑1 (MQ=255)
              ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGAGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                                      <  2:4215963/150‑1 (MQ=37)
                ggCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGt                                                                                                                                    >  2:5341912/1‑150 (MQ=255)
                   cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTAACCCCCAAAGCGTCAAT‑GAATTTTACGACGCATTAGCCGGGaaa                                                                                                                                 >  1:4484483/1‑150 (MQ=255)
                    tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAAc                                                                                                                                <  1:1996649/150‑1 (MQ=255)
                     aCAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTTCCCGCGCATGCGTCATT‑GAATTTGCCGGCGCCTTAGCGGGTaaaac                                                                                                                               >  2:1663163/1‑148 (MQ=255)
                      cAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAAccc                                                                                                                              <  1:4063022/150‑1 (MQ=255)
                       aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                                <  1:4728946/131‑1 (MQ=255)
                       aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                                >  2:1989606/1‑131 (MQ=255)
                       aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                                <  1:1989606/131‑1 (MQ=255)
                       aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                                >  2:4728946/1‑131 (MQ=255)
                         cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCt                                                                                                                           <  2:4166837/150‑1 (MQ=255)
                          gccgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACCGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCCACGCATTCGCGGGTAAACCCGCTc                                                                                                                          <  2:4597795/148‑1 (MQ=255)
                                    aGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAggg                                                                                                                <  2:1748011/150‑1 (MQ=255)
                                         aaTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGAATTCGCGGGTAAACCCGCTCCTACACAGGGGCGtt                                                                                                           <  2:1642017/150‑1 (MQ=255)
                                         aaTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTGCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGtt                                                                                                           <  2:3268618/150‑1 (MQ=255)
                                          aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTg                                                                                                          <  2:1770353/150‑1 (MQ=255)
                                          aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTg                                                                                                          <  2:1770174/150‑1 (MQ=255)
                                           tGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTCCCGCTGACGGGTGGGCGC‑GGT‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGc                                                                                                         <  1:4932882/150‑1 (MQ=255)
                                            gCCCGTTCGACACCGGGCGAAGTTTTGGGGGTGTATGCGGGGGTTC‑CCTACAGCTTCCGGTGACGGGTGGGCGC‑GGT‑TTTACCCGCGAATGCGTATTT‑GAATTTGCCGACGAATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGca                                                                                                        <  2:1430412/150‑1 (MQ=255)
                                            gCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGca                                                                                                        <  2:4329158/150‑1 (MQ=255)
                                              ccGTTCGACACCGGGCGATGGTTTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGcaca                                                                                                      <  1:4953213/150‑1 (MQ=255)
                                               cGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GGTCCCCGCGCATGCGGCCTG‑GCATTTGCCCCCGCCGGCCCGGGGAAACCCGCGCCGACACCGGGGCGGTGCAcag                                                                                                     >  2:3573696/1‑150 (MQ=255)
                                                gTTCGACACCGGGCGATGGATTGGGCGTGTATGCCGGGGTTT‑CCTCCAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagc                                                                                                    <  2:3735245/150‑1 (MQ=255)
                                                 ttCGACACCGGGCTAAGTATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTCCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCTAATGCTTCATT‑TAATTGGCCCACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                                   <  2:622517/150‑1 (MQ=38)
                                                      cacCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTTACCCCCGAAT‑CGTCATT‑GAATTTTCCCACGCAGTCGCGGGGAAACCCCCGCCTACACAGGGGGGGTGGACAGCAGGGtt                                                                                              >  1:3597797/1‑150 (MQ=255)
                                                       acCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTAACCCCGAATTCGTCATT‑GTATTTTGCGACCCATTCGCCGGTAAAACCCCTCCAAAAAAGGGGCGGTTCCCAGCAGGCTTa                                                                                             >  2:3405312/1‑150 (MQ=255)
                                                       acCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTACCCGCGACAGCGGCCAT‑GCAGTTTGCGCCGCAGTCGCGGGGGAACCCGCGCCGACACCGGGGCGGTGCCCAGCCggggtt                                                                                             >  2:331027/1‑146 (MQ=255)
                                                        ccGGGCGATGGATTTGGCGTTTATGCCGGGTTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGG‑GGT‑TTTACCCGCGAATGCTTCATT‑GAATTTGCCGACGCTTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAg                                                                                            <  1:801518/150‑1 (MQ=255)
                                                          gggCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGCGGGG‑TTAACCCCCCAAT‑CGTCATT‑GAAATTTCCGACCCAATAGCGGGTAAAACCCCAACAACACACGGGCGGTGCACAGCCGGGATAggg                                                                                          >  2:4205153/1‑150 (MQ=255)
                                                          gggCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTAACCCCCAAAGCGGGAAT‑GTATTTTGCCACCCAATAGCGGGTAAAACCCCACCAAAACACGGGCGGTGCACAGCACGATTAggg                                                                                          >  2:4205146/1‑150 (MQ=255)
                                                              gacggATTGGGCGTGTCTGCCGGGGTGT‑CCTCCCGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTc                                                                                      <  2:583681/147‑1 (MQ=38)
                                                               ttgtttttGGCGTTTTTGCGGGGTTTT‑CCTACTGCTTCCGGTGCGGGGGGGGGGG‑GGT‑TTTACCCGCTAATGCTTCTTT‑TATTTTGCCGACGAATGGGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                                     <  1:1119245/145‑1 (MQ=255)
                                                                    tttGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTTACCGCGAATTCGTCATT‑GAATTTTCCGACGCATTCGCGGGTAAACCCGCGCCTAAAAAGGGGCGTTGCACAGCAGGGTTAGGGGTTGAGACCg                                                                                >  2:329190/1‑150 (MQ=255)
                                                                      tGGCGCGTAGGCCGGGGTTT‑CCAACCGCGCCCGCTGACGGGTGGGCGC‑GGG‑TGGCCCCGCGAATGCGACATG‑GAATTTGCCGACGCATGCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                              <  2:5097005/150‑1 (MQ=255)
                                                                       ggCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTTCCCCCGAAAGCGGCCTTGGAATTTGCC‑ACGCAATAGCGGGTAAAACCCCGCCGACACAAGGGCGGTGCACCCCCGGATTAGGGGTTGAGAcccccc                                                                             >  2:1187038/1‑150 (MQ=255)
                                                                        gCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTAACCCCGAAAGCGTCATT‑GAATTTTGCGACCCATTAGCGGGTAAAACCGCACCTAAAAAAGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcccccca                                                                            >  2:3502255/1‑149 (MQ=255)
                                                                        gCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTAACCCCGAATTGGGGATT‑TAATTTTGCGACGCAATGGCGGGTTAACCCCCACCAAAACACGGGCGGTTCACAGCAGGATTATGGGTTAAGAACccccc                                                                            >  2:3502333/1‑149 (MQ=255)
                                                                          gtgtATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTACCCGCGAATGGGGCCTT‑TGGTTTGGCGGCGCATTCGCGGGGTAACCCCCTCCGACACAGGGGCGGTGCACCGCCGGGTTAGGGGTTGAGACCCCCCGGt                                                                          >  2:4798119/1‑150 (MQ=255)
                                                                            gtATGCGGGGGTTT‑CCTACCGCTTCCGCTGTGGGGTGGGCGC‑GGG‑TTTACCCGCGATTGCTTATTT‑TATTTTGCCGACGCTTTGGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                        <  1:4205153/150‑1 (MQ=255)
                                                                            gtATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                        <  2:435105/150‑1 (MQ=255)
                                                                            gtATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                        <  2:4972664/150‑1 (MQ=255)
                                                                             tttggcGGGGTTTGCCTACCGCTGCGGCTGCGGGGGGGGCGC‑GGG‑TTTCCCCGGTATTGCTTATTT‑GATTTGGCCGACGATTTCGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                        <  1:4205146/145‑1 (MQ=255)
                                                                              ttGCGGGGTTTT‑CCTCCTGCTTCCGCTTCGGGGTGGGGGG‑GGT‑TTTCCCGGCTAGTGCTTCATT‑TATTTTGCCGCCGCTTTGGGGGTAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                                      <  1:4674234/149‑1 (MQ=255)
                                                                              tggcgggggTTT‑CCTACCGCTTCCGCTGAGGGGTGGGCGG‑GGT‑TTTCCCCGCTAATGCGTATTT‑TAATTGGCCGACGCATTCGGGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                                      <  1:3237869/145‑1 (MQ=255)
                                                                              aTGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                                      <  1:4955390/150‑1 (MQ=255)
                                                                               tGCCGGGGTTT‑CCTACAGCTGCCGCTGACGGGTGGGCGC‑GGG‑TTTACCGGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAAc                                                                     <  2:4656945/150‑1 (MQ=255)
                                                                                     ggTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑TAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                               <  2:4707938/150‑1 (MQ=255)
                                                                                     ggTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                               <  2:4354594/150‑1 (MQ=255)
                                                                                        tt‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCAt                                                            <  2:1590183/150‑1 (MQ=255)
                                                                                           ccTACTGCTTCCGCTGACGGGTGGGCGG‑GGG‑TTTACCCGCTATTGCGTAATT‑TAATTTGCCGACGCTTTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGa                                                          <  1:3502255/150‑1 (MQ=255)
                                                                                           ccTACAGCTCCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCTAATGCGTCATT‑TATTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGa                                                          <  1:3502333/150‑1 (MQ=255)
                                                                                            atACTGCTTCGGGTGCGGGGTGGGGGG‑GGT‑TTTACCCGCTATTGCTTCTTT‑TATTTTGCCGCCGATTTCGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                         <  1:189902/149‑1 (MQ=255)
                                                                                                  cTTCCGCTGAGGGGTGGGCGC‑GGG‑TTTACCCGCGATTGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                                   <  2:2017390/150‑1 (MQ=255)
                                                                                                    tCCGCTGCGGGGGGGGCGC‑GGG‑TTGACCCGCGAAGGCGTCATT‑GAATTTGCCGACGCTTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                                 <  1:1187038/150‑1 (MQ=255)
                                                                                                      cggggccGGGTGGGCGC‑GGG‑TTTACCCGCGATTGCGTCTTT‑GAATTTGCCGCCGCTTTGGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                                               <  2:2274546/144‑1 (MQ=255)
                                                                                                         ggACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGAATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGa                                            <  2:1587192/149‑1 (MQ=255)
                                                                                                                  ggCGC‑GGG‑TTTACCCGCGAATGCGTAATT‑TATTTTGCCGACGCTTTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                                   <  1:3281242/150‑1 (MQ=255)
                                                                                                                   gcgc‑GGG‑TTTACCCGCTAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                                  <  2:436750/150‑1 (MQ=255)
                                                                                                                   gcgc‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                                  <  2:437016/150‑1 (MQ=255)
                                                                                                                   gcgc‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                                  <  2:436768/150‑1 (MQ=255)
                                                                                                                     ggggggttttACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                                  <  2:2341525/143‑1 (MQ=255)
                                                                                                                     gc‑gGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCa                                >  1:2443995/1‑150 (MQ=255)
                                                                                                                          g‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                                          <  2:1633197/104‑1 (MQ=255)
                                                                                                                          g‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                                          >  1:1633197/1‑104 (MQ=255)
                                                                                                                          g‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                            <  2:1689060/150‑1 (MQ=255)
                                                                                                                                cccGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                       <  2:2443995/150‑1 (MQ=255)
                                                                                                                                 ccGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  1:2396191/1‑150 (MQ=37)
                                                                                                                                 ccGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  2:237113/1‑150 (MQ=37)
                                                                                                                                 ccGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATCACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  2:5112792/1‑150 (MQ=37)
                                                                                                                                 ccGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  1:5369882/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  2:2874174/1‑150 (MQ=255)
                                                                                                                                  cgcgAATGCGTCTTT‑GAATTGGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     <  1:4296739/150‑1 (MQ=255)
                                                                                                                                  cgcgAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     >  1:5266800/1‑150 (MQ=255)
                                                                                                                                  cgcgAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     <  2:821310/150‑1 (MQ=255)
                                                                                                                                   gcgACTGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACgggcatt                    <  2:2145972/150‑6 (MQ=255)
                                                                                                                                     gAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTc                  >  2:1947230/1‑150 (MQ=255)
                                                                                                                                                  gAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGc      >  1:1423268/1‑150 (MQ=255)
                                                                                                                                                     tttGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTg   >  1:343885/1‑150 (MQ=255)
                                                                                                                                                     tttGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTg   >  1:1031708/1‑150 (MQ=255)
                                                                                                                                                      ttGCCGACGCATTCGCGGGTACCCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTATGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGc  <  2:1829990/150‑1 (MQ=255)
                                                                                                                                                      ttGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                                          <  1:843712/78‑1 (MQ=255)
                                                                                                                                                      ttGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                                          >  2:843712/1‑78 (MQ=255)
                                                                                                                                                         |                                                                                                                                                  
CAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGC  >  NC_002947/4677163‑4677458

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