Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,807,784 C→T 6.8% R800R (CGC→CGT)  PP_3364 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,807,7840CT 6.8% 173.3 / 2.2 59R800R (CGC→CGT) PP_3364two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base T (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

TCAACGATGCGCAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTGCTTGAGCTGGGGT  >  NC_002947/3807636‑3807926
                                                                                                                                                    |                                                                                                                                                 
tcAACGATGCGCAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACGAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGc                                                                                                                                                   >  5:117133/1‑149 (MQ=255)
  aaCGATGCGCAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCg                                                                                                                                                 <  8:289032/149‑1 (MQ=255)
    cGATGCGCAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCgg                                                                                                                                                >  4:278459/1‑148 (MQ=255)
     gATGCGCAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCgg                                                                                                                                                >  7:73912/1‑147 (MQ=255)
       tGCGCAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAgg                                                                                                                                            >  2:274874/1‑149 (MQ=255)
          gcAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGgttc                                                                                                                                         >  3:317344/1‑149 (MQ=255)
           cAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACATGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGcc                                                                                                                                                  <  1:263510/139‑1 (MQ=255)
           cAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCAGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGcc                                                                                                                                                  >  2:263510/1‑139 (MQ=255)
            aaaTCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGgttccg                                                                                                                                       >  6:66373/1‑149 (MQ=255)
               tCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGgttccgtt‑c                                                                                                                                   <  1:226165/149‑1 (MQ=255)
                 aaGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCTATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGgttccgtt‑cc                                                                                                                                  <  6:139219/148‑1 (MQ=255)
                       gTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGc                                                                                                                           >  7:264060/1‑149 (MQ=255)
                                gCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTa                                                                                                                  >  7:22361/1‑149 (MQ=255)
                                  cAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATg                                                                                                                <  2:58794/149‑1 (MQ=255)
                                     cTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCtt                                                                                                                              >  8:323947/1‑132 (MQ=255)
                                     cTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCtt                                                                                                                              <  7:323947/132‑1 (MQ=255)
                                                 aGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCa                                                                                               <  5:66373/149‑1 (MQ=255)
                                                     aTGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTcc                                                                                           <  1:274874/149‑1 (MQ=255)
                                                     aTGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTcc                                                                                           >  3:193647/1‑149 (MQ=255)
                                                      tGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTccc                                                                                          >  2:40410/1‑149 (MQ=255)
                                                       gCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTccc                                                                                          >  7:348656/1‑148 (MQ=255)
                                                        cAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCACCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGa                                                                                        <  7:106131/149‑1 (MQ=255)
                                                         aaCTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAg                                                                                       <  8:22361/149‑1 (MQ=255)
                                                              cTTGGCAGCCACGAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCAc                                                                                    <  6:117133/147‑1 (MQ=255)
                                                              cTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACtt                                                                                  >  7:114987/1‑149 (MQ=255)
                                                               ttGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTg                                                                                 >  6:40407/1‑149 (MQ=255)
                                                                 ggCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTgc                                                                                <  8:264060/148‑1 (MQ=255)
                                                                   cAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTgcagc                                                                             <  4:193647/149‑1 (MQ=255)
                                                                      ccACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCAtg                                                                          >  8:11166/1‑149 (MQ=255)
                                                                                  cgcgGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTcgcgc                                                              <  1:29920/149‑1 (MQ=255)
                                                                                     ggCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCAtt                                                                                             <  3:304621/115‑1 (MQ=255)
                                                                                     ggCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCAtt                                                                                             >  4:304621/1‑115 (MQ=255)
                                                                                           aCTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTCGTTCGCCGGCAGGTGCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGCGCACTTGCAGCATGTGACCGTCGCGCGCAAGccg                                                      >  2:304777/1‑147 (MQ=255)
                                                                                             tGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTa                                                    >  8:322345/1‑148 (MQ=255)
                                                                                             tGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTAt                                                   >  8:197704/1‑149 (MQ=255)
                                                                                             tGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTAt                                                   >  3:191772/1‑149 (MQ=255)
                                                                                               gATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTa                                                    <  7:322345/146‑1 (MQ=255)
                                                                                               gATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTAt                                                   <  5:40407/147‑1 (MQ=255)
                                                                                               gATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATg                                                  <  8:348656/148‑1 (MQ=255)
                                                                                                aTGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGAc                                                <  8:73912/149‑1 (MQ=255)
                                                                                                aTGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGAc                                                <  8:114987/149‑1 (MQ=255)
                                                                                                                cGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGTCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTccc                                                                                          <  4:6597/91‑1 (MQ=255)
                                                                                                                cGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGTCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTccc                                                                                          <  4:11979/91‑1 (MQ=255)
                                                                                                                cGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGTCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTccc                                                                                          >  3:6597/1‑91 (MQ=255)
                                                                                                                cGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGTCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTccc                                                                                          >  3:11979/1‑91 (MQ=255)
                                                                                                                cGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCt                                                                                                                   <  4:248077/68‑1 (MQ=255)
                                                                                                                cGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCt                                                                                                                   >  3:248077/1‑68 (MQ=255)
                                                                                                                  gAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGGGTAGAAGCATAAGCATTACCGAGCACTTGAAGCATGTGACCCTCTCGCGCAAGCCTTATGACGCCAATGACATTATgt                                >  1:115332/1‑148 (MQ=255)
                                                                                                                     ccTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGc                           >  1:114739/1‑149 (MQ=255)
                                                                                                                             tCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCaaa                    <  7:11166/148‑1 (MQ=255)
                                                                                                                             tCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCaaa                    <  2:114739/148‑1 (MQ=255)
                                                                                                                             tCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCaaa                    <  1:40410/148‑1 (MQ=255)
                                                                                                                                 tAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTTGCCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCt                >  1:49337/1‑149 (MQ=255)
                                                                                                                                 tAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTg               >  8:9823/1‑149 (MQ=255)
                                                                                                                                    cGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTGCtt            <  3:278459/149‑1 (MQ=255)
                                                                                                                                      gATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGcc                 <  2:90869/142‑1 (MQ=255)
                                                                                                                                      gATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGcc                 >  1:90869/1‑142 (MQ=255)
                                                                                                                                           gCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTGCTTGAGCTgg     >  6:363000/1‑149 (MQ=255)
                                                                                                                                            cTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑TCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTGCTTGAGCTggg    <  5:363000/149‑1 (MQ=255)
                                                                                                                                               gTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTGCTTGAGCTgggg   <  7:197704/147‑1 (MQ=255)
                                                                                                                                               gTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTGCTTGAGCTgggg   <  2:49337/147‑1 (MQ=255)
                                                                                                                                               gTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTGCTTGAGCTGGGGt  >  5:70308/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCAACGATGCGCAAATCAAGGACGTGATCACCGCCAGCTCCGAGGACCAAGCCATGCAACTGCTTGGCAGCCACAAGCCGGACGCGGCAATACTGGATGTCAACCTGGGCACCGGAACCTCCATTTCCGTAGCGGATGAGCTGGTTCGCCGGCAGGTTCCGTT‑CCTCTTTGCCACCGGCTATGGCGACG‑‑GCATCAGCATTCCCGAGCACTTGCAGCATGTGCCCGTCGCGCGCAAGCCTTATGACGCCAATGCCATTCTGGCCAGCCTGCAAAGCCTGCTTGAGCTGGGGT  >  NC_002947/3807636‑3807926

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