Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,677,313 C→T 9.0% intergenic (+44/+60) 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,3130CT 9.0% 275.7 / 13.4 100intergenic (+44/+60)PP_4137/chrRouter membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand):  ref base C (27/63);  new base T (4/5);  total (32/68)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGC  >  NC_002947/4677167‑4677458
                                                                                                                                                      |                                                                                                                                                 
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:3763219/150‑1 (MQ=255)
        cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTcgc                                                                                                                                        <  1:4707577/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:443606/150‑1 (MQ=37)
          ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                                      <  2:3402003/150‑1 (MQ=37)
          ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                                      <  2:2634065/150‑1 (MQ=37)
          ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                                      <  2:5300967/150‑1 (MQ=255)
          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‑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:4728946/1‑131 (MQ=255)
                   aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                                >  2:1989606/1‑131 (MQ=255)
                   aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                                <  1:1989606/131‑1 (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:1770174/150‑1 (MQ=255)
                                      aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTg                                                                                                          <  2:1770353/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:4972664/150‑1 (MQ=255)
                                                                        gtATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                        <  2:435105/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                                                                          >  1:1633197/1‑104 (MQ=255)
                                                                                                                      g‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                                          <  2:1633197/104‑1 (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)
                                                                                                                             ccGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCGGGTAAACCCGCGCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  2:988838/1‑150 (MQ=37)
                                                                                                                             ccGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  2:2874174/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  1:5369882/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)
                                                                                                                                                      |                                                                                                                                                 
GACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGC  >  NC_002947/4677167‑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: ‑