Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,841 A→T 6.4% E22V (GAA→GTA) ‡ 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,8410AT 6.4% 186.9 / 2.4 62E22V (GAA→GTA) ‡mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base T (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGC  >  NC_002947/635696‑635985
                                                                                                                                                 |                                                                                                                                                
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:419468/149‑1 (MQ=255)
         gACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACcg                                                                                                                                      <  6:311760/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                                                                                                                         <  4:146387/123‑1 (MQ=255)
                                                gCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctct                                                                                                                         >  3:146387/1‑123 (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                                                                           <  7:261799/148‑1 (MQ=255)
                                                                     gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                           <  6:192005/148‑1 (MQ=255)
                                                                         ggTCGACCATGCATATTCACTTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATc                                                                      >  7:9217/1‑149 (MQ=255)
                                                                            cGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGATTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGt                                                                                                              >  8:32658/1‑106 (MQ=255)
                                                                             gACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTa                                                                        <  8:202116/143‑1 (MQ=255)
                                                                             gACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTa                                                                        >  7:202116/1‑143 (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)
                                                                                                                                                 |                                                                                                                                                
ACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGC  >  NC_002947/635696‑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: ‑