Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,494,824 C→T 8.1% V388M (GTG→ATG)  cobJ ← precorrin‑3B C17‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,494,8240CT 8.1% 475.4 / 22.4 161V388M (GTG→ATG) cobJprecorrin‑3B C17‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (68/80);  new base T (5/8);  total (73/88)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

ACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGC  >  NC_002947/5494675‑5494971
                                                                                                                                                     |                                                                                                                                                   
aCCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGAcaccac                                                                                                                                                     >  2:610815/1‑150 (MQ=255)
     cctgGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACAccaccgcca                                                                                                                                                >  1:5401651/1‑150 (MQ=255)
      ctgGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACAccaccgccac                                                                                                                                               <  2:2151371/150‑1 (MQ=255)
       tgGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACAccaccgccacc                                                                                                                                              >  1:3399111/1‑150 (MQ=255)
        gggCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACAccaccgccaccg                                                                                                                                             >  1:2030795/1‑150 (MQ=255)
          ggggcggTTGCCAGCGAGGCGTATACCCCCGGCAAGATCTGCAGGTCAACGCGTTGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGaa                                                                                                                                           <  1:4624324/148‑1 (MQ=255)
          gcggcgGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGaa                                                                                                                                           <  2:5401651/150‑1 (MQ=255)
          gcggcgGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGaa                                                                                                                                           <  2:3727493/150‑1 (MQ=255)
           cggcggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAAc                                                                                                                                          >  1:2145854/1‑150 (MQ=255)
           cggcggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCAGGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAAc                                                                                                                                          <  2:639752/150‑1 (MQ=255)
              cggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGc                                                                                                                                       <  2:3338689/150‑1 (MQ=255)
                   gCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTgg                                                                                                                                  <  2:828044/150‑1 (MQ=255)
                   gCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTgg                                                                                                                                  <  2:1535439/150‑1 (MQ=255)
                   gCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTgg                                                                                                                                  <  1:2502383/150‑1 (MQ=255)
                     cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGAGTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                                >  1:1775767/1‑150 (MQ=255)
                      aGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgcc                                                                                                                               <  2:1483688/150‑1 (MQ=255)
                        cGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgccgc                                                                                                                             <  1:3039002/150‑1 (MQ=255)
                          aGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCa                                                                                                                           >  1:4647302/1‑150 (MQ=255)
                              ggATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTc                                                                                                                       >  1:261437/1‑150 (MQ=255)
                                  aCCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAg                                                                                                                   <  2:4556908/150‑1 (MQ=255)
                                    ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                                 >  2:275326/1‑150 (MQ=255)
                                    ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                                 >  1:3436424/1‑150 (MQ=255)
                                    ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                                 >  1:4386758/1‑150 (MQ=255)
                                     cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGaa                                                                                                                    <  1:5326261/146‑1 (MQ=255)
                                     cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGaa                                                                                                                    >  2:5326261/1‑146 (MQ=255)
                                     cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCa                                                                                                                >  2:1826349/1‑150 (MQ=255)
                                     cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCa                                                                                                                >  1:1704259/1‑150 (MQ=255)
                                                 tGCGGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                    >  2:2689907/1‑150 (MQ=255)
                                                   cAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTg                                                                                                  >  1:1183014/1‑150 (MQ=255)
                                                     ggTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCa                                                                                                <  1:665237/150‑1 (MQ=255)
                                                      gTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAt                                                                                               >  1:4501816/1‑150 (MQ=255)
                                                          aCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATctct                                                                                           <  2:3753171/150‑1 (MQ=255)
                                                                tGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCgg                                                                                     <  2:2145855/150‑1 (MQ=255)
                                                                      gTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                               <  2:577064/150‑1 (MQ=255)
                                                                      gTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                               <  2:4890140/150‑1 (MQ=255)
                                                                       tGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCgg                                                                              <  2:5031421/150‑1 (MQ=255)
                                                                        ggggTCGCTGGAGTCGTGCACTGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCGGCATCTCTTCGCGGTTGTCGgt                                                                             <  2:5260100/150‑1 (MQ=255)
                                                                           gTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGTCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgca                                                                          <  2:3940599/150‑1 (MQ=255)
                                                                           gTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGTCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgca                                                                          <  2:1181975/150‑1 (MQ=255)
                                                                           gTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgca                                                                          >  2:2301258/1‑150 (MQ=255)
                                                                             cGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagt                                                                        >  2:2497782/1‑150 (MQ=255)
                                                                             cGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagt                                                                        <  2:4030809/150‑1 (MQ=255)
                                                                              gcgggATTCGTGCACTGCTTCGAGCACCGCAGCGGCCAGGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtg                                                                       <  2:4881227/147‑1 (MQ=255)
                                                                              gCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtg                                                                       >  1:2572546/1‑150 (MQ=255)
                                                                               cTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgc                                                                      >  1:5447617/1‑150 (MQ=255)
                                                                                tGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                     <  2:4158177/150‑1 (MQ=255)
                                                                                tGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                     <  1:232764/150‑1 (MQ=255)
                                                                                  gATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTgcagc                                                                   <  2:2440567/150‑1 (MQ=255)
                                                                                  gATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTgcagc                                                                   >  2:1878873/1‑150 (MQ=255)
                                                                                       gTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCAccga                                                              <  2:27398/150‑3 (MQ=255)
                                                                                        tGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTgg                                                             >  1:4926199/1‑150 (MQ=255)
                                                                                           aaTGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcg                                                          >  1:3956622/1‑150 (MQ=255)
                                                                                              gCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccg                                                       <  1:5137799/150‑1 (MQ=255)
                                                                                               cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                      >  2:4702637/1‑150 (MQ=255)
                                                                                               cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                      >  2:2668780/1‑150 (MQ=255)
                                                                                               cTTCGAGCACCGCAGAGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                               >  1:3046958/1‑125 (MQ=255)
                                                                                                ttCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCcgcg                                                     >  2:4909461/1‑150 (MQ=255)
                                                                                                   gAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGt                                                                                 <  1:4725159/119‑1 (MQ=255)
                                                                                                     gCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCACCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACg                                                >  2:3484468/1‑150 (MQ=255)
                                                                                                      cACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACgg                                               <  2:3399111/150‑1 (MQ=255)
                                                                                                      cACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACgg                                               <  2:2394070/150‑1 (MQ=255)
                                                                                                        ccGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgc                                             >  2:1821500/1‑150 (MQ=255)
                                                                                                          gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcca                                           >  2:3181006/1‑150 (MQ=255)
                                                                                                          gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcca                                           <  1:3699410/150‑1 (MQ=255)
                                                                                                           cAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgccag                                          >  1:5055794/1‑150 (MQ=255)
                                                                                                           cAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgccag                                          <  1:1874248/150‑1 (MQ=255)
                                                                                                                    tGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCa                                 >  1:3018642/1‑150 (MQ=255)
                                                                                                                     ggCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCAc                                >  2:2072823/1‑150 (MQ=255)
                                                                                                                     ggCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCAc                                <  1:5302910/150‑1 (MQ=255)
                                                                                                                        gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                    >  1:4442048/1‑79 (MQ=255)
                                                                                                                        gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTa                             >  1:333463/1‑150 (MQ=255)
                                                                                                                         aaaaCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTaa                            <  2:5203002/150‑1 (MQ=255)
                                                                                                                         aaaaCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTaa                            <  2:4395782/150‑1 (MQ=255)
                                                                                                                            aCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGtt                         >  1:1095312/1‑150 (MQ=255)
                                                                                                                            aCGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGtt                         <  1:1928505/150‑1 (MQ=255)
                                                                                                                             cGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGttt                        <  2:1095312/150‑1 (MQ=255)
                                                                                                                             cGCCCGGATCGCCTGACGATACCATCACCACATCCCGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGttt                        <  2:3139159/150‑1 (MQ=255)
                                                                                                                             cGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGttt                        <  1:4077252/150‑1 (MQ=255)
                                                                                                                              gCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                       >  2:1236431/1‑150 (MQ=255)
                                                                                                                               cccGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTca                      >  1:477503/1‑149 (MQ=255)
                                                                                                                                ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTgg                                                             >  1:525683/1‑110 (MQ=255)
                                                                                                                                ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTgg                                                             <  2:525683/110‑1 (MQ=255)
                                                                                                                                ccGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTcat                     >  2:2997130/1‑148 (MQ=255)
                                                                                                                                ccGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTcat                     >  1:4104977/1‑148 (MQ=255)
                                                                                                                                 cGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATa                    <  1:3595887/150‑1 (MQ=255)
                                                                                                                                       gCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAg              >  1:4144341/1‑150 (MQ=255)
                                                                                                                                        ccTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGc             <  2:2189310/150‑1 (MQ=255)
                                                                                                                                        ccTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGc             <  2:2374939/150‑1 (MQ=255)
                                                                                                                                        ccTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGc             <  2:3113598/150‑1 (MQ=255)
                                                                                                                                         cTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGca            <  2:83706/150‑1 (MQ=255)
                                                                                                                                         cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca            >  2:4636378/1‑150 (MQ=255)
                                                                                                                                         cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca            >  1:2093177/1‑150 (MQ=255)
                                                                                                                                         cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca            <  2:1704259/150‑1 (MQ=255)
                                                                                                                                         cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca            >  1:1270046/1‑150 (MQ=255)
                                                                                                                                           gacgacACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGcaca          <  1:275326/150‑1 (MQ=255)
                                                                                                                                                caccacCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                       <  2:364642/132‑1 (MQ=255)
                                                                                                                                                 accaccGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCg    >  2:4736684/1‑150 (MQ=255)
                                                                                                                                                 accaccGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCg    >  1:2147756/1‑150 (MQ=255)
                                                                                                                                                   caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGc  >  1:5159363/1‑150 (MQ=255)
                                                                                                                                                   caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTACGTAACCGAGCACATCTTCGGc  >  1:2945029/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGC  >  NC_002947/5494675‑5494971

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