Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,840 G→T 7.3% E22* (GAA→TAA) ‡ mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947635,8400GT 7.3% 193.5 / 2.1 62E22* (GAA→TAA) ‡mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base G (25/32);  new base T (3/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CATTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGC  >  NC_002947/635692‑635985
                                                                                                                                                    |                                                                                                                                                 
cATTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAgg                                                                                                                                                   <  8:248799/149‑1 (MQ=255)
    aCCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAc                                                                                                                                                <  2:366490/148‑1 (MQ=255)
      cAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTg                                                                                                                                             <  2:39659/149‑1 (MQ=255)
        ggggCGACTGAGCGAAACTGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCTTGCATATTCACATTCTCTGTATTTGGGGCACTTTGATGGGTTCGGTGGCGGTGCTGGCCAAGGAACTTGGc                                                                                                                                           <  2:196450/149‑1 (MQ=255)
        ggggCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGc                                                                                                                                           <  1:213111/149‑1 (MQ=255)
         gggCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGcc                                                                                                                                          >  4:90726/1‑149 (MQ=255)
          ggCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCa                                                                                                                                         <  5:307226/149‑1 (MQ=255)
          ggCGACTGAGCGAAACAGCGGGAGCGATTTGTCGCGCCCTTGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACGATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCa                                                                                                                                         <  5:287654/149‑1 (MQ=255)
             gACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACcg                                                                                                                                      <  6:311760/149‑1 (MQ=255)
             gACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACcg                                                                                                                                      <  6:419468/149‑1 (MQ=255)
                tGAGCGAAACAGCGGGTGGGATTTGTCGCGCCCTGGGCTACCCTGCCTCAATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGt                                                                                                                                   >  8:81442/1‑149 (MQ=255)
                   gCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCAc                                                                                                                                >  8:261799/1‑149 (MQ=255)
                   gCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCAc                                                                                                                                >  5:261830/1‑149 (MQ=255)
                   gCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGGGCTGGCCAAGGAACTTGGCCACCGCGTCAc                                                                                                                                <  1:264537/149‑1 (MQ=255)
                       aaCAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGc                                                                                                                            >  5:192005/1‑149 (MQ=255)
                         cAGCGTGGGCGATTAGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTCCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctc                                                                                                                          <  5:227944/149‑1 (MQ=255)
                         cAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctc                                                                                                                          <  8:149415/149‑1 (MQ=255)
                            cGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGa                                                                                                                       <  7:211803/149‑1 (MQ=255)
                                   gATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAAc                                                                                                                >  3:392329/1‑149 (MQ=255)
                                        gTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCt                                                                                                            <  6:231770/148‑1 (MQ=255)
                                         tCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTAt                                                                                                          >  6:83978/1‑149 (MQ=255)
                                                    gCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctct                                                                                                                         >  3:146387/1‑123 (MQ=255)
                                                    gCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctct                                                                                                                         <  4:146387/123‑1 (MQ=255)
                                                     cTACCCTGCCTCAATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGc                                                                                              <  7:81442/149‑1 (MQ=255)
                                                        cccTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTATGACGCCAACGTCTATCCCCCGATGAGCAcc                                                                                           >  7:365111/1‑149 (MQ=255)
                                                              cTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTc                                                                                     >  5:26996/1‑149 (MQ=255)
                                                               tCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTc                                                                                     >  6:131669/1‑148 (MQ=255)
                                                               tCGATGTTTTGCGAGGTCGAACATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGAGGTGCTGGCCAAGTAACTTGGCAACCGCGTCACCGGCTCTGACGCCAACGACTATCCCCCGATGAGCACCCACCTc                                                                                     >  3:27548/1‑148 (MQ=255)
                                                                cGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCg                                                                                    >  4:281150/1‑148 (MQ=255)
                                                                 gATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGaa                                                                                  <  8:282269/149‑1 (MQ=255)
                                                                       ttGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAg                                                                            >  8:323231/1‑149 (MQ=255)
                                                                         gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                           <  6:192005/148‑1 (MQ=255)
                                                                         gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                           <  7:261799/148‑1 (MQ=255)
                                                                             ggTCGACCATGCATATTCACTTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATc                                                                      >  7:9217/1‑149 (MQ=255)
                                                                                cGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGATTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGt                                                                                                              >  8:32658/1‑106 (MQ=255)
                                                                                 gACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTa                                                                        >  7:202116/1‑143 (MQ=255)
                                                                                 gACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTa                                                                        <  8:202116/143‑1 (MQ=255)
                                                                                 gACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGaa                                                                   <  5:36261/148‑1 (MQ=255)
                                                                                 gACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAAc                                                                  >  1:356587/1‑149 (MQ=255)
                                                                                    cATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGa                                                               <  3:90726/149‑1 (MQ=255)
                                                                                     aTGCATATTCACTTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGAc                                                              <  8:9217/149‑1 (MQ=255)
                                                                                        cATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGACAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGAccc                                                            <  1:279366/148‑1 (MQ=255)
                                                                                               acaTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTAt                                                    <  3:281150/149‑1 (MQ=255)
                                                                                                              gCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCt                                      <  4:392329/148‑1 (MQ=255)
                                                                                                                ggCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGa                                   <  5:83978/149‑1 (MQ=255)
                                                                                                                             gTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGAc                      >  1:215890/1‑149 (MQ=255)
                                                                                                                             gTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGAc                      >  3:124023/1‑149 (MQ=255)
                                                                                                                              ttCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGAcc                     >  4:205327/1‑149 (MQ=255)
                                                                                                                                 gCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtgg                  <  5:131669/149‑1 (MQ=255)
                                                                                                                                 gCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtgg                  >  3:186437/1‑149 (MQ=255)
                                                                                                                                     gCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtg                <  7:323231/147‑1 (MQ=255)
                                                                                                                                     gCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtg                <  8:180350/147‑1 (MQ=255)
                                                                                                                                     gCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtg                <  6:261830/147‑1 (MQ=255)
                                                                                                                                         tGCTGGCCAAGGAACTTGTCCACCGCGTCACCGGCTCTGACGCCAAAGTCTATCCCCCGATGAGCACCCAGCTCGAATCCCAGGGTATCGAACTGACCCAAGTCTATGACCCGTCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATc          <  4:27548/149‑1 (MQ=255)
                                                                                                                                          ggtGGCCAAGTTCCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGCGAACACAGCTCGAAGCCCAGGGTATCGAACTGCCCCAAGGCGATGACCCGGCCAAGCTGGATCCGGCGCCAGACCTGGTGGTCATCg         >  1:321318/3‑149 (MQ=255)
                                                                                                                                           gtGGCCAAGTTCCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCgg        >  5:157747/2‑149 (MQ=255)
                                                                                                                                           cTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCgg        <  4:124023/149‑1 (MQ=255)
                                                                                                                                              gCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcg                            >  2:387021/1‑126 (MQ=255)
                                                                                                                                              gCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcg                            <  1:387021/126‑1 (MQ=255)
                                                                                                                                                cAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACg   >  5:284733/1‑149 (MQ=255)
                                                                                                                                                 aaGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                           <  7:51415/124‑1 (MQ=255)
                                                                                                                                                 aaGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                           >  8:51415/1‑124 (MQ=255)
                                                                                                                                                 aaGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCAAAGGCTATGACCCGGCCCCGCTGGATCCGGCGCCAGACCTTGTGGTCATCGGCAACGc  >  3:310063/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CATTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGC  >  NC_002947/635692‑635985

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