Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,531,231 T→C 7.6% intergenic (+33/+15) lpxC → / ← PP_1344 UDP‑3‑O‑[3‑hydroxymyristoyl] N‑acetylglucosamine deacetylase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,531,2310TC 7.6% 189.7 / 4.8 66intergenic (+33/+15)lpxC/PP_1344UDP‑3‑O‑[3‑hydroxymyristoyl] N‑acetylglucosamine deacetylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/27);  new base C (3/2);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

GAACAACCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGC  >  NC_002947/1531084‑1531378
                                                                                                                                                    |                                                                                                                                                    
gcacaacCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACACTCTAGTTCATTAAGTCCACCtt                                                                                                                                                     >  4:278432/3‑149 (MQ=255)
 aacaacCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTc                                                                                                                                                    <  4:271360/149‑1 (MQ=255)
   caacCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCgg                                                                                                                                                  >  1:152296/1‑149 (MQ=255)
    aacCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCggg                                                                                                                                                 >  2:178811/1‑149 (MQ=255)
    aacCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑ATGTTGCGGCAGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCggg                                                                                                                                                 >  2:21514/1‑149 (MQ=255)
      ccAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCggg                                                                                                                                                 <  2:23262/147‑1 (MQ=255)
      ccAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTg                                                                                                                                               >  2:15791/1‑149 (MQ=255)
      ccAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATCAAGGCCACCTTCGGGTg                                                                                                                                               >  4:28390/1‑149 (MQ=255)
      ccAGCTGTTGCGCAAGCTGATTGCTGAAAACGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCggg                                                                                                                                                 >  1:23262/1‑147 (MQ=255)
               gcgcAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCttttt                                                                                                                                       <  2:152296/148‑1 (MQ=255)
                         aTTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATgg                                                                                                                                  <  5:6705/143‑1 (MQ=255)
                         aTTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATgg                                                                                                                                  >  6:6705/1‑143 (MQ=255)
                         aTTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCCGAAGGTGGCCttaatgaa                                                                                                                                    <  5:43011/141‑7 (MQ=255)
                         aTTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCCGAAGGTGGCCttaatgaa                                                                                                                                    >  6:43011/1‑135 (MQ=255)
                             cTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTgg                                                                                                                        <  5:108030/149‑1 (MQ=255)
                              tGAAACCGATGCCT‑GGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGc                                                                                                                      >  2:186760/1‑149 (MQ=255)
                                aaaCCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGc                                                                                                                      <  2:173690/148‑1 (MQ=255)
                                      aTGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGgcc                                                                                                               >  7:233289/1‑149 (MQ=255)
                                          cTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGgccgcc                                                                                                            >  4:137033/1‑148 (MQ=255)
                                           tGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCaa                                                                                                          >  1:269974/1‑149 (MQ=255)
                                                   ggTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTc                                                                                                    >  5:20697/1‑147 (MQ=255)
                                                    gTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCCACTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTc                                                                                                    >  1:70553/1‑147 (MQ=255)
                                                    gTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGgcc                                                                                                               <  1:124207/135‑1 (MQ=255)
                                                    gTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCa                                                                                                  <  3:137033/148‑1 (MQ=255)
                                                    gTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAg                                                                                                 >  4:178910/1‑149 (MQ=255)
                                                        cTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCg                                                                                             <  2:70553/149‑1 (MQ=255)
                                                           tCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCa                                                                                          >  8:274351/1‑149 (MQ=255)
                                                              aaGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCCCCTTCGGGTGGCCGTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACg                                                                                        >  1:235451/1‑148 (MQ=255)
                                                                              aCCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCg                                                                           <  4:287358/145‑1 (MQ=255)
                                                                              aCCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCg                                                                           >  3:287358/1‑145 (MQ=255)
                                                                               ccGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAAt                                                                      <  3:178910/149‑1 (MQ=255)
                                                                                       gTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCAc                                                                                         <  4:102507/122‑1 (MQ=255)
                                                                                       gTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCAc                                                                                         >  3:102507/1‑122 (MQ=255)
                                                                                           aTGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCt                                                                                                                          <  1:146317/85‑1 (MQ=255)
                                                                                           aTGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCt                                                                                                                          >  2:146317/1‑85 (MQ=255)
                                                                                                   cTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTcgcg                                                    >  5:157170/1‑148 (MQ=255)
                                                                                                    tGTTGCGGCCGTGTAAGTTAGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGACGCAAACCGTTCCAGTGCGGCAAGCAGTTTGGGGTCGGTAATGACCTcgg                                                              >  3:206249/1‑137 (MQ=255)
                                                                                                     gTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAg                                                 >  5:150854/1‑149 (MQ=255)
                                                                                                        gCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGct                                              >  5:243170/1‑149 (MQ=255)
                                                                                                           gCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCTGTAATGCCCTCGGCGGTGTCGCGCAGGCtctc                                           <  6:243170/149‑1 (MQ=255)
                                                                                                             cGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCgg                                         <  5:210420/149‑1 (MQ=255)
                                                                                                                    gTTCGAACACTCTCTAGTTCATCAAGGCCACCTTCGGGTTGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCt                                  <  3:28390/149‑1 (MQ=255)
                                                                                                                        gAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGtt                               <  6:157170/148‑1 (MQ=255)
                                                                                                                          acacTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGGCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGa                            <  3:278432/149‑1 (MQ=255)
                                                                                                                              tctctAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACAGTACCAGTGCGGCAGGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCCAA‑AGc                       >  4:219651/1‑149 (MQ=255)
                                                                                                                                tctAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTc                     >  4:110118/1‑149 (MQ=255)
                                                                                                                                  tAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCTCTGTTCGAA‑AGCTc                     <  1:21514/147‑1 (MQ=255)
                                                                                                                                    gTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCcgcg                                    <  8:235791/131‑1 (MQ=255)
                                                                                                                                    gTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCcgcg                                    >  7:235791/1‑131 (MQ=255)
                                                                                                                                     ttCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGtt                               >  7:130038/1‑135 (MQ=255)
                                                                                                                                     ttCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGtt                               <  8:130038/135‑1 (MQ=255)
                                                                                                                                     ttCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCg                >  4:54436/1‑149 (MQ=255)
                                                                                                                                      taaaaaaGGCCACCCGAAGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAg                                                 <  1:79693/111‑1 (MQ=255)
                                                                                                                                      taaaaaaGGCCACCCGAAGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAg                                                 >  2:79693/6‑116 (MQ=255)
                                                                                                                                       cATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGc               <  7:274351/148‑1 (MQ=255)
                                                                                                                                       cATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCa              >  1:121823/1‑149 (MQ=255)
                                                                                                                                        aTTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCa              <  1:186760/148‑1 (MQ=255)
                                                                                                                                        aTTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCAt             <  2:269974/149‑1 (MQ=255)
                                                                                                                                        aTTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGAGCT‑TTCGAACAGCGCGGCCgaga             <  8:233289/149‑4 (MQ=255)
                                                                                                                                         ttAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCGAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCGTg            >  8:300053/1‑149 (MQ=255)
                                                                                                                                         ttAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTc                              >  4:140726/1‑132 (MQ=255)
                                                                                                                                         ttAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTc                              <  3:140726/132‑1 (MQ=255)
                                                                                                                                           aaGGCCACCCGAAGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGcc          >  6:67887/1‑149 (MQ=255)
                                                                                                                                            aGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGcc          <  6:20697/148‑1 (MQ=255)
                                                                                                                                             ggCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTCGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGATTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCAGTTCGAA‑AGCTCGGCCGCATGCCCg        >  6:158468/1‑149 (MQ=255)
                                                                                                                                                   cTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGc  >  6:82707/1‑149 (MQ=255)
                                                                                                                                                   cTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGGCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGc  >  6:262376/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAACAACCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGC  >  NC_002947/1531084‑1531378

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: ‑