Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,494,819 C→T 6.1% V389V (GTG→GTA cobJ ← precorrin‑3B C17‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,494,8190CT 6.1% 984.4 / 31.9 312V389V (GTG→GTAcobJprecorrin‑3B C17‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (152/141);  new base T (13/6);  total (165/147)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 319 total aligned reads displayed.

GGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTC  >  NC_002947/5494670‑5494968
                                                                                                                                                     |                                                                                                                                                     
 gTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGaca                                                                                                                                                      >  1:4478685/1‑150 (MQ=255)
    cACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGAcacca                                                                                                                                                   <  2:425148/150‑1 (MQ=255)
      cctgcctgGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACaccacc                                                                                                                                                 <  2:3452964/150‑1 (MQ=255)
           ctgGGCGGCGGTTGCCATCGAGGCGGATTCCCCCGGCATGTTCTGCATGTCAACGCGATGCCTGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACAccaccgccac                                                                                                                                            <  2:4173427/150‑1 (MQ=255)
               gcggcgGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGaa                                                                                                                                        >  1:2968135/1‑150 (MQ=255)
                     gTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCt                                                                                                                                  >  1:846669/1‑150 (MQ=255)
                       tGCCAGCGAGGCGTATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACCCCACCGCCACCGAACGGCCTTg                                                                                                                                <  2:279519/150‑1 (MQ=255)
                        gCCAGCGAGGCGGATACCCCCGGCATGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTACGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTgg                                                                                                                               <  2:8522532/150‑1 (MQ=255)
                        gCCAGCGAGGCGGATACCCCCGGCAAGCTCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTgg                                                                                                                               <  1:630232/150‑1 (MQ=255)
                          cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCCGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCAGGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                             <  2:441650/150‑1 (MQ=255)
                          cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                             >  2:184437/1‑150 (MQ=255)
                          cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                             >  2:8941162/1‑150 (MQ=255)
                          cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                             >  1:1002056/1‑150 (MQ=255)
                               tggCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCTGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCa                                                                                                                        <  1:3303176/149‑1 (MQ=255)
                               aGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCa                                                                                                                        >  2:4944166/1‑150 (MQ=255)
                                 gCGGAGACCCCCGGCAAGATCTGCCGGTCACCGCGATGCCCGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGt                                                                                                                      <  1:6462881/150‑1 (MQ=255)
                                    gATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCg                                                                                                                   >  2:728668/1‑150 (MQ=255)
                                    gATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCg                                                                                                                   >  1:6691003/1‑150 (MQ=255)
                                    gATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGCCCTTGGGCCGCCAGTTCg                                                                                                                   >  1:2701167/1‑150 (MQ=255)
                                      tACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGTGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGaa                                                                                                                 <  1:1183712/150‑1 (MQ=255)
                                      tACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGaa                                                                                                                 <  2:2328285/150‑1 (MQ=255)
                                      tACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGaa                                                                                                                 <  1:4764329/150‑1 (MQ=255)
                                       aCCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCCGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAg                                                                                                                <  2:7653521/150‑1 (MQ=255)
                                         ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAATGc                                                                                                              >  1:2307033/1‑150 (MQ=255)
                                         ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                              <  2:1191062/150‑1 (MQ=255)
                                         ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                              <  1:8352400/150‑1 (MQ=255)
                                         ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                              >  2:6447818/1‑150 (MQ=255)
                                         ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                              >  2:398938/1‑150 (MQ=255)
                                         ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                              >  1:616636/1‑150 (MQ=255)
                                           ccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCAt                                                                                                            <  1:3547367/150‑1 (MQ=255)
                                           ccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCAt                                                                                                            >  1:2012456/1‑150 (MQ=255)
                                             ggCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGACCGGCCTTGGGCCGCCAGTTCGAAGGCATGt                                                                                                          <  2:258064/150‑1 (MQ=255)
                                                aaGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCgg                                                                                                       <  2:5562260/150‑1 (MQ=255)
                                                aaGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCgg                                                                                                       <  2:459803/150‑1 (MQ=255)
                                                 aGATCTGCAGGTCAACGCGAGGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCggg                                                                                                      <  1:8942383/150‑1 (MQ=255)
                                                  gATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgc                                                                                                     >  1:5603507/1‑150 (MQ=255)
                                                   aTCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcg                                                                                                    >  1:4919263/1‑150 (MQ=255)
                                                     cTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcg                                                                                                  >  2:5050082/1‑150 (MQ=255)
                                                      tGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                 >  1:7060500/1‑150 (MQ=255)
                                                      tGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                 <  1:4010885/150‑1 (MQ=255)
                                                          ggTCAACGCGATGCCTGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCa                                                                                             <  1:6602351/150‑1 (MQ=255)
                                                           gTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAt                                                                                            >  1:681808/1‑150 (MQ=255)
                                                             cAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAtct                                                                                          >  2:5005513/1‑150 (MQ=255)
                                                              aaCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAtctc                                                                                         >  2:6859800/1‑150 (MQ=255)
                                                              aaCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAtctc                                                                                         >  2:3341201/1‑150 (MQ=255)
                                                              aaCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAtctc                                                                                         >  1:3176522/1‑150 (MQ=255)
                                                              aaCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAtctc                                                                                         >  1:5907439/1‑150 (MQ=255)
                                                                cgcgATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCtt                                                                                       >  1:7850228/1‑150 (MQ=255)
                                                                       ccAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGtt                                                                                >  1:2836958/1‑150 (MQ=255)
                                                                       ccAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGtt                                                                                <  1:8190040/150‑1 (MQ=255)
                                                                           gTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                            <  2:3944882/150‑1 (MQ=255)
                                                                           gTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                            <  1:7390447/150‑1 (MQ=255)
                                                                            tGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCgg                                                                           <  2:1462966/150‑1 (MQ=255)
                                                                             ggggTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCTCCGCCACCGAACGGCCTTGGGCCGCCAGTTCGCAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgt                                                                          <  2:6800209/150‑1 (MQ=255)
                                                                             ggggTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgt                                                                          <  1:9034366/150‑1 (MQ=255)
                                                                              gggTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtg                                                                         <  1:4359551/150‑1 (MQ=255)
                                                                              gggTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtg                                                                         <  2:7507277/150‑1 (MQ=255)
                                                                                gTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgca                                                                       <  2:2088688/150‑1 (MQ=255)
                                                                                 tCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcag                                                                      >  2:3045614/1‑150 (MQ=255)
                                                                                     tGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                  <  2:7187104/150‑1 (MQ=255)
                                                                                     gggATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                  <  2:4894867/149‑1 (MQ=255)
                                                                                         ttCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgcacgcc                                                              >  1:3693880/1‑146 (MQ=255)
                                                                                         ttCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCAc                                                              >  1:7666840/1‑150 (MQ=255)
                                                                                          tCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCAcc                                                             >  2:5872385/1‑150 (MQ=255)
                                                                                            gTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTg                                                           <  2:5741921/150‑1 (MQ=255)
                                                                                             tGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTgg                                                          >  2:3136015/1‑150 (MQ=255)
                                                                                             tGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTgg                                                          >  2:8307326/1‑150 (MQ=255)
                                                                                                aaTGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcg                                                       >  1:1488457/1‑150 (MQ=255)
                                                                                                  tGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgcc                                                     <  2:2521949/150‑1 (MQ=255)
                                                                                                   gCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccg                                                    >  1:7240566/1‑150 (MQ=255)
                                                                                                    cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                   <  2:7796300/150‑1 (MQ=255)
                                                                                                     ttCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCcgcg                                                  >  2:3118001/1‑150 (MQ=255)
                                                                                                     ttCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCcgcg                                                  <  1:7449855/150‑1 (MQ=255)
                                                                                                        gAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGaa                                               >  1:9066111/1‑150 (MQ=255)
                                                                                                        gAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGaa                                               <  2:4500972/150‑1 (MQ=255)
                                                                                                              cGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcc                                         >  2:6288129/1‑150 (MQ=255)
                                                                                                              cGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcc                                         >  2:1687306/1‑150 (MQ=255)
                                                                                                               gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGTAGCACCTGGTCGCCGCGGAACGGgcca                                        <  2:8835117/150‑1 (MQ=255)
                                                                                                                cAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgccag                                       >  2:658025/1‑150 (MQ=255)
                                                                                                                  gCGGCCTTGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgccagcc                                     <  1:3636953/150‑1 (MQ=255)
                                                                                                                      ccATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgc                                 <  2:5577289/150‑1 (MQ=255)
                                                                                                                      ccATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgc                                 >  1:7129146/1‑150 (MQ=255)
                                                                                                                       cATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgcg                                >  2:2768455/1‑150 (MQ=255)
                                                                                                                            caaataCGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCAcat                           >  1:5570857/6‑148 (MQ=255)
                                                                                                                             gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTa                          >  1:5270167/1‑150 (MQ=255)
                                                                                                                             gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTa                          >  2:5003534/1‑150 (MQ=255)
                                                                                                                               aaaCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAg                        >  1:7498353/1‑150 (MQ=255)
                                                                                                                                 aCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGtt                      <  2:743259/150‑1 (MQ=255)
                                                                                                                                   gCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                    <  2:4771949/150‑1 (MQ=255)
                                                                                                                                   gCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                    >  2:7856871/1‑150 (MQ=255)
                                                                                                                                   gCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                    >  2:3241412/1‑150 (MQ=255)
                                                                                                                                   gCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                    <  1:2625703/150‑1 (MQ=255)
                                                                                                                                     ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGt                  >  2:4128882/1‑150 (MQ=255)
                                                                                                                                            gCCTGCCGACACCACCGCCACCGAACGGCCTTGGGCCGCCCGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCCCCTGGTCGCCGCGGAACGGGCCCGCCATGCGCACGTACGTGTCGTCACCGAg           >  2:285196/1‑150 (MQ=255)
                                                                                                                                             ccTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGc          >  2:2712344/1‑150 (MQ=255)
                                                                                                                                             ccTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGc          >  1:1971685/1‑150 (MQ=255)
                                                                                                                                               tGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGcac        >  1:869000/1‑150 (MQ=255)
                                                                                                                                               tGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGcac        >  1:1538430/1‑150 (MQ=255)
                                                                                                                                               tGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGcac        >  2:759765/1‑150 (MQ=255)
                                                                                                                                                gacgacACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGcaca       <  2:6116735/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTC  >  NC_002947/5494670‑5494968

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: